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STUOIES  OF  COLLISIONAL  AND  NONLINEAR  RADIATIVE 
PROCESSES  FOR  DEVELOPMENT  OF 
UV  AND  XUV  SOURCES 

ABSTRACT 

Studies  of  the  interaction  of  matter  with  high  intensity  radiation  (>  101* 
W/cm2)  are  leading  to  the  observation  of  new  physical  phenomena  and  the  production 
of  new  classes  of  highly  excited  matter.  These  recently  discovered  processes 
involve  atoms,  molecules,  solids,  and  dense  plasmas.  A  central  factor  in  the 
discussion  is  the  ability  to  generate  very  high  levels  of  electronic  excitation  in  a 
manner  which  enables  the  system  to  remain  kinetically  cold  for  the  time  scale  of 
the  interaction.  This  conventionally  paradoxical  situation,  in  alliance  with  the 
existence  of  a  new  high-field  mode  of  channeled  propagation,  is  highly  conducive 
to  coherent  x-ray  generation  under  a  rather  wide  range  of  circumstances. 
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I.  INTRODUCTION 


The  controlled  concentration  of  power  in  materials  is  the  fundamental  issue 
for  the  creation  of  bright  sources  of  radiation  in  the  x-ray  range.  If  amplification 
is  desired,  this  issue  acquires  a  heightened  prominence,  since  gain  at  these  wavelengths 
(1-10  A)  requires  truly  prodigious  energy  deposition  rates1  spatially  organized  in  a 
high-aspect-ratio  volume  of  matter.  The  confluence  of  advances  in  femtosecond 
laser  technology,  the  results  of  a  range  of  physical  measurements,  and  theoretical 
analysis  all  suggest  that  electromagnetic  coupling,  in  the  strong-field  regime  (E>> 
e/a  I),  may  be  capable  of  establishing  the  demanding  conditions  necessary  for  such 
amplification.2  In  addition,  there  is  a  growing  body  of  evidence  that  fundamentally 
new  forms  of  highly  excited  matter  can  be  produced  with  strong-field  interactions. 
Therefore,  subjects  of  marked  significance  concerning  the  generation  of  radiation 
at  x-ray  wavelengths  with  strong-field  coupling  are  (1)  mechanisms  leading  to 

high  levels  of  electronic  excitation,  particularly  those  involving  ionization  and  the 
excitation  of  inner-shell  states,  (2)  the  characteristics  of  matter  which  govern  the 

efficiency  and  specific  power  of  x-ray  generation,  and  (3)  the  physical  processes 

determining  the  propagation  of  intense  waves  in  plasmas.  These  three  basic 
questions,  and  others  related  directly  to  them,  are  the  subjects  of  this  report. 

II.  GENERAL  DISCUSSION  OF  RESEARCH 
A.  Laser  Technology 

Advances  in  femtosecond  lasers  are  extending  the  exploration  of  multiphoton 
interactions  well  into  the  regime  for  which  the  external  field  is  greater  than  an 
atomic  unit  (e/a *).  The  performance  projected  for  ultraviolet  rare  gas  halogen 
technology  is  currently  being  realized2-11  and  new  near-infrared  solid  state  systems, 
such  as  Ti:Ala03,  are  under  vigorous  development.  Both  technologies12  should 

reach  a  field  strength  of  -  100  e/a|  with  instruments  that  produce  an  output 
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energy  of  ~  1  J.  With  these  experimental  tools,  intensities  approaching  I0ai 
W/cma  will  be  available,  a  value  comparable  to  that  produced  by  radiating  matter 
under  vigorous  thermonuclear  conditions. 

The  KrF*  (248)  nm  source  used  for  the  studies  discussed  herein  is  schematically 
shown  in  Fig.  (1).  This  instrument  has  an  output  pulse  width  of  =  600  fs,  a  typical 
output  pulse  energy  in  the  300-400  mj  range,  and  good  focusabiiity,9  nominally 
within  a  factor  of  two  of  the  diffraction  limit.  For  example,  with  a  simple  f/10 
CaF2  lens,  peak  intensities  of  ~  1017  W/cma  are  produced  even  though  such  a  lens 
exhibits  appreciable  spherical  aberration.  This  intensity  corresponds  to  a  peak 
electric  field  somewhat  above  one  atomic  unit.  The  application  of  more  sophisticated 
focusing  systems  can  clearly  generate  substantially  greater  field  strengths. 9,10 

A  recent  modification  has  incorporated  a  Ti:AlaO,  amplifier  into  the  system 
for  amplification  of  the  radiation  at  745  nm.  This  change  has  resulted  in  approximate¬ 
ly  a  40-fold  increase  in  the  248  nm  output  energy  from  the  crystals  to  a  value  of 

-  40  pJ.  This  enables  the  generation  of  a  significantly  more  intense  seed  beam  to 

drive  the  power  amplifier.  The  use  of  Ti:AlaO,  to  provide  radiation  at  745  nm 
greatly  enhances  the  performance  of  the  overall  laser  system.  In  addition, 

modifications  of  the  power  amplifier  will  be  made  this  spring  that  are  expected  to 
raise  the  output  energy  to  -  1  J  or,  perhaps,  somewhat  above. 

A  more  detailed  characterization  of  the  pulse  shape13  for  this  system  is 
becoming  available  with  the  use  of  a  technique  involving  a  third  order  intensity 
correlation14.  The  data  shown  in  Fig.  (2)  illustrate  the  character  of  the  pulse 
shape  observed  at  point  A  in  Fig.  (1).  One  of  the  properties  is  a  clearly  observed 
asymmetry  in  the  profile  of  the  pulse,  in  our  case  a  waveform  that  rises  more  rapidly 

than  it  falls.  The  pulse  width  measured  is  x  -  260  fs. 
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Fig.  (1):  Schematic  showing  the  configuration  of  the  ultrahigh-ln-intensity  KrF* 
laser  system  used  in  the  studies  discussed  in  the  text. 


4 


Measurement  of  248  nm  pulse  width  taken  at  point  (A)  In  Fig.  (1)  using 
4-wave  mixing  process,  (a)  Measured  third-order  autocorrelation  function 
of  260  fs  pulse  showing  clear  asymmetry,  (b)  Temporal  profile  of  pulse 
derived  from  (a). 


Intensity  (Arbitrary  Units) 


O'  ■  1  1  1  t  i  1  i.i 

0  6  12  18  24  30  36  42  48  54 

Radial  Coordinate  (pm) 


Fig.  (3):  Spatial  intensity  distribution  of  the  248-nm  seed  beam  using  an  f/7 
optical  system  measured  for  a  4  2  cm  diameter  beam.  The  diameter  at 
the  l/e  point  is  6.6  pm  a  value  which  is  ~  1  5-fold  the  theoretical  limit 
(4.2  pm).  The  dips  are  the  2  pm  grid  markers. 
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The  ability  to  concentrate  the  pulse  energy  in  a  small  focal  zone  is  another 
crucial  property  necessary  for  the  generation  of  high  field  strengths.  Fig.  (3) 
illustrates  the  focal  spot  measured  at  point  B  in  Fig.  (1)  with  an  f/7  focusing 
system.  The  measured  diameter  of  the  focal  zone  is  ~  6.6  pm,  a  value  which 
corresponds  to  -  1.5-fold  the  diffraction  limit. 

This  subpicosecond  KrF*  (248  nm)  laser  system  has  been  used  to  explore  the 
coupling  of  intense  radiation  to  atoms,  molecules,  and  solids.  The  findings  of  these 
experimental  studies  are  discussed  in  the  following  sections. 

B.  Studies  of  Ion  Production 

Studies  of  collision-free  ionization  have  been  a  valuable  source  of  information 
on  multiphoton  processes15-18.  In  particular,  such  experiments  give  a  direct  and 
unambiguous  measure  of  the  scale  of  the  energy  transfer  occurring  between  the 
radiation  field  and  the  free  atomic  or  molecular  target.  Although  most  studies 
have  been  conducted  with  atoms,  recent  work  has  begun  to  examine  molecules, 
including  certain  polyatomic  systems.  The  results  obtained  with  both  classes  of 
materials  are  discussed  below. 

1.  Atoms 

Studies  of  atomic  ion  production  have  examined  several  aspects  of  the  multiphoton 

19-22 

interaction  including  (1)  the  intensity  dependence.  (2)  the  atomic  number  (Z) 

dependence,2,17  (3)  the  frequency  dependence,18,23  and  (4)  the  dependence  upon 

pulse  width.21  Important  related  data  concerning  the  mechanism  of  ionization  has 

24-26 

also  come  from  numerous  photoelectron  studies  ,  work  which  has  accounted 

for  the  influence  of  the  optical  stark  effect.27-30 
21  22  31—37 

Several  models  ’  ’  have  been  used  to  interpret  the  extant  data  on 

collision-free  multiphoton  ionization  of  the  rare  gases.  If  the  customary  distinction 
founded  on  the  Keldysh  parameter31  is  made,  which  separates  the  multiphoton  (Y 
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>>  1)  and  tunneling  (Y  <<  1)  regimes  based  on  the  magnitudes  of  the  ionization 
and  ponderomotive  potentials,  theories  have  been  developed  for  both  the  former31- 
33  and  the  latter,21 -34_37  respectively. 

A  tunneling  picture  has  recently  been  applied  to  data  concerning  threshold 
ionization  of  the  rare  gases  in  the  infrared21  (1.06  pm)  and  ultraviolet22  (248  nm) 
spectral  regions.  For  the  infrared  study,  the  data  all  conformed  to  the  conventional 
tunneling  regime  (Y  <  1),  and  good  agreement  was  found  between  the  experimental 
results  and  the  theory.  However,  since  the  ponderomotive  potential  varies  as  the 
square  of  the  wavelength,  the  data  for  the  ultraviolet  experiments  corresponded  to 
conditions  significantly  outside  (Y  >  1)  of  the  validity  of  the  normal  tunneling 

analysis.  Interestingly,  although  the  ultraviolet  study22  involved  the  Keldysh  parameter 
in  the  range  1  <  Y  <  8,  reasonable  agreement  with  the  pure  tunneling  formulation 
was  still  found  for  the  threshold  intensities  occurring  in  the  1013  -  1018  W/cma  range. 

Figure  (4a)  presents  the  experimental  data  on  the  threshold  ionization  intensities 
of  the  rare  gases  observed  with  248-nm  irradiation.38  Clearly,  a  meaningful 
statement  of  a  threshold  intensity  requires  some  discussion  of  the  ionization  rate 
defining  the  limit  (threshold)  of  observation.  In  these  experiments,  the  ionization 
rate  at  threshold  was  estimated38  to  be  -  2  x  10*  s'*.  However,  for  all  the 
models  considered  below,  calculations  show  that  only  a  factor  of  two  change  in 
intensity,  either  way  around  this  threshold  value,  alters  the  ionization  rate  by  at 
least  two  orders  of  magnitude.  Thus,  it  is  not  necessary  to  know  the  threshold 
ionization  rate  to  high  accuracy  in  order  to  define  the  threshold  intensity  within 
relatively  narrow  limits. 

Two  important  features  in  Fig.  (4a)  are  readily  apparent.  Firstly,  there  is  no 
obvious  functional  dependence  of  the  threshold  ionization  intensity  1^  on  the 
ionization  potential  Ep,  and  second,  a  systematic  lowering  of  I^h  with  increasing 


atomic  number  exists.  These  two  salient  characteristics  have  been  basic  aspects 
of  the  data  on  muitiphoton  ionization  originating  with  some  of  the  earliest 
observations17  of  these  processes. 

The  simplest  model  of  tunneling  ionization  consists  of  a  one-dimensional  (1  — D) 
coulomb  potential.21  In  an  external  static  electric  field,  an  electron  in  this  1-D 
atom  sees  a  finite  potential  barrier  to  ionization  whose  width  and  height  depend 

on  the  strength  of  the  field.  The  threshold  intensity  It^  for  ionization  can  then 

be  defined  as  the  intensity  at  which  the  potential  barrier  is  reduced  to  the  ionization 
potential  of  the  atom,  the  situation  allowing  for  classically  allowed  escape  of  the 

electron.  The  result21  of  this  picture  is  the  relationship 

Ith  =  cEp/ 1 28ire®Z2,  (1) 

where  Ep,  2,  c,  and  e  are  the  ionization  potential,  the  charge  of  the  resulting 
ion,  the  speed  of  light,  and  the  charge  of  an  electron,  respectively. 

This  model  has  two  inherent  defects  which  tend  to  cancel  each  other  resulting 
in  a  fairly  accurate  prediction  of  It^.  First,  when  the  barrier  is  lowered  to  the 
ionization  potential,  the  barrier  is  completely  removed  producing  ionization  rates 
characteristic  of  atomic  time  scales  (~101‘  s'1).  Since  a  rate  of  this  magnitude 

is  much  higher  than  the  applicable  threshold  rate,  the  resulting  value  of  the 
threshold  ionization  intensity  tends  to  be  overestimated.  However,  in  contrast  to 

the  1-D  picture,  the  lowering  of  the  potential  barrier  to  the  ionization  potential 
for  a  three  dimensional  (3D)  atom  at  the  intensity  given  by  Eq.(l),  will  occur  in 

only  a  single  spatial  direction.  In  all  other  directions  the  potential  barrier  will  be 

higher,  a  situation  leading  to  an  underestimation  of  the  threshold  ionization  intensity. 

Eq.(l)  reveals  two  important  scaling  relationships  despite  these  corrections. 
It  shows  that  besides  a  dependence  on  Ep,  1^  also  depends  inversely  on  Z2.  This 
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latter  scaling  factor  can  be  accounted  for  by  simply  considering  the  combined  variable 

Ith*Z2.  In  addition,  the  variable  Ith*Z2  varies  in  direct  proportion  to  E^. 

Clearly,  a  plot  of  l^Z2  versus  Ep  should  exhibit  these  basic  scaling  relationships 

and  Fig.  (4b)  shows  such  a  representation  for  the  data  illustrated  in  Fig.  (4a).  The 

experimental  data  are  seen  to  fall  almost  perfectly  on  a  straight  line  with  a  slope 

of  four  and  the  systematic  dependence  on  atomic  number  has  been  significantly 

reduced.  Furthermore,  this  agreement  for  the  threshold  intensities  occurs  over  a 

span  of  five  orders  of  magnitude  in  the  parameter  Ith-Z2.  Nevertheless,  the 

values  of  Eq. (1 )  appear  consistently  too  high  and  there  remains  a  more  complex 

detailed  dependence  on  the  atomic  number. 

More  sophisticated  formulations  of  this  process  have  been  developed.  For 

example,  an  improved  model  involves  a  3-D  atom  and  the  calculation  of  the  tunneling 

34  35 

rate  through  the  potential  barrier  at  an  arbitrary  intensity.  ’  The  resulting 
ionization  rate  for  a  static  electric  field  in  this  theory39  is  given  by 

Wst(Es)  =  4alo(2Ep)5/2Es~1exp(  -2/3(2Ep)3/2Es_1  ).  (2) 

for  (2Ep)3//2  >  >  Es,  where  Es  is  the  static  field  in  atomic  units  (e/a2),  Ep  the 

ionization  potential  in  atomic  units  (e2/a0),  and  w0  the  atomic  frequency  in  atomic 

units  (4.1  x  I0ie  s'1).  In  order  to  compute  the  ionization  rate  generated,  Eq.(2) 

can  be  time-averaged  over  one  cycle  of  the  field.  This  time  average  can  be  done 
exactly  in  terms  of  the  K0  modified  Bessel  function.40  However,  in  the  limit 

O  / O  Og 

(2Ep)  >>ES,  the  expression  for  the  average  reduces  to  a  simpler  form  given  by 

Wac  =  (3/ir)1/2E1/2(2Epr3/4Wst(E),  (3) 
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where  £  is  the  peak  of  the  alternating  field  in  atomic  units. 

Eq.(3)  can  be  further  improved  to  include  non-hydrogenic  systems.  The  static 

36  37 

ionization  rate  in  this  orbital  picture  is  given  by 
(2i  +  1)(i,+  I  m  | )! 

wst-  Cjk  Ep  — - — - 

2|m|(|m|)l(l-|m|)! 

-1  /2 

where  n*  is  the  effective  principle  quantum  number  [n*=Z(2Ep)  '  ],  L  the  orbital 

angular  quantum  number,  and  m  the  magnetic  quantum  number.  Cn*.  is  a  numerical 
constant  on  the  order  of  two.  The  approximate  time  average  of  Eq.(4)  is  also 
given  by  Eq.(3). 

The  last  treatment  that  we  consider  is  the  well  known  Keldysh  theory.^1  Since 
the  Y  parameter  is  in  the  intermediate  regime,  1  <  Y  <  8,  a  limiting  form41  of  the 
theory  is  not  valid  and  the  full  form  is  necessary. 

Figures  (5a)  through  (5c)  compare  the  experimental  data  for  He,  Ar,  and  Xe 
from  Fig.  (5)  with  the  results  of  the  four  models  represented  by  Eq.(1),  Eq.(3) 
using  the  static  ionization  rates  for  the  hydrogenic  and  nonhydrogenic  systems, 
and  the  Keldysh  theory,  respectively.  Several  features  are  manifest.  First,  both 
the  3-D  simple  atom  picture  and  the  Keldysh  theory  agree  well  with  the  experimental 
results  for  helium.  However,  they  become  consistently  worse  for  the  heavier 
atoms  particularly  at  the  higher  charge  states.  This  trend  makes  sense  as  these  models 
assume  exact  hydrogenic  potentials,  a  poor  approximation  for  the  heavier  atoms. 
Second,  as  noted  above,  the  1-D  model  uniformly  gives  too  high  a  threshold 
intensity  and  tends  to  be  least  successful  in  describing  the  neutral  species.  Third, 
the  3-D  complex  model,  expressed  by  Eq.(4),  consistently  gives  reasonable  agreement 
over  the  entire  range  of  atomic  number  and  charge  states  studied.  Fourth,  the  single 


^2(2Ep)3/2E  1|f1 


-|m|-1 


exp(-2/3(2Ep)3/2E  \ 
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Fig  (4):  (a)  Threshold  ionization  intensities  of  the  rare  gases  with  subpicosecond 

24d-nm  irradiation,  (b)  Same  data  replotted  with  Eq.  (1)  to  show  the 

scaling  relationships.  All  intensities  are  peak  values.  [See  Fig.  (5)  for 

the  uncertainty  in  the  intensities). 
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Figure  5:  Threshed  Ionization  intensities  for  (a)  helium,  (b)  argon,  and  (c)  xenon 
from  Fig.  (4)  along  with  calculated  results  from  various  theories.  The 
i-O  Modei  corresponds  to  Eq.  (i),  the  3-0  Simple  Atom  to  Eqs.  (2)  and 
(3),  the  3-0  Complex  Atom  corresponds  to  Eqs.  (3)  and  {♦).  The  Keldysh 
line  was  calculated  from  the  equations  In  Ref.  (19).  All  intensities  are 
peak  values.  The  relative  uncertainty  in  the  intensities  is  indicated  by 
an  error  bar  in  the  legend. 
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point  which  deviates  considerably  from  all  models  is  the  threshold  ionization  of 
neutral  xenon  shown  in  Fig. (5c).  However,  this  anomaly  clearly  arises  from  resonant 
multiphoton  ionization,  since  there  is  a  very  close  two-photon  resonance  (5p  -  6p) 
with  248-nm  radiation. 

The  three  tunneling  models  considered  here  are  all  quasi-static  theories, 
equivalently,  Y  is  vanishingly  small.  However,  as  mentioned  above,  the  ultraviolet 
data  all  correspond  to  a  regime  in  which  Y  is  significantly  greater  than  one. 
Although  a  full  tunneling  ionization  model  for  a  coulomb  potential  and  an  arbitrary 
Y  has  been  solved,36  the  form  of  the  solution  is  extraordinarily  complicated.  The 
consideration  of  a  6-function  potential  results  in  a  considerably  simpler  expression 
for  the  ionization  rate36  and,  from  this  model,  the  effect  of  a  variation  in  Y  on 
the  threshold  intensities  has  been  estimated.  The  finding  is  that,  up  to  a  value 

of  Y  =  3,  the  threshold  intensities  decrease  by  less  than  a  factor  of  two.  Therefore, 
the  dependence  is  not  particularly  strong,  but  the  tendency  is  for  a  decrease  in 
the  threshold  intensity. 

The  experimental  measurement  of  threshold  ionization  intensities  is  a  relatively 
simple  procedure,  but  gives  results  which  are  rather  insensitive  to  the  actual 
ionization  rates.  Theoretical  threshold  values  are  correspondingly  insensitive  to 
the  details  of  a  given  model  on  the  same  basis.  It  follows  that  it  is  difficult  to 
use  threshold  data  to  evaluate  the  merits  of  various  theories  of  ionization,  unless 

very  accurate  measurements  are  made.  On  the  other  hand,  the  measurement  of 

threshold  intensities  provides  a  practical  and  easy  method  of  determining  focused 

intensities  to  a  reasonable  accuracy. 

Three  main  conclusions  have  emerged  from  studies  of  this  nature  on  the 
production  of  atomic  ions.  First,  even  using  248-nm  irradiation,  threshold  ionization 
intensities  can  be  rather  accurately  described  by  a  relatively  elementary  model  of 
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tunneling  ionization,  at  least  for  intensities  i  101*  W/cm2.  Second,  while  a  1-0 
coulomb  potential  model  qualitatively  gives  the  basic  scaling  laws  for  tunneling 
ionization,  a  3-D  complex  atom  picture  is  required  to  give  uniformly  better  quantitative 
agreement  over  the  entire  range  of  materials  and  charge  states  studied.  Third, 

the  measurement  of  threshold  ionization  intensities  is  rather  insensitive  to  the 
details  of  the  theoretical  models  and,  therefore,  has  limited  usefulness  in  testing 
alternative  models  of  ionization.  Finally,  on  the  basis  of  the  scaling  revealed  in 
these  studies,  we  note  that  with  the  use  of  intensities  in  excess  of  101*  W/cm2, 
which  are  now  achievable,  it  should  be  possible  to  remove  completely  the  4d  shell 
from  xenon  producing  krypton-like  Xe11*.  The  total  energy  invested  in  ionization 
for  this  case  would  be  -  3.5  keV.  Species  with  this  level  of  ionization  are  appropriate 
for  the  production  of  radiation  in  the  soft  x-ray  range.  Furthermore,  considering 
that  these  ions  can  be  kinetically  cold,  this  is  a  rather  unusual  state  of  highly 

ionized  matter.42 

2.  Molecules 

Studies  of  molecular  ionization  have  produced  additional  insights  into  the 

43—51 

dynamics  of  multiphoton  processes.  Energetic  coulomb  explosions  have  been 

studied  having  kinetic  energies  of  the  fragments  with  several  tens  of  eV  per 

particle,  evidence  for  molecular  inner-shell  excitation  has  been  observed,  and 
molecular  atomic  site-dependent  energy  deposition  has  been  detected.  Furthermore, 
these  measurements  have  given  rise  to  the  hypothesis  that  appropriately  designed 
molecules  may  enable  the  selective  production  of  highly  excited  ionic  fragments 

suitable  for  the  generation  of  radiation  in  the  x-ray  range52.  This  possibility  is 
discussed  in  Section  III. A  below. 
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a.  Atomic  Ion  Production 


Many  interesting  results  have  come  from  the  study  of  a  simple  triatomic 
system  (NaO)  whose  structure  is  linear  and  asymmetric  (C^y).  With  irradiation  at 
248  nm  at  an  intensity  of  ~  3  x  101*  W/cm*,  both  atomic  and  molecular  (diatomic) 
ionic  fragments  are  seen.  Furthermore,  since  some  of  the  studies52  have  been 
performed  with  isotopic  material  (lHNl*Nl*0),  it  has  been  possible  to  distinguish 
the  original  molecular  sites  of  all  the  fragments  observed  in  simple  measurements 
of  the  ions. 

The  study  of  energetic  atomic  fragments  from  NaO  can  reveal  considerable 
information  on  the  process  of  molecular  multiphoton  ionization.  In  particular,  it 
is  possible  to  measure  the  kinetic  energies  of  the  products  arising  from  a  molecular 
coulomb  explosion  with  a  suitably  arranged  time-of-flight  apparatus.47  As  an 
example,  we  discuss  below  the  results  of  recent  measurements  of  atomic  fragments 
produced  in  ,'*Nal*0  conducted  under  conditions  for  which  the  fragment  energy 
distributions  have  been  determined. 

Figure  (6)  illustrates  the  observed  energy  distribution  for  ‘“N3*  produced 
from  Of  course,  this  experiment  does  not  distinguish  between  the  two 

nonequivalent  nitrogen  positions  and  both  nitrogen  sites  in  the  molecule  contribute 
to  the  signal.  Two  aspects  of  the  data  are  manifest,  namely,  the  presence  of 
significant  structure  in  the  distribution  and  the  scale  of  the  maximum  energy 
observed.  In  consideration  of  the  latter,  the  arrow  in  Fig.  (6)  represents  a 
kinetic  energy  of  ~  77.5  eV,  the  value  which  represents  the  endpoint  of  the 
observable  signal.  Interestingly,  77.5  eV  is  also  the  ionization  potential53  of  ground 
state  N3*  ions.  This  may  be  a  coincidence,  but  a  possible  implication  of  these 
data  is  that  a  mechanism  exists  which  effectively  limits  the  kinetic  energy  of  the 
fragment  to  its  corresponding  ionization  potential.  Molecular  crossings  to  other  state 
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Fig.  (6).  The  kinetic  energy  distribution  of  l*N**  fragments  observed  in  preliminary 
studies  from  irradiated  at  248  nm  with  ~  600  fs  pulses  at  a 

peak  intensity  of  -  3  x  10“  w/cma.  The  location  of  the  arrow  (♦) 
corresponds  to  a  kinetic  energy  for  the  fragment  of  -  77.5  eV,  the 
value  of  the  ionization  potential  of  N**. 
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involving  differently  charged  species,  in  this  case  the  N1**  +  e"  channel,  could 
provide  the  necessary  mechanism.  In  this  way,  sufficiently  energetic  molecular 
dissociation  may,  through  the  action  of  a  crossing,  be  directed  into  channels 
representing  either  electronic  excitation  or  ionization.  This  process  could  be 
imagined  as  one  involving  electron  promotion  in  a  manner  similar  to  that  developed 
by  Fano  and  Lichten54  for  the  description  of  energetic  ion-atom  collisions. 
Basically,  new  channels  including  NH*  could  be  opening  with  the  consequence  that 
the  N3*  yield  is  depressed.  Related  aspects  of  this  type  of  interaction  are  considered 
in  Section  III. A. 

The  N3*  kinetic  energy  distributions,  as  measured  from  both  N2  and  N20, 
are  distinguished  in  another  way,  which  may  be  related  to  this  observation.  From 
both  N2  and  N20  nitrogen  ions  with  a  charge  as  high  as  Ns*  are  detected.55 
Among  these  ions,  however,  the  N3*  species  clearly  are  clearly  the  most  energetic. 
This  aspect  of  the  observed  behavior  remains  to  be  understood. 

A  rough  appraisal  of  the  kinetics  of  molecular  dissociation  for  the  N3*  data 
shown  in  Fig.  (6)  is  also  informative.  For  this  estimate  we  will  assume  that  the 
dissociation  occurs  as  a  two-body  event  of  the  form 

(N20)(3  +  q)+  - ►  N3+  +  (NO)9*  (5) 

with  the  participation  of  a  (NO)^*  fragment.  In  the  ground  state56,  N20  has 
equilibrium  N-N  and  N-0  spacings  of  1.128  A  and  1.184  A,  respectively.  If  we 
assume  that  the  effective  location  of  the  (q  +  )  charge  on  the  NO  fragment  is  located 
at  its  midpoint,  then  the  N3*  ion  is  at  a  distance  of  r0  -  3.25  ae.  The  total  coulomb 
energy  release  Ec  is  then  given  approximately  by 
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(6) 


The  full  kinetic  energy  Ek  associated  with  the  dissociation,  given  the  equal  magnitudes 
of  the  momenta  for  the  two  fragments,  is 


in  which  Efsj3*  corresponds  to  the  measured  N3*  kinetic  energy  with  MN  and  MNq 
representing  the  masses  of  the  N3*  and  (NO)0!*  fragments,  respectively.  For  EN3* 
■  77.5  eV,  E|<  =  114  eV.  Equating  E^  =  Ec,  we  obtain  q  =  4.6,  which  we  interpret 
as  q  =  5  for  the  NO1**  fragment.  Since  the  NOq*  fragment  involves  two  atoms,  an 
effective  value  of  q  =  5  for  the  diatomic  fragment  does  not  appear  unreasonable 
for  conditions  that  would  generate  N3*.  Obviously,  considerable  further  study  will 
be  necessary  to  fully  ascertain  the  dynamics  of  these  complex  dissociations. 

An  important  basic  consideration  relates  to  these  molecular  data.  In  the 

example  examined  in  relation  to  Fig.  (6),  the  nascent  system  produced  is  a  species 

of  the  form  (NaO)qo  with  a  charge  q„  representing  several  electrons,  perhaps, 

with  a  value  of  q0  -  8  as  estimated  above.  Moreover,  this  multiply  charged 

system  is  produced  by  an  interaction  causing  a  negligible  transfer  of  momentum  to 

the  molecule.  Core  state  excitation  followed  by  Auger  decay  can  produce  multiple 

ionization,  but  generally  that  process  generates  for  this  molecule  only  doubly 
57  58 

ionized  material  ’  with  perhaps,  a  small  fraction  of  triply  ionized  ions.  Basically, 
the  maximum  level  of  ionization  that  can  be  attained  by  the  Auger  mechanism  for 
molecules  composed  of  atoms  having  relatively  low  atomic  numbers  is  modest. 
Energetic  heavy  ion  collisions59,60  can  generate  considerably  higher  levels  of 
ionization,  but  this  interaction  generally  imparts  a  very  substantial  recoil  momentum 
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to  the  constituents  of  the  molecular  targets®2.  However,  the  combination  of  (1)  a 
recoil-free  interaction  and  a  (2)  high  level  of  ionization  can  be  quite  simply  achieved 
by  the  strong-field  multiquantum  mechanism.  This  method,  therefore,  enables  the 
unique  production  of  forms  of  kinetically  cold  highly  excited  matter42  that  are 
essentially  impossible  to  produce  by  other  available  methods®2.  Therefore,  these 
results  give  evidence  that  entirely  new  forms  of  excited  material  can  be  generated 
by  strong-field  coupling.  Further  discussion  of  this  important  point  is  presented 
in  Section  III.C.2. 

b.  Molecular  Site-Specific  Energy  Deposition 

42  47  51 

Previous  experimental  studies  ’  of  N}  and  NaO  suggest  that  strong- 
field  coupling  may  produce  atomic  site-selective  energy  deposition  in  molecules. 
Such  a  process  would  represent  the  multiphoton  analogue  of  the  site-specific  excitation 

6 1  64 

that  can  be  achieved  with  the  absorption  of  x-rays.  This  could  be  observed, 

for  example,  with  the  nitrogen  atoms  in  NaO,  a  linear  molecule  with  the  structure 
NNO.  The  two  nitrogen  atoms  sites  in  NsO  are  nonequivalent;  one  resides  on  the 
end  of  the  molecule  while  the  other  is  located  in  the  central  position. 

In  preliminary  experiments,  the  comparison  of  the  fluorescence  spectra 
exhibited  by  Na  and  NaO  in  the  70-nm  to  80-nm  range,  in  relation  to  the  behavior 
of  certain  NJ*  lines  involving  lsa2p3  doubly  excited  levels,  has  furnished  specific 
evidence  indicating  this  type  of  site-selective  behavior.  In  order  to  verify  this 
hypothesis,  isotopic  studies  of  ionic  fragmentation  have  been  performed  to  more 
fully  examine  these  processes.  For  example,  since  lsN,l*NO  and  1HN,sNO  can  be 
readily  synthesized®®  from  isotopically  substituted  NHhNO,,  the  phenomenon  of 
site-selective  energy  deposition  can  be  sensitively  and  easily  detected  by  comparative 
studies  of  ion  production  in  l4Nll,N0  and  the  isotopic  variants.  The  same  isotopic 
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material  was  also  used  to  facilitate  the  measurements  described  below  in  Section 
II.B.2.C  concerning  the  production  of  diatomic  ionic  fragments. 

Initial  experiments  using  l,,NlsNl*0,  have  confirmed  the  expected  site 
dependence  from  measurements  of  the  charge  state  spectrum  of  the  nitrogen  ion 
fragments.  A  clear  demonstration  of  this  dependence  is  illustrated  in  Fig.  (7). 
For  the  range  of  flight  times  involving  the  singly  charged  species  N*  and  0*  (13 
<  M/q  ^  17),  normal  (ll*Nl,*Nl*0)  material  produces  the  signal  shown  in  Fig.  (7a). 

In  this  case,  it  is  impossible  to  distinguish  the  N’  ions  coming  from  the  two 
different  molecular  sites. 

In  considering  the  data  shown  in  Fig.  (7),  it  should  be  understood  that  the 
shapes  of  the  ion  signals  are  determined  by  both  the  instrumental  response  and 
the  kinetic  energies  of  the  fragments.  The  individual  ion  signals  are  actually 
composed  of  two  components  which  are  not  fully  resolved  in  these  relatively  low 
energy  resolution  data.  One  component  arises  from  ions  initially  formed  with  velocities 
parallel  to  the  time-of-flight  direction  while  the  other  arises  from  ions  whose  initial 
velocity  is  anti-parallel  to  that  vector.  The  two  peaks  are  of  unequal  height  and 
width  due  to  the  detailed  nature  of  the  trajectories  which,  in  turn,  are  determined 
by  the  extraction  and  acceleration  voltages  used  in  the  time-of-flight  apparatus47. 
Naturally,  the  total  ion  yield  is  given  by  the  integral  of  the  signal. 

The  same  region  of  the  time-of-flight  spectrum  is  shown  in  Fig.  (7b)  for 
l>*NlsNl*0.  In  this  case,  a  small  quantity  of  H0Ar  was  present  and  the  '40Ar3*  signal 
is  clearly  seen.  The  argon  peak  can  be  used  to  calibrate  t.ie  time-of-flight  axis 
and  give  a  direct  measure  of  the  instrumental  width  determined  by  the  apparatus. 

In  this  regard,  note  the  narrow  and  symmetric  shape  of  the  H0Ar3*  signal. 

Three  distinct  signals  arising  from  the  lHNlsNl,0  are  clearly  represented  in 
Fig.  (7b).  Plainly,  the  total  yield  of  the  lsN*  iors  is  approximately  one  half  of 
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(a)  lW*0 


(b) 


l"Nl*NuO 


Fig.  (7)  Charge  state  fragments  observed  from  N.O  in  the  range  13  < 
(a)  Spectrum  from  ‘-N^hT'O.  (b)  Spectrum  from  ‘"Nl*NuO.  See  text  for 


M/q  «  17 
discussion 
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the  l"N*  yield.  Although  several  different  detailed  pathways  can  lead  to  the  observed 

distribution,  the  outcome  is  that  the  total  probability  of  producing  singly  charged 

ions  depends  upon  the  molecular  site.  We  note  that  the  widths  of  l5N  and  lHN 

signals,  which  indicate  the  range  of  kinetic  energies  of  the  fragments,  are  also 

quite  different,  with  the  former  exhibiting  a  considerably  narrower  value.  This 

conforms  to  the  expectation  that  the  atom  in  the  center  (1SN)  would  develop,  in 

the  process  of  fragmentation,  a  relatively  smaller  kinetic  energy  than  the  atoms 

50  66 

on  the  ends  of  the  molecule.  ’  Similar  general  characteristics,  in  regard  to 
both  the  signal  strengths  and  widths,  are  also  seen  in  the  N3*  channel  in  comparison 
of  the  ^N3*  and  1SN3*  signals. 

c.  Diatomic  Ion  Yields 

Informative  behavior  is  also  exhibited  by  the  observation  of  diatomic  fragments 

such  as  Na  and  NO*.  The  data  on  these  ions  for  lHNlBNl*0  are  shown  in  Fig. 

(8).  In  this  case,  the  lsNlsO*  yield  appears  to  be  slightly  greater  than  the 

l,*NlsN*  yield.  This  result  provides  an  interesting  comparison  with  the  yield  ratio 

of  these  fragments  formed  by  other  known  processes.  For  example,  it  is  considerably 

67  68 

different  than  that  observed  by  photoionization  ’  of  NaO  at  30.4  nm  and  also 

appears  to  differ  from  the  ratio  observed  in  other  multiphoton  studies,®®  conducted 

under  somewhat  different  conditions  of  irradiation,  at  least  for  the  NaOa*  — ►  N* 

+  NO*  channel.  However,  the  observed  ratio  is  rather  close  to  that  seen  when 

monochromatic  soft  x-rays  (-  405  eV)  are  used  to  excite  NaO.  With  a  quantum 

energy  of  ~  405  eV,  the  2a  core  electron  on  the  central  nitrogen  atom  is  selectively 

excited  to  an  empty  3ir  orbital,  thereby  initially  concentrating  the  deposited 

57 

energy  at  the  molecular  midpoint.  A  similar  conclusion  follows  when  the  is  — * 

2p*  orbital  is  excited  on  the  O  atom  at  ~  536  eV.  Interestingly,  it  is  also  known 

57  64 

that  Auger  spectra  in  NaO  are  sensitive  to  the  site  of  excitation  Finally, 
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the  comparison  of  the  results  in  Fig.  (7)  with  those  obtained  from  electron  scattering70 

is  at  present  uncertain,  since  the  latter  reveal  an  unexpected  difference  in  the 

57 

NO*/N*  ratio  in  relation  to  the  findings  observed  with  soft  x-ray  excitation 

Overall,  the  comparison  on  the  basis  of  the  fragment  yield  ratio  of  the  multiphoton 

data  with  the  results  of  other  experiments  on  NaO  suggests  that  the  experimental 

findings  are  consistent  with  a  picture  of  the  interaction  involving  a  significant 

level  of  localization  of  the  deposited  energy  in  the  NO  region  of  the  molecule. 

Considerable  further  work  will  be  necessary  to  fully  test  this  hypothesis. 

The  appreciable  yield  of  NO*  gives  rise  to  a  further  interpretation71.  The 

NO*  system  represents  a  triple  bond  ,  whereas,  initially  in  the  NaO  system,  the 

considerably  stronger  bond  occurs  between  the  two  nitrogens.  This  implies  the 

movement  of  a  significant  electron  density  along  the  axis  of  the  molecule  in  the 

mechanism  leading  to  the  formation  of  the  NO*  fragment.  Given  the  large  number 

72 

of  molecular  state  crossings  that  are  expected  to  exist  at  the  level  of  excitation 

pertinent  to  such  a  transition,  it  seems  unlikely  that  this  flow  of  charge  would 

occur  adiabatically  and  channels  leading  to  the  production  of  excited  states  may 

be  significant.  The  application  of  a  sufficiently  strong  external  field  could  be 

effective  in  generating  such  a  motion.  Considerable  further  analysis  of  these 

questions,  particularly  in  light  of  known  information  on  the  photodissociation  of 

72  73 

highly  excited  molecules,  '  is  clearly  necessary. 

C.  Radiation  Studies 

Over  the  past  several  years,  many  processes  involving  the  production  of 

radiation  have  been  examined.  These  studies  have  involved  measurements  of 

51  74  75  74  76  77  7fi 

fluorescence,  ’  ’  harmonic  generation  in  gases,  ’  '  parametric  processes, 

and  x-ray  generation  from  solid  matter.70-0^  In  the  sections  following  directly 

below,  we  discuss  certain  aspects  of  this  area  of  research. 
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10.233  10.732  11-231 

Flight  Time  (ysec) 


Fig.  (8):  Charge  stata  spectrum  of  diatomic  ions  (1>,N1,N)*  and  (l*Nv*0)’  observed 
from  l*Nl#Nl,0.  For  discussion  see  text. 
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1.  Fluorescence  from  Atomic  Ions 

The  experimental  arrangement  used  to  conduct  the  studies  of  the  radiation 
produced  in  gases  is  illustrated  in  Fig.  (9).  The  source  initially  used  for  irradiation 
of  the  target  gases  was  KrF*  (248-nm)  laser  system7  that  produced  pulses  having 
a  maximum  energy  of  -  20  mJ  and  a  minimum  pulse  duration86  of  ~  350  fs.  A  23- 
cm  (f/10)  focal-length  lens  focused  the  248-nm  laser  radiation  into  a  gaseous  target 
produced  by  a  modified  Lasertechnics  pulsed  gas  jet.  With  this  optical  system, 
the  intensity  in  the  focal  region  was  estimated  to  be  in  the  range  of  ~  10IS-  I0ie 
W/cm2.  Subsequently  experiments  utilizing  the  more  powerful  source  illustrated  in 
Fig.  (1)  enabled  the  achievement  of  focal  intensities  approximately  one  order  of 
magnitude  higher. 

The  gas  jet  was  mounted  11.5  cm  in  front  of  the  entrance  slit  of  a  2.2-m 
grazing-incidence  spectrometer  (McPherson  Model  247)  equipped  with  a  600-line/mm 
gold-coated  spherical  grating  blazed  at  120  nm.  A  single-stage  microchannel  plate 
with  a  phosphored  fiber-optic  anode  served  as  the  detector.  Mounting  of  the 
detector  tangentially  to  the  Rowland  circle  of  the  spectrometer  allowed  the  observation 
of  radiation  between  7.5  and  80  nm.  The  FWHM  resolution  and  accuracy  of  the 
spectrometer-detector  system  were  typically  -  0.1  nm.  The  pulsed-gas  valve,  which 
was  modified  to  allow  a  backing  pressure  up  to  ~  5  x  10”  Torr  (1000  psi),  was 
typically  operated  in  the  range  of  ~  10*  to  -  3  x  10“  Torr  with  a  pulse  repetition 
rate  of  2  Hz.  Background  pressures  in  various  parts  of  the  apparatus  are  indicated 
in  Fig.  (9).  The  laser  was  focused  to  a  position  a  few  hundred  micrometers  above 
the  nozzle  tip  which  had  a  diameter  of  0.5  mm  and  a  throat  depth  of  1  mm.  The 
estimated  gas  density  in  the  interaction  region  was  ~  101*  cm'3. 

In  the  studies  employing  the  gas  jet,  fluorescence  shorter  than  10  nm  from 
orbitally  excited  ions  of  Ar,  Kr,  and  Xe  has  been  observed.  A  region  of  the  spectrum 
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Radiation  Detection  System 


Fig.  (9):  Schematic  of  pulsed-gas  jet  and  spectrometer-detector  assembly  used  in 
studies  of  emitted  radiation.  Typical  background  pressures  are  indicated. 
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obtained  for  Kr  is  shown  in  Fig.  (10),  which  illustrates  two  general  features  of 
these  observed  spectra.  They  are  (1)  identified  radiation  from  high  charge  states 
and  (2)  an  abundance  of  unidentified  lines. 

Table  I  summarizes  the  main  properties  of  the  observed  fluorescence  in  Ar,  Kr, 
and  Xe.  In  all  three  materials  a  large  number  of  emissions  was  seen.  Although 
specific  identification  was  often  possible,  as  illustrated  in  Fig.  (10)  and  Table  I,  a 

considerable  fraction  of  the  observed  transitions  could  not  be  associated  with  the 
known  spectra  of  any  ionic  species.  It  seems  possible  that  multiply  excited  levels 
and/or  core-excited  states  could  account  for  some  of  these  anomalous  spectral 
features.  At  present,  the  information  on  the  properties  of  such  states,  particularly 
for  ionic  spectra,  is  quite  meager. 

The  large  aperture  KrF  (248  mm)  laser  system9  shown  in  Fig.  (1)  has  also 

been  used  for  spectroscopic  studies  of  gaseous  targets.  With  this  system,  the 

intensity  produced  is  somewhat  greater  than  that  used  for  the  lines  observed  in 

Table  I  and  radiation  from  the  L-shell  of  argon  originating  from  Ar**  has  been 

observed.  The  spectrum  of  argon  shown  in  Fig.  (11)  clearly  shows  the  presence 

of  the  well  known53  2s2p*  — ►  2ss2p*  doublet  along  with  the  4p  - ►  3s  line  of 

Ar7*.  It  is  also  noted  that  the  production  of  ground  state  Ar9*  from  Ar9*,  by 

the  removal  of  a  2p  electron,  requires  a  minimum  of  422  eV  and  that  the  total 

energy  investment  needed  to  produce  Ar9*  from  the  neutral  atom  exceeds  a 

1 08 

kilovolt.  Analysis  of  these  data  showed  that  the  most  consistent  explanation  is 
that  the  laser  pulse  creates  a  plasma,  through  a  process  not  yet  thoroughly  understood, 
containing  Ar7*,  Ar**,  and  Ar9*  with  an  effective  electron  temperature  of  -  20 
eV.  This  temperature  is  sufficient  to  excite  the  Ar7*  states  and  the  first  excited 
state  of  Ar9*,  but  no  levels  in  Ar**.  The  conclusion  is  that  the  lines  observed 
in  Ar7*  are  fully  consistent  with  a  collisional-radiative  model  of  the  plasma  with 
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■  Kr7  + 


4Pj/a  -  5Sl/a  0  Unidentified 

29.71  nm 


Wavelength  (nm) 


Fig.  (10):  A  section  of  the  Kr  spectrum  near  30  nm  illustrates  identified  features 

from  Kr7*  and  unidentified  lines  of  comparable  strength. 
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Intensity  (arb.  units) 


Fig.  (ii):  Spectrum  of  Ar  showing  features  from  Ar7*  and  Ar**  observed  in 

a  gaseous  target  with  a  subpicosecond  248  nm  pulse  produced  by 
the  I0a  cm3  aperture  laser  system  illustrated  in  Fig.  (1). 
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Table  I 

Main  Properties  of  the  Fluorescence  Observed 


from 

Ar,  Kr,  and  Xe 
Intensity  of 

with  248-nm  Radiation 
-  io“  -  10“  W/cma 

at  an 

Material 

Ar 

Kr 

Xe 

Ionic  charge 
states  observed 

4  +  ,  5  +  , 

6  +  ,  7  + 

6  +  ,  7+, 

8  + 

6  +  .  7+, 

8  + 

Wavelength  range 
of  emission  (nm) 

12.0-47.1 

9.9-45.3 

9.8-48.1 

Typical 

identified 

transition 

7  + 

4f  — ►  3d 

7  + 

5p  — ►  4s 

7+ 

6s  —  5d 

7  + 

4d*5s5p  — ►  4d 

Unidentified 

emissions 

Many 

Many 

Many 

References 

53,85-88 

53,89-95 

96-103 

* 
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an  electron  temperature  of  -  20  eV.  However,  in  this  picture  the  mechanism  of 

formation  of  the  Ar9*  ions  is  not  explained,  since  independent  ion  production 

studies  in  the  collision-free  regime  did  ngt  show  production  of  Ar9*  under  the 

same  conditions  of  irradiation.  It  is  conjectured  that  the  third  harmonic  radiation 
76  77 

at  83  nm  '  produced  in  the  plasma  may  be  playing  a  role  in  the  appearance  of 
the  Ar9*  transitions.  Therefore,  although  much  has  been  accomplished  to  understand 
the  properties  of  these  rare  gas  plasmas,  considerable  work  remains  to  be  done  in 
unraveling  the  details  of  the  observed  spectra,  particularly  in  connection  with 
radiation  originating  from  the  higher  charge  states. 

2.  Fluorescence  from  Ionic  Molecules 

Intense  fluorescence  from  an  intact  molecular  ion  (Na2*)  has  been  observed51 
after  irradiation  of  Na  at  an  intensity  of  -  101*  W/cma.  The  emission  was  strong 
and  somewhat  unexpected,  since  the  probability  for  survival  of  intact  molecular 
species  might  be  presumed  to  be  small  under  such  conditions.  Certain  molecular 
systems,  nevertheless,  are  detected  in  ion  charge  state  spectra,  as  shown  in  Fig. 
(8).  These  finding  are  providing  hints  on  the  details  of  the  complex  nature  of  the 
molecular  coupling  to  intense  radiation. 

The  experimental  results  on  fluorescence  from  Na2*  indicate  the  excitation  of 
inner  electrons  by  a  direct  multiquantum  process.  The  specific  system  illustrating 
this  process  is  the  production  of  2Cg  holes  in  Na,  the  formation  of  which  have 
been  detected  in  two  different  kinds  of  experimental  studies.  The  initial  hint 
came  from  the  measurement47  of  kinetic  energy  distributions  of  ionic  fragments 
produced  in  Na  with  subpicosecond  248-nm  radiation  which  indicated  the  presence 
of  the  charge  asymmetric  channel 

Na2*  — ►  N3*  +  N  (8) 
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resulting  from  multiphoton  ionization.  Interestingly,  in  experiments  performed 

with  subpicosecond  600-nm  radiation,  although  some  differences  in  the  branching 

of  the  fragment  charge  distributions  appear,46,47  additional  evidence  was  also 
109 

indicated  for  reaction  (8). 

Since  other  experiments116  examining  soft  x-ray  produced  fragmentation  of  Na 
showed  that  process  (8)  corresponds  to  the  production  of  a  two-hole  state  having 
an  excitation  energy  of  -  70  eV,  the  obvious  implication  was  that  the  laser  excitation 
produced  a  similar  configuration  of  molecular  excitation.  The  states111,112  implicated 
in  the  x-ray  studies110  of  Na.  which  has  the  ground  state  configuration  (icrg)a(iau)a- 
(2ag)s(2cJu):,(nru)'*(3cg)at  were  lEg[2ag*\  3cg*1].  lru(2ag'*,  2ffu'1],  and  liy2ffg’\ 
liru‘l],  all  involving  a  2ag  hole.  We  note  that  the  2ag  orbital  has  a  binding 

energy  of  37.7  eV,  a  value  considerably  higher  than  the  2ou  (l8.7eV),  nru  (15.6 
eV).  and  3ag  (l6.7eV)  orbitals.  Therefore,  the  Naa*  states  having  a  2ag  vacancy 

and  correlating  to  the  Na*  +  N  asymptote  have  an  energy  substantially  above  the 
lowest  manifold  of  states  in  Naa*  which  are  associated  with  the  N*  +  N*  limit. 

We  describe  below  the  observation  of  a  very  intense  band  of  radiation  at  - 
55.8  nm  which  appears  to  be  fundamentally  related  to  the  mode  of  excitation 
producing  the  charge  asymmetric  fragmentation  represented  by  process  (8).  These 
experiments  were  performed  at  an  intensity  of  -  10l*  W/cma  with  a  pulse  width 
of  ~  600  fs,  conditions  close  to  those  pertaining  to  the  ion  studies  reported  in 
Section  ll.B.2. 

The  55.8-nm  emission  was,  at  the  outset,  highly  unusual  in  one  obvious 

112 

respect.  Although  the  spectrum  of  Na  has  been  studied  for  nearly  a  century 

with  many  techniques  of  excitation,  it  was  not  possible  to  find  any  previous 

53  112  114 

report  ’  ’  of  such  a  characteristic  emission  originating  from  material  composed 

solely  of  nitrogen  atoms  or  molecules  in  either  neutral  or  ionic  form.  From  this 
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fact  alone,  it  appears  that  the  observation  of  the  55.S-nm  band  originates  from  the 

excitation  of  a  new  class  of  states  and  that  the  strong  production  of  these  levels 

is  intimately  connected  with  the  basic  nature  of  the  multiquantum  mechanism 

involved  in  the  coupling  to  the  molecule.  This  line  of  reasoning  and  its  implications 

in  regard  to  the  generation  of  x-rays  and  the  comparison  with  collisional  processes 

are  explored  below  in  Sections  III. A  and  III.C,  respectively. 

The  hypothesis  that  emerged  from  the  analysis  of  the  55.8-nm  radiation  is 

that  it  originates  from  a  quasi-bound  3TTU  level  of  N22*  having  a  (20g'l,1iru‘l) 

configuration  which  undergoes  a  bound-bound  transition  terminating  on  the  D  3TTg 

51 

level  associated  with  the  N*  +  N*  asymptote  In  addition,  the  previous  studies 

47 

of  ionic  fragmentation,  combined  with  the  comparison  of  the  results  under 

comparable  experimental  conditions  with  N20,  indicate  that  the  20g  vacancy  is 

produced  by  a  direct  multiphoton  process  in  Na  leading  to  an  excited  state  of  the 

type  (N2*  +  N),  namely,  a  level  with  evident  charge-transfer  character.  Therefore, 

the  findings  concerning  the  fluorescence  are  fully  consistent  with  the  results 

47 

independently  derived  from  the  ion  studies 

Overall,  several  transitions  were  assigned  arising  from  the  production  of  the 
20g  vacancy.  Qualitatively,  the  relative  strengths  of  the  molecular  features  could 

be  accounted  for  by  the  overlap  of  the  wavefunctions  involved  in  the  transitions. 
Basically,  the  transition  producing  the  55.8-nm  radiation  is 

2og'1  —  2ou'1  +  Y,  (9) 

a  case  which  has  a  very  favorable  overlap,  a  fact  that  can  be  seen  directly  by 
visualization  of  the  wavefunctions  of  those  two  orbitals.1 15,1 16 

The  spectrum  in  the  region  ~  56  nm  is  shown  in  Fig.  (12).  Particularly 

significant  from  the  standpoint  of  the  identification  is  the  measured  width  of  the 
molecular  feature,  with  respect  to  the  widths  of  known  N*  lines,  and  the  accompanying 


34 


presence  of  a  vibrational  progression  (v  =  0,  1,  2).  The  molecular  energy  level 

diagram  derived  from  the  analysis51  is  illustrated  in  Fig.  (13).  Finally,  Table  II 

summarizes  the  assignments  of  the  molecular  features  observed. 

This  spectroscopic  evidence  carries  with  it  an  important  implication.  Combined 

47 

with  the  results  stemming  from  the  ion  studies,  the  indication  is  that  the  2ag 

vacancy  is  produced  by  a  direct  multiphoton  process  in  Na  leading  to  an  excited 

state  of  (N3*  +  N)  character.  Indeed,  in  referring  to  Table  II,  the  presence  of 

one  transition  (59.15  nm)  implies  that  the  2ag  electron  can  be  removed  without  the 

additional  loss  of  any  valence  electrons.  For  this  to  occur,  the  removal  of  the 

outer  more  weakly  bound  electrons  must  be  sufficiently  slow  so  that  the  interaction 

can  communicate  the  energy  needed  to  excite  the  2cg  orbital  before  the  outer  shells 

26  43  11 7—1 1 9 

are  ionized.  Controversy  has  existed  on  this  point  for  some  years.  ’  ’ 

There  is,  however,  some  more  recent  theoretical  analysis  indicating  that  such  a 

suppression  of  the  ionization  rate  may  be  able  to  occur  at  sufficiently  high 

.  120-123 

frequency. 

3.  X-Ray  Generation  from  Solid  Targets 

Solid  materials  offer  the  possibility  of  combining  a  high  particle  density  with 

a  very  high  level  of  electronic  excitation,  if  the  deposition  of  energy  in  the  solid 

can  be  made  to  occur  sufficiently  rapidly.  High  brightness  subpicosecond  light  sources 
of  the  type  shown  in  Fig.  (1)  now  enable  such  rapid  excitation  to  be  achieved. 

Given  the  fact  that  the  studies  of  ionization  described  in  Section  II. B  indicated 
total  energy  investments  of  ~  1  keV  per  atom,  at  least  for  the  heavier  materials, 

the  immediate  suggestion  is  that  a  properly  irradiated  solid  could  serve  as  a 

powerful  source  of  x-rays  in  the  kilovolt  region. 

Consider  a  slab  of  unit  area  with  density  p  and  thickness  6  which,  upon 

excitation  with  intense  subpicosecond  radiation,  produces  radiation  in  the  kilovolt 
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Fig.  (12):  The  spectrum  observed  in  Na  in  the  56-nm  region  produced  with 

intense  subpicosecond  {-  600  fs)  248-nm  irradiation.  Two  groups 
of  molecular  features  are  identified,  including  the  main  transition 
at  55  6-nm,  as  well  as  several  atomic  lines.  The  density  of  N, 
used  was  -  101*  cm'1 
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Fig.  (13):  Potential  energy  curves  for  N,  for  selected  states  of  Interest. 

Two  observed  transitions  involving  a  20  hole  are  Indicated.  Data 
for  the  curves  were  taken  from  a  variety  of  sources  (Refs.  109-1 11) 
as  well  as  the  author's  own  estimates.  Dissociation  limits  are 
given  for  reference. 
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TABLE  II 


Assignment  of  Observed  Transitions  to 
States  of  Molecular  Nitrogen  Ions 


Observed 

Wavelength 

(nm) 

Energy 

(eV) 

Identification* 

Vibrational* 
Levels  (v-V) 

Electronic 

Configuration** 

67.57 

18.35 

N22*:A'  3nu  -  G  3TTg 

0  -  1 

2og  1  lVl"2au'1  1lru'1 

59.15 

20.96 

Njl*:X'  2Ig-  A  2nu 

0-0 

V'iv‘ 

57.16 

21.69 

N22*:A'  3TTu  -  E  3ng 

0  -  1 

2tJg‘l  IV1"20^1  1lfu‘l 

57.02 

21.74 

n22*:A'  3nu-  E3ng 

0-0 

2og  1  1V-V  1lrU*1 

55.84 

22.20 

N22*:A'  3TTu  -  D  3flg 

0-0 

20g  1  1iru‘l -2au'1  1iru‘l 

55.21 

22.46 

n22*:A'  3nu-  d  3ng 

1  -  0 

20g-l1V1-2Cu-llVl 

54.76 

22.64 

N22*:A'  3nu  -  D  3T?g 

2-0 

2°g’1  l^u'1  -  20U  1  11TU'1 

48.56 

25.53 

N22*:A'  3nu-  A3Zg* 

0-0 

V’V-  *** 

*  Using  calculations  from  Ref.  (Ill)  and  authors'  estimates 

**  Determined  from  Ref.  (109) 

***  Lower  state  not  included  in  Ref.  (109) 
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spectral  range.  In  order  to  set  the  scale,  we  want  to  obtain  a  simple  upper  bound 

on  the  achievable  x-ray  intensity.  Since  the  time  of  excitation  is  short,  as  a 

first  approximation,  we  will  ignore  the  hydrodynamical  expansion,  a  process  which 
will  occur  roughly  at  a  rate  of  -  1-10  A/fs.  An  upper  limit  on  the  x-ray  intensity 
Ixm  can  then  be  simply  written  as 

h<i)xp6 

Ixm - .  00) 

2Tx 

in  which  ftaix  represents  a  kilovolt  quantum  and  tx  is  the  effective  radiative  rate 

of  the  excited  material.  For  the  present  discussion,  we  will  set  p  as  solid  density 
(~  1.7  x  10”  cm'3)  and  take  as  the  lower  limit  for  xx  a  typical  radiative  rate  for 

kilovolt  transitions.  For  the  latter  will  use  the  range  2  x  10'1*  $  tx  ^  10'13  s 

as  a  rough  measure.124  If  we  take  for  6  the  recently  measured  value85  of  the  depth 
of  energy  penetration  of  -  250-300  nm  for  MgF2/Si02  targets  irradiated  with  248 

nm  at  an  intensity  of  -  3  x  10l*  W/cm2,  as  described  in  Appendix  A,  Eq.  (2)  can 

be  evaluated.  With  these  parameters  we  find  3  x  101"  ^  Ixm  ^  2  x  10IS  W/cm2. 

125 

This  simple  upper  limit  on  Ixm  can  now  be  compared  to  the  experimental  result, 
as  described  in  Appendix  B. 

An  experimental  concern  in  studies  of  the  interaction  of  intense  radiation 
with  solid  matter  is  the  influence  of  any  low  intensity  prepulse  on  the  true  state 
of  the  irradiated  system.  As  discussed  in  Appendix  C,  it  has  been  experimentally 
4ound126  that  the  interference  of  the  prepulse  can  be  entirely  eliminated  at  248  nm 
by  the  use  of  a  target  composed  of  ultraviolet  transmitting  materials.  For  example, 
MgF2,  CaF2,  and  BaF2  exhibit  essentially  no  absorption  up  to  relatively  high  intensities 
and  studies  have  demonstrated  that  these  values  are  sufficiently  high  to  prevent 
any  significant  influence  by  the  prepulse.  This  enables  the  intense  subpicosecond 
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energy  to  couple  directly  to  virgin  solid  matter.  ’ 

The  experimental  configuration  used  for  the  measurement  of  the  x-ray  emission 
from  solid  BaFa  targets  is  described  fully  in  Appendix  B.  The  spectrum  observed 
from  BaFa  irradiated  at  a  peak  intensity  of  ~  0.3-1 .0  x  to17  W/cm2  exhibited 
strong  emission  in  the  9-13  A  range  as  shown  in  Fig.  (14).  The  positions  of 
three  M-shell  thresholds  for  Xe  (Z  =  54)  are  also  indicated  for  reference.  In 
addition  to  the  spectrum,  calibration  of  the  film  revealed  that  -  2-4  mJ  of  x-rays 
were  produced  on  each  shot.  These  measurements  showed,  assuming  isotropic 
emission  from  the  plasma,  that  the  experimentally  radiated  intensity  Ir  was 
approximately  Ir  ~  (0.5  x  1.0)  x  10IS  W/cm2.  We  see  immediately  that  this  value 
falls  in  the  range  given  by  the  estimate  of  the  upper  bound  Ixm.  The  measured 
performance,  therefore,  appears  very  close  to  the  maximum  possible.  Although 
such  a  source  is  incoherent,  a  radiation  rate  of  this  magnitude,  coupled  with  the 
relatively  narrow  spectrum  shown  in  Fig.  (14),  constitutes  a  source  of  extremely 
high  spectral  brightness  that  can  produce  unusually  high  excitation  rates  in  absorbing 
material.  We  consider  the  case  of  xenon  as  an  example,  a  material  whose  M-shell 
absorptions,  as  shown  in  Fig.  (14),  are  very  well  matched  to  the  measured  BaFa 
spectrum. 

Consider  the  response  of  xenon  atoms  located  close  to  the  x-ray  source. 
The  production  rate  of  M-shell  vacancies  is  given  approximately  by 
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in  which  (-  2  x  10  cm2)  denotes  the  absorption  cross  section  of  xenon 
at  -  1  keV.  From  Eq.  (11)  we  obtain  ~  10ihs’\  a  rate  which  is  not  far 
below  the  Auger  decay  rates  of  such  M-shell  vacancies.  For  example,  an  Ma 
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Fig.  (14):  Spectrum  obtained  from  BaF,  at  an  intensity  of  -  0.3-1. 0  x  1017 

W/cma.  The  positions  of  the  Mp  Ma,  and  M,  thresholds  of  Xe 
are  indicated  for  reference.  For  details  see  the  text  and  Appendix 
B. 
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vacancy  has  a  theoretically  estimated  width  ’  in  xenon  of  r(Ma)  =  4.83  eV,  a 
value  which  corresponds  to  a  decay  time  t  -  0.4  fs.  The  probability  Pa  of  producing 
a  direct  double  vacancy  is  then  on  the  order  of 

P*  ~  ($Mt)S  ~  10'*  •  (12) 
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Since  such  double  M-shell  vacancy  states  would  be  expected  to  have  Auger 

132 

patterns  significantly  different  from  both  single  vacancies  or  double  vacancies 

produced  sequentially,  a  fraction  of  this  magnitude  (-  10'*)  may  be  sufficiently 

large  for  experimental  detection.  With  very  few  exceptions,  such  shell-specific 

1 33— 

double-vacancy  states  constitute  a  virtually  unexplored  class  of  excited  matter. 
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III.  DISCUSSION  OF  NEW  POSSIBILITIES 

The  production  of  bright  x-ray  sources,  as  stated  in  the  Introduction,  is 
centered  on  the  controlled  and  rapid  production  of  x-ray  emitters  in  an  appropriate 
spatial  volume.  In  this  context  three  concepts  are  currently  beginning  evaluation. 
One  involves  the  behavior  of  highly  excited  molecules,  and  is  related  to  the 
discussions  in  Sections  II.B.2,  II.C.2,  the  second  involves  the  studies  of  solids,  as 
described  in  Section  II. C. 3,  and  the  final  top’c  concerns  modes  of  electromagnetic 
propagation  occurring  in  the  strong-field  regime.  While  not  fully  examined  and 
speculative,  on  the  basis  of  the  experimental  and  theoretical  results  currently 
available,  these  topics  clearly  represent  fertile  and  important  new  areas  of  active 
research. 

A.  Molecular  Processes  for  X-Ray  Generation 

The  results  described  in  Sections  II.B.2  and  II.C.2  above  lead  directly  to  the 
formulation  of  a  new  principle  of  molecular  design  for  x-ray  emitting  systems. 
The  specific  findings  that  lead  to  this  conclusion  are  summarized  in  Table  III. 


The  properties  of  the  strong-field  coupling  compiled  in  Table  III  are  fundamental 
to  the  production  of  coherent  radiation  of  x-ray  wavelengths  primarily  because 


they  enable  the  combination  of  (1)  very  highly  excited  forms  of  matter  with  (2) 
the  characteristics  of  dense  cold  materials.  This  aspect  of  the  coupling52  is 
considered  further,  particularly  in  comparison  with  collisional  process.52  in  Section 
III.B  of  Appendix  D.  Such  circumstances  are  exceptionally  conducive  to  the 
generation  of  amplification  at  x-ray  wavelengths.  It  is  the  achievement  of  this 
conventionally  paradoxical  condition,  which  can  be  made  to  occur  in  the  molecular 
frame,  that  permits  the  controlled  generation  of  the  very  high  rate  of  excited 
state  production  necessary  for  x-ray  amplification.  Moreover,  since  molecular  systems 
are  utilized,  specific  molecular  designs  can  be  created  whose  properties  selectively 
govern  the  flow  of  energy  into  particular  excited  states.  The  resulting  wavelength 
of  emission  will  thereby  be  dependent  upon  the  detailed  structure  of  the  molecular 
system. 

These  considerations  are  now  illustrated  within  the  framework  of  a  particular 

example.  Since  we  desire  a  high  energy  deposition  localized  on  a  particular  atom. 

the  results  of  the  work  described  above  indicate  that  this  may  be  accomplished 

by  placing  a  relatively  heavy  atom  on  a  line  between  two  lighter  systems.  We 

now  assume  that  the  axis  defined  in  this  way  coincides  with  the  polarization  vector 

of  the  incident  radiation.  In  order  to  take  advantage  of  the  anisotropy  of  the 
44  46  47  50 

interaction  '  '  ,  we  now  place  two  atoms,  one  on  either  side  adjacent  to  the 
heavy  atom  on  a  line  perpendicular  to  the  original  axis.  In  this  form,  the  system 
has  the  shape  of  a  cross.  Since  it  is  inconvenient  to  orient  molecules,  we  now  add 
two  more  atoms  extending  on  either  side  of  the  central  atom  out  of  the  plane  of 
the  cross.  If,  for  simplicity,  we  let  the  six  lighter  systems  be  equivalent,  the 
result  is  an  octahedral  molecule  (Of,  symmetry)  with  the  heavy  atom  in  the  center. 
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TABLE  III 


Features  of  Strong-Field  Multiquantum 
Interaction  Which  Lead  to  Molecular  X-Ray  Emitter  Design 


Property  of  Interaction  References  to  Discussion 


(1).  Fields  e  >>  e/a \ 

(1). 

Ref. 

(9)  and  Ref.  (10) 

can  be  readily  produced 
(2).  Rapid  energy  transfer  to 

(2). 

Ref. 

(47)  and  Ref.  (51) 

molecule,  At  ~  X/c, 

Ar  small  -  0.1  ^ 

(3).  Coupling  exhibits  strong 

(3). 

Ref. 

(44),  Ref.  (46),  and 

anisotropy,  mainly  parallel 

Ref. 

(50),  and  Ref.  (107) 

to  molecular  axis 

(4).  Substantial  e'  displacements  Ax 

(4). 

Ref. 

(47)  and  Ref.  (51) 

driven  for  e  >  >  e/a£  Ax  >  r0 
(5).  Inner  molecular  shells  can  be 

(5). 

Ref. 

(51) 

directly  excited  (e.g.  Na(2Cg)]. 

(6).  Site  selective  molecular  energy 

(6). 

Ref. 

(42)  and  Section 

deposition  (N20/N2). 

2.b 

above. 

Specific  examples  would  be  SF#,  MoFt,  and  UF„  all  of  which  are  convenient  gas 
phase  species  at  room  temperature.  A  molecule  with  symmetry,  of  course,  will 
always  present  an  axis  involving  the  two  lighter  atoms  with  the  central  heavy 
atom  which  is  favorably  aligned  with  the  polarization  vector  of  the  incident  radiation. 

We  now  consider,  as  a  specific  example,  the  behavior  of  a  molecule  (SFs)  as 
a  result  of  intense  irradiation  in  light  of  the  considerations  described  above.  In 
order  to  perform  this  analysis,  we  introduce  a  simplification  that  enables  the 
system  to  be  regarded  as  a  diatomic  molecule.  Hence,  we  will  assume  that  the 
central  S  atom  interacts  primarily  with  a  single  F  atom.  A  slight  breaking  of  the 
symmetry  occurring,  for  example,  through  vibrational  excitation  could  be 
envisioned  as  the  physical  basis  for  this  dynamical  choice.  Although  we  anticipate 
that  the  remaining  F  atoms  will  influence  the  molecular  behavior,  we  expect  that 
many  of  the  important  features  of  the  dynamics  will  be  revealed  by  the  simpler 
diatomic  picture. 
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Consider  the  SF  diatom  in  the  charge  state  SF17*  as  a  simple  illustration. 
Oue  to  the  ability  to  produce  site-selective  energy  deposition  (point  (6)  in  Table 
III),  we  will  assume  that  the  ionization  is  distributed  as  Sia*  -  F*\  For  ground 
state  atoms,  this  gives  the  configurations  Sia*(lsa2sa)  and  F**(lsa2sa),  both  of 
which  give  ‘So  terms  and,  consequently,  a  single  lZ*  molecular  state.  Furthermore, 
since  the  energy  transfer  producing  the  ionization  occurs  very  rapidly  (point  (2) 
in  Table  III,  also  see  Section  III.B  in  Appendix  D),  we  take  the  initial  internuclear 
separation  of  the  SF17*  system  as  equivalent  to  the  S-F  internuclear  distance 
occurring  in  the  neutral  parent  SF,  molecule.  This  distance  is  r0  =  1.56  A.  The 
point  (A)  representing  this  situation  can  now  be  placed  on  the  molecular  energy 
level  diagram  shown  in  Fig.  (15)  along  with  the  corresponding  lZ*  coulomb  state. 
From  this  initial  configuration,  point  (A),  the  system  will  begin  to  evolve  along 
the  ‘Z*  curve.  Furthermore,  since  essentially  no  momentum52  was  communicated 
to  atoms  in  the  molecule  during  the  strong-field  multiquantum  excitation,  the 
dissociative  motion  along  the  lZ*  curve  begins  with  a  vanishing  velocity  (v  =  0). 
As  the  diatomic  system  dissociates,  curve  crossings  will  be  encountered,  such  as 
the  one  labeled  (B),  shown  arising  from  the  (Su*)*  +  F**  curve. 

Several  important  aspects  of  the  process  are  illustrated  in  the  diagram  of 
Fig.  (15).  They  are  (1)  a  crossing  to  the  ground  state  Su*+  F**  configuration 
is  ngl  accessible,  (2)  the  crossing  leading  to  excited  (Sa*)*  by  electron  transfer 
is  reached,  since  r;  >  r#,  (3)  the  velocity  at  point  (B)  can  be  very  low,  and  (4) 
the  dynamics  of  the  crossing  will  be  governed  only  by  radial  couplings. 

1 38 

The  behavior  at  such  crossings  is  well  described  in  the  Landau-Zener  picture 
The  probability  P  of  going  though  the  crossing  and  remaining  on  the  original 
curve  is  given  by 

P  aexpf-(il-)( - l¥ui - )j,  (13) 

[  \2irvA  d/dR^-E,)  /J 
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in  which  R  is  the  internuciear  coordinate,  v  is  the  velocity,  Via  represents  the 
matrix  element  for  the  radial  coupling,  and  Et(R)  and  Ea(R)  denote  the  energies 
of  the  participating  states.  We  note  immediately  that 

Mm  P  =  0,  (14) 

v-*0 

so  that  efficient  crossing  is  expected  to  the  (S“*)*  +  F**  curve,  which  generates 
highly  excited  (S11*)*  ions,  since  v  is  small.  Furthermore,  such  crossings  confer 
substantial  selectivity  on  the  state  (S11*)*  produced  through  the  properties  of  the 
radial  matrix  element  Vla.  For  example,  since  the  interaction  is  radial,  no  angular 
momentum  is  coupled  to  the  electrons  and  crossings  are  only  permitted  to  other 

.  54 

curves  of  the  same  (I  )  symmetry.  This  is  a  powerful  selection  rule  for  the 

selective  channeling  of  energy  into  a  restricted  set  of  excited  quantum  states. 

Experimental  evidence  has  shown  that  such  mechanisms  can  be  quite  effective  in 

139 

selectively  generating  excited  states  from  both  molecular  and  quasi-molecular/coi- 
140 

lisional  systems.  Finally,  for  the  states  shown  in  Fig.  (15)  representing  the 
(S11*)*  levels,  the  corresponding  wavelength  is  Xx  =  40  A  (-  300  eV). 

The  most  important  aspect  of  Fig.  (15)  is  the  fact  that  the  nature  of  the 
particle  interactions  is  such  that  no  single  step  enables  a  direct  transition  from 
the  initially  excited  configuration  to  the  around  state  of  the  atomic  system  of 
interest  (Su*).  Therefore,  the  dynamics  of  the  excited  material  force  the  energy 
to  flow  in  an  orderly  pathway  that  leads  inevitably  to  the  excited  states  of  the  system. 
For  the  free  ion,  the  only  available  downward  channel  involves  radiation.  In 
summary,  a  dominant  energy  pathway  in  the  mechanism  of  molecular  relaxation 
of  the  excited  system  is  established  which  contains  a  segment  that  can  be  traversed 
only  by  the  emission  of  radiation.  When  this  general  circumstance  is  present, 
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Fig.  (15):  Partial  potential  energy  diagram  of  Sq/  +  Fqa*  (qx  +  q,  -  17) 

system.  Note  that  the  electron  transfer  crossing  at  rj  is  such 
that  r;  >  r0.  For  Su*,  Xx  =  40  A  (~300eV).  The  ordinate  is  not 
drawn  to  a  linear  scale. 
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stimulated  emission  from  the  resulting  highly  nonequilibrium  distribution  of  states 

can  be  readily  generated.  This  situation  is  directly  analogous  to  that  occurring  in 
141 

excimer  media 

B.  Hard  X-Ray  Generation  from  Solid  Targets 

As  discussed  in  Section  II. C. 3  and  Appendix  B,  very  intense  emission  in  the 
9-13  A  range  arising  from  M-shell  transitions  in  Ba  ions  was  observed  from  BaF2. 
Indeed,  the  strength  of  the  observed  emission  was  so  great  that  it  was  only 
slightly  below  an  estimate  of  the  upper  bound  possible.  The  model  used  to  develop 
the  estimate  of  the  upper  bound  was  both  simple  and  extreme.  Basically,  it 
required  that  every  atom  in  the  radiating  volume  radiate  a  kilovolt  quantum  in 

the  radiative  lifetime  of  the  Ba  ion.  This  type  of  picture,  obviously,  represents 
an  extraordinary  level  of  excitation.  In  such  circumstances,  it  is  natural  to  inquire 
if  excitation  of  deeper  lying  states  might  also  be  occurring,  albeit,  at  a  significantly 
lower  rate. 

For  Ba,  the  next  level  of  excitation  would  correspond  to  the  L— shell,  the 

edges128  of  which  lie  in  the  5-6  keV  region  (~  2  A).  In  order  to  check  for  the 
presence  of  such  radiation  in  this  range,  we  equipped  the  spectrometer  with  a  Si 
(111)  crystal  and  examined  the  region  -  2  A  experimentally.  The  spectrum  produced 
in  this  preliminary  study  is  shown  in  Fig.  (16).  The  film  shows  a  definite  exposure 
in  the  region  (shaded)  -  2  A.  We  note  that  the  calibration  of  the  wavelength 
may  be  uncertain  as  much  as  -  0.5  A.  The  recorded  strength  of  the  exposure  is  - 

1-2%  of  the  kilovolt  (9-13  A)  spectrum  shown  in  Fig.  (14).  Although  the  detected 

emission  lies  in  the  anticipated  spectral  region,  the  lineshape  is  not  of  the  expected 
form  and  considerable  further  work  needs  to  be  done  to  clarify  the  origins  and 
properties  of  this  observed  radiation.  Obvious  points  of  interest  are  the  scaling 
relationships  on  intensity  and  atomic  number  (Z)  governing  the  production  of 
emission  in  the  multi-kilovolt  range. 
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C.  General  Consideration  of  Multiquantum  Strong-Field  Excitation 

Two  important  conclusions  concerning  the  character  of  strong-field  multiquantum 
coupling  stem  from  the  discussion  given  in  sections  above.  They  involve  (t)  the 
ability  to  selectively  excite  inner  electrons  and  (2)  the  general  nature  of  the 
strong-field  coupling  to  molecules. 

1.  Inner-Electron  Excitation 

As  described  in  Section  II. C. 2,  the  observation  of  new  band  radiation  at  55.8 
nm  from  transitions  involving  inner-orbitals  in  NT  and  N2*,  produced  by  a  direct 
multiphoton  process  involving  N2  in  reactions  of  the  form 

nY  +  N2  (X  lZg)  -  N;  (20g*1)  +  e*  (15) 

and 

nY  +  N2(X  lEg)  -  N2a*  (2<Jg'\  1ttu'1)  +  e*  +  e‘,  (16) 

provides  strong  evidence  that  deeply  bound  electrons,  in  this  case  2ag,  can  be 
directly  excited117,119  by  strong-field  processes  without  extensive  loss  of  more  weakly 
bound  outer  valence  electrons.  It  is  expected  that  this  feature  of  the  interaction 
will  manifest  itself  in  many  materials.  Conventional  pictures  of  multiphoton 

processes,  particularly  those  involving  the  sequential  emission  of  outer  electrons 
from  the  system,  do  not  provide  a  basis  for  understanding  these  results.  Moreover, 
although  N2  has  been  one  of  the  most  common  subjects  of  spectroscopic  analysis 
for  nearly  a  century,  it  is  striking  that  the  55.8  nm  band  was  one  of  the  strongest 
features  in  the  spectrum,  but  had  never  previously  been  observed.  These  facts 
suggest  the  presence  of  an  important  new  feature  in  the  mechanism  of  coupling. 
We  now  examine  that  aspect  with  the  use  of  a  simple  approximate  model  of  the 
interaction. 

2.  General  Features  of  Strong-Field  Coupling 

In  addition  to  providing  specific  data  on  N2  and  its  ions,  the  information 
obtained  from  the  studies  of  ion  production47  and  fluorescence  described  above  in 
Section  II. C. 2  leads  to  an  important  general  hypothesis  concerning  the  basic  character 
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P'9-  (16):  Spectrum  of  BaF,  in  the  region  ~  2  A.  The  wavelength  scale  may 

be  uncertain  in  the  amount  of  -  0.5  A.  The  visible  exposure  on 
the  film  corresponds  to  the  shaded  zone.  The  peak  at  -  4.8  A 
was  caused  by  damage  from  the  film  holder. 
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of  molecular  multiphoton  interactions  in  the  strong-field  (E  >  e/a02)  regime. 
Through  comparison  with  electron  and  ion  collisional  processes,  it  will  be  shown 
that  there  is  a  large  class  of  excitations  which  can  be  strongly  excited  only  by 

strong-field  coupling.  Therefore,  new  states  of  excitation  uniquely  produced  by 
the  strong-field  mode  of  interaction  are  generally  expected  for  all  molecular 
materials  regardless  of  structure  or  composition.  Furthermore,  it  will  be  seen 
that  the  excited  states  produced  naturally  by  this  form  of  interaction  would  commonly 
embody  highly  asymmetric  charge  distributions. 

The  comparison42  with  charged-particle  collisions  can  be  accomplished  in  two 
ways,  both  of  which  express  different  aspects  of  the  dynamics  of  coupling.  One 
examines  the  pattern  of  the  electric  force  developed  during  the  collision  while 

the  other  represents  the  encounter  graphically  in  terms  of  variables  describing 
(1)  the  maximum  electric  field  em  produced  and  (2)  the  energy  transfer  Et 
communicated  to  the  molecular  system  during  the  collisional  event. 

We  now  describe  a  very  simple  approximate  picture  in  order  to  present  this 

comparison.  Consider  the  ion-molecule  (N2)  collisional  process  illustrated  in  Fig. 

(17a)  involving  an  ion  having  charge  Ze  and  velocity  v  colliding  with  impact  parameter 
b.  We  assume  that  the  ion  follows  a  classical  straight-line  trajectory  approximately 
perpendicular  to  the  molecular  axis.  As  the  collision  proceeds,  the  electric  field 
from  the  ion  experienced  by  the  molecule  varies  appreciably  in  both  magnitude  and 
direction.  Only  at  point  B  on  the  trajectory  is  the  field  aligned  aiong  the  molecular 
axis.  Otherwise  the  direction  of  the  force  on  the  molecular  electrons  deviates  con¬ 
siderably  from  the  molecular  axis.  The  situation  can  be  made  more  favorable  for 
the  excitation  of  electronic  motions  along  the  molecular  axis  by  going  to  large 
impact  parameter  (b  >  >  r0),  but  then  the  maximum  collision  field 

e  =  —  (17) 

Em  b2 

is  greatly  reduced. 
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The  electric  field  of  a  short-duration  (Tp)  intense  pulse  of  radiation  polarized 
along  the  molecular  axis  represents  a  qualitatively  different  type  of  interaction. 
In  this  case,  the  resulting  force  field  acts  primarily  to  excite  electronic  motions  along 
the  molecular  axis.  For  a  sufficiently  short  pulse,  nominally  in  the  ip  <  1  ps  range, 
the  effect  of  molecular  rotation  for  thermally  excited  levels  at  normal  temperatures 
is  negligible  for  almost  ail  molecular  materials.  This  can  be  seen  by  noting  that 
molecular  rotational  frequencies  vr  are  given  approximately  by  =  2BJ  for  rotational 
constant  B  and  angular  momentum  J.  With  the  exception  of  the  hydrides,  if  we 
take  B  ~  1  cm'1  as  a  typical  value,  vrxp  <<  1  even  for  J  =  10. 

Further  comparison  is  shown  in  Fig.  (17b)  for  the  case  involving  a  parallel 
orientation  of  the  ionic  trajectory  and  the  molecular  axis.  The  main  effect  in 
this  orientation  is  a  force  transverse  to  the  molecular  axis.  Although  the  field 
can  have  an  appreciable  component  along  the  molecular  axis  in  the  limit  b  -  0, 
this  situation  has  a  vanishing  probability  of  collision.  In  addition,  as  shown 
below,  such  close  encounters  with  charged  particles,  particularly  ions,  differ  from 
the  strong-field  coupling  in  another  important  way,  since  they  generally  impart  a 
substantial  recoil  energy  to  the  molecular  system. 14^  In  sharp  contrast,  the  strong- 
field  coupling,  although  capable  of  exerting  a  force  with  a  magnitude  representative 
of  ion  collisions,142  is  extremely  ineffective  in  transferring  kinetic  energy  directly 
to  the  target  molecule  or  its  constituent  atoms. 

The  differences  between  the  charged-particle  collisional  processes  and  the 
strong-field  radiative  interaction  can  be  further  understood  by  determining  the 
relationship  between  kinetic  energy  transferred  to  the  target  system  and  the 
strength  of  the  peak  electric  field  occurring  in  the  interaction.  To  find  this 
relationship,  these  two  processes  can  be  parameterized  by  the  peak  field  strength, 
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Fig.  (i  7a):  Geometry  of  ion-molecule  (N,)  collisions  showing  the  changing  electric 

fields  experienced  by  the  target  during  the  colllsional  event,  (a)  Ionic 
trajectory  perpendicular  to  molecular  axis;  (b)  ionic  trajectory  parallel 
to  molecular  axis. 
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em,  and  the  effective  duration  of  the  interaction,  t.  For  a  given  field  strength, 
we  set  the  interaction  times  equal  for  the  two  processes  so  that  the  comparison 
is  made  for  a  fixed  transition  probability.  In  the  case  of  the  strong-field  process, 
the  studies  of  ion  fragmentation47  indicated  that  the  interaction  time  can  be 
taken  as  approximately  comparable  to  a  period  of  the  electromagnetic  wave  t  = 
X/c,  while  the  corresponding  field  strength  is  given  by  em  =  (4irl/c)1/2  in  which  I 
is  the  laser  intensity. 

For  the  charged  particle  collisions,  we  consider  a  process  of  the  form 

AZl+  +  B  — ►  AZl+  +  BZa  +  +  Zse‘  (18) 

in  which  a  projectile  ion  A  of  charge  Zx  collides  with  a  neutral  system  B  at  rest 

producing  a  recoiling  ion  BZa  +  with  the  collateral  generation  of  Z2  electrons. 

The  recoil  energy  Et  of  the  target  system  B,  for  a  projectile  with  energy  E  at 

1 43  1 44 

impact  parameter  b,  is  given  ’  approximately  by 

ZX%VMA 

Et  =  -  ,  (19) 

EMgb2 

with  Ma  and  Mg  representing  the  masses  of  the  particles  A  and  B,  respectively, 
and  e  denoting  the  electronic  charge.  The  time-scale  of  the  collisional  interaction 
will  be  given  by 


The  combination  of  Eqs.  (17),  (19),  and  (20)  gives  the  desired  relationship  between 
Et  and  em,  namely, 


Interestingly,  the  energy  transfer  Et  is  independent  of  the  nature  of  the  charged 


projectile.  The  final  connection  between  the  strong-field  interaction  and  the 
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charged-particle  collision  is  obtained  by  setting  the  time  scale  of  the  collision 
equal  to  X/c  which  yields 


where  f  is  the  maximum  field  strength  in  atomic  units. 

The  validity  of  Eq.  (22),  for  the  simple  collisional  picture  used,  requires  that 


<<  1 


(23) 


MAvA 


where  Aq  is  the  momentum  transfer  occurring  in  the  collision  and  vA  is  the 
velocity  of  the  projectile.  Using  Eqs.  (17),  (21),  and  (22),  we  find  that 


If  we  take  Z22 /Zx  ~  1  as  a  characteristic  example,  Eq.  (23)  is  satisfied  for  f  ^  30. 
Actually,  for  field  strengths  larger  than  ~  30  atomic  units,  the  energy  transfer 
will  actually  be  greater  than  that  given  by  Eq.  (22).  Thus,  Eq.  (22)  represents  a 
lower  bound  for  charged-particle  collisions.  Furthermore,  since  the  characteristics 
of  the  projectile  do  not  appear  in  Eq.  (22),  this  relationship  is  also  valid  for 
electron  collisions. 

We  now  examine  the  magnitude  of  other  physical  processes  which  can  impart 
energy  transfer  to  atomic  or  molecular  systems.  Generally,  for  the  electron 
collisions,  the  emission  of  Auger  electrons  will  impart  a  recoil  energy  exceeding 
that  arising  from  the  incident  electron,  the  magnitude  of  which  can  be  readily 
estimated  from  the  kinematics  of  two-body  decay.  If  we  assume  that  an  electron 
with  mass  me  is  emitted  with  an  Auger  energy  Ea,  then 


In  order  to  estimate  an  upper  bound,  we  will  assume  an  energy  scale  for  Ea  on 
the  order  of  the  K-edge.  In  this  case,  for  light  materials  in  which  the  Auger 
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yields  arising  from  a  K-vacancy  are  high,  a  system  such  as  Ne  (hwK  =  870  eV), 
gives  Et  *  24  meV,  a  value  almost  exactly  equal  to  the  mean  thermal  energy  at 
300  K. 

For  heavy  materials,  such  as  uranium,  since  radiative  yields  dominate  the 
decay  of  K-shell  vacancies,  the  recoil  developed  is  given  approximately  by 


Et 


(^K)a 

2Mc3 


(26) 


where  fiaiK  is  the  emitted  photon  energy.  For  U,  Eq.  (26)  gives  Et  =  30  meV,  a 
magnitude  close  to  that  arising  from  Auger  decay  in  light  elements.  L— shell 
Auger  decay  from  U  (-  21  keV)  gives  Et  =  46  meV,  so  it  is  also  in  the  same 
range.  Less  energetic  excitations  arising  from  electron  scattering  will,  of  course, 
produce  smaller  recoil  energies.  Overall,  for  a  wide  range  of  scattering  conditions 
and  materials,  the  recoil  energies  induced  by  electron  scattering  will  fall  in  the 
range  -  3  x  10'*  ^  Et  ^  5  x  10*3  eV. 

The  recoil  energy  imparted  to  an  ion  during  multiphoton  ionization  is  naturally 
expected  to  be  small,  since  the  incident  photons  carry  very  little  momentum  For 
a  given  number  N  of  quanta  absorbed  with  energy  fiu,  the  recoil  is  expressed  as 


(Nhu)3 

Et  =  2Mc3 


(27) 


the  multiquantum  counterpart  of  Eq.  (26).  From  the  previous  discussion,  for  a 
total  energy  such  that  Nhw  =  's  *ess  l^an  30  me^  *or  materials. 

Two  additional  mechanisms  can  communicate  recoil  energy  to  the  ion  produced 
by  multiphoton  ionization.  They  are  (1)  the  interaction  of  the  field-driven  ionized 
electron  with  the  ion,  a  process  equivalent  to  inverse-bremsstrahlung  in  the  field 
of  the  ion,  and  (2)  the  subsequent  emission  of  Auger  electrons  similar  to  that 
discussed  above  in  relation  to  electron  collisions.  Since  the  latter  can  be  associated 
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with  essentially  the  same  bounds  as  the  electron  scattering  case,  we  now  consider 
the  former. 

During  the  inverse-bremsstrahlung  process  a  photon  (hu)  is  absorbed  with  the 
resulting  energy  appearing  in  the  kinetic  energy  of  the  electron-ion  pair.  Assuming 
the  ion  to  be  structureless,  this  is  a  form  of  elastic  electron  scattering  for  which 
the  resulting  recoil  can  be  approximately  estimated  from  Eq.  (19)  provided  that  it 
is  possible  to  determine  appropriate  values  for  the  electron  energy  E  and  the 
impact  parameter  b.  In  terms  of  the  particle  energy,  we  can  associate  the  average 

energy  of  driven  motion  with  the  ponderomotive  potential.  For  radiation  at  248 

nm,  an  intensity  of  1016  W/cm2  gives  an  equivalent  ponderomotive  potential  of  57 
eV 

Measurements47  of  ion  energy  distributions  for  Na2*  ions  give  an  upper  bound 
on  Ej  =  60  meV  and,  therefore,  a  lower  limit  on  the  effective  value  of  b.  For 
Na2*  irradiated  at  248  nm  for  ~  600  ps  at  an  intensity  of  -  10l*  W/cm2,  with 

=  1,  Za  =  2,  E  =  57  eV,  and  M^  =  me,  Eq.  (19)  gives  b  >  0.1  a0.  This  value  of  the 

impact  parameter  is  expected  to  be  a  considerable  underestimate,  since  the  experimental 
bound  on  Et  was  instrumentally  limited.47  We  now  scale,  for  a  fixed  value  of  this 

impact  parameter,  to  fully  ionized  uranium  (Za  =  92)  and  let  I  =  1021  W/cm2  (em 

=  170),  a  procedure  which  gives  Et  <  1  meV.  Even  accounting  for  multiple  electron 
emission,  most  of  which  occurs  for  rather  low  values  of  Za,  it  is  unlikely  that 
the  recoil  will  be  above  the  experimental  bound  established  for  Na2\  namely, 
within  a  factor  of  two  of  the  thermal  energy  for  T  =  300  K.  This  places  an  upper 
bound  on  the  kinetic  energy  transfer  associated  with  the  strong-field  coupling  to 
be  on  the  order  of  those  arising  in  the  case  of  electron  scattering. 

We  observe  that  any  ion  will  also  experience  a  ponderomotive  potential  in 

the  high  intensity  wave.  For  a  proton  irradiated  with  248  nm  radiation  at  an 
intensity  of  1021  W/cm2,  the  ponderomotive  potential  is  -  3.1  keV.  This  energy, 
however,  is  fully  returned14®’146  to  the  radiation  field  during  the  course  of  the 
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pulse  so  that  no  net  displacement  of  any  molecular  ions,  regardless  of  charge  or 
mass,  will  occur  as  a  result  of  the  ponderomotive  potential. 

This  analysis  leads  to  the  following  general  conclusions.  Ion  collisions  can 
present  strong  forces  (em  >>  1)  and  large  recoil  energies  while  electron  collisions 
generate  relatively  weak  forces  (em  $  1)  and  correspondingly  small  recoil  effects. 
Furthermore,  both  of  these  scattering  processes  are  ineffective  in  exciting  organized 
electronic  motions  along  specific  molecular  axes.  In  comparison,  since  the  short- 
pulse  strong-field  excitation  (1)  can  generate  a  high  force  strength  (em  >>  1),  (2) 
communicates  a  recoil  energy  which  is  very  small  in  relation  to  molecular  bond 
strengths,  and  (3)  readily  excites  electronic  motions  along  specific  molecular 

directions,  it  represents  a  unique  domain  of  interaction. 

These  respective  regions  of  interaction  can  be  conveniently  displayed  in  the 
em  -  Et  plane  shown  in  Fig.  (18).  In  this  picture,  the  strong-field  multiquantum 
coupling  (em  >  1)  clearly  represents  a  distinct  regime  of  interaction  that  is  totally 
unachievable  with  charged  particle  collisions  and  permits  strong  electronic  forces 
to  be  applied  to  a  system  whose  molecular  frame  can  remain  intact.  It  is  this 

situation,  in  combination  with  the  dynamical  features  illustrated  in  Fig.  (17).  that 
enables  the  strong-field  coupling  to  efficiently  produce  new  modes  of  excitation  in 
molecular  materials.  The  region  corresponding  to  the  experiments  on  Na,  as 

shown  in  Fig.  (18),  clearly  falls  in  this  new  zone.  It  is  this  feature  of  the  interaction 
which  can  account  for  the  otherwise  anomalous  appearance  of  the  20g  vacancies 

and  the  55.8  nm  emission  described  in  Section  II.C.2.  Although  Na  has  served  as 
a  specific  example,  such  processes  are  expected  to  occur  generally  for  ail  materials 
regardless  of  structure  or  composition.  Given  the  nature  of  this  form  of  excitation, 
it  is  particularly  well  suited  for  the  direct  generation  of  states  involving  a  very 
large  charge  asymmetry,  such  as  that  represented  in  Na  by  Eq.  (8). 
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0.  Strong-Field  Propagation 


A  fundamentally  new  regime  of  electromagnetic  propagation  is  expected  to 

1 42— 

arise  in  plasmas  for  subpicosecond  radiation  at  a  sufficiently  high  intensity 
1 53 

Among  the  important  new  phenomena  that  may  occur  is  a  mode  of  channeled 
propagation,  a  key  process  for  the  generation  of  coherent  radiation  in  the  x-ray 
region.14^  When  a  very  intense  pulse  of  radiation  enters  a  cold  medium,  multiphoton 
ionization  locally  produces  a  plasma,  and  the  following  qualitative  behavior  involving 
a  relativistic/charge-displacement  mechanism  can  develop.  The  conversion  of  electrons 
from  bound  to  free  states  by  the  ionization  induces  a  strong  reduction  in  the 
refractive  response  of  the  atomic  material.  The  resulting  plasma  exhibits  a  further 
decrease  in  index  from  the  free-electron  component  collaterally  produced.  If  the 
ambient  electron  density  n0  is  such  that  the  plasma  frequency  up  =  /4irn0eJ/m  is 
below  the  frequency  of  the  radiation  u,  an  interesting  mode  of  channeled  propagation 
could  develop.  Indeed,  calculations  indicate  that  the  presence  of  both  relativistic 
motion  and  charge-displacement  can  lead  to  stable  self-focusing,  as  discussed  in 
Appendix  E. 

The  influence  of  the  charge-displacement  can  be  seen  in  the  following  way. 
For  a  sufficiently  short  pulse  (~  100  fs),  the  ions  remain  spatially  fixed  while  the 
relatively  mobile  electrons  are  expelled  by  the  ponderomotive  potential  from  the 
central  high-intensity  region  of  the  beam.  A  state  of  equilibrium  can  then  be 
established  between  the  ponderomotive  and  the  electrostatic  force  densities  owing 
to  the  net  charge  displacement.  Since  the  electrons,  which  embody  a  negative 
contribution  to  the  index,  are  expelled,  an  on-axis  region  of  relatively  high 
refractive  index  is  formed  which  can  support  channeling  of  the  laser  beam. 

Preliminary  analysis  concerning  the  channeling  indicates  (1)  that  the  Compton 

intensity  (~  4.5  x  10**  W/cm*  at  X  =  248-nm)  establishes  the  scale  for  the  prominence 

of  this  phenomenon  and  (2)  that  the  frequency  scaling  strongly  favors  the  use  of 

1 49 

the  shortest  wavelength  available  at  high  power  A  channel  radius  of  several 
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Energy  Transfer  (eV) 


100 


Fig.  (18):  Display  of  tha  comparison  of  the  energy  transfer  versus  electric  field 

strength  for  charged  particle  collisions  and  strong-field  interactions, 
including  the  influence  of  Auger  events  and  Inverse-bremsstrahlung. 
Clearly,  strong-fleid  interactions  lie  in  a  region  totally  inaccessible 
by  the  coilisional  processes.  The  region  of  the  N,  experiments  discussed 
in  the  text  is  shown.  The  parameters  for  the  Ion  boundary  stemming 
from  Eq.  (22)  correspond  to  X  *  248  nm  and  the  other  data  conforming 
to  the  N,  case. 
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wavelengths  also  emerges  from  this  analysis.  Moreover,  recent  dynamical  calcula- 
159 

tions  support  the  conclusion  that  such  modes  of  propagation  can  be  formed  under 
stable  conditions.  A  description  of  the  stability  that  arises  from  the  interplay  of 
relativistic  and  charge-displacement  mechanisms  is  given  in  Appendix  E.  The  ability 
to  stably  propagate  high  intensity  radiation  in  this  manner  would  be  a  fundamental 
advance  in  the  control  and  concentration  of  power  in  materials,  the  key  question 
relevant  to  the  production  of  amplification  in  the  x-ray  range.  Furthermore,  it 
facilitates  the  spatial  control  necessary  to  produce  a  high  brightness  output. 

Recent  experiments  have  led  to  the  observation  of  this  mode  of  propagation, 
as  described  in  Appendix  F.  A  comprehensive  theoretical  treatment  is  given  in 
Appendix  G.  The  importance  of  these  results  on  the  scaling  of  x-ray  amplifiers  is 
profound.  This  influence  is  described  in  Appendix  H. 

IV.  SUMMARY  AND  CONCLUSIONS 

The  use  of  recently  developed  high  brightness  subpicosecond  light  sources 
for  the  study  of  the  interaction  of  matter  with  strong  electromagnetic  fields  is 
leading  to  the  observation  of  new  phenomena  associated  with  new  classes  of 
excited  states  of  matter  and  new  modes  of  electromagnetic  propagation.  One 
important  consequence  of  these  findings  is  the  generation  of  powerful  x-ray 
sources  associated  with  properly  controlled  plasma  environments. 

This  report  has  discussed  several  aspects  related  to  the  general  question  of 
the  controlled  deposition  of  energy  at  high  power  levels  in  materials,  the  central 
issue  in  the  development  of  bright  x-ray  sources.  The  common  factor  linking  the 
various  topic  is  the  ability  to  produce  extremely  high  levels  of  electronic  excitation 
in  a  manner  which  enables  the  system  to  remain  kineticallv  cold  for  the  time 
scale  of  the  interaction.  For  the  situations  considered,  this  time  interval  is 
generally  less  than  a  picosecond.  It  is  the  mastery  and  control  of  this  class  of 
physical  processes  that  will  inevitably  lead  to  a  new  generation  of  extremely 
bright  sources  in  the  x-ray  range. 
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The  energy  penetration  depth  characteristic  of  the  interaction  of  intense  subpicosecond 
( -  600  fs )  ultraviolet  ( 248  nm )  laser  radiation  with  solid  density  material  has  been 
experimentally  determined.  This  was  accomplished  by  using  a  series  of  ultraviolet 
transmitting  targets  consisting  of  a  fused  silica  (Si02)  substrate  coated  with  an  80-600  nm 
layer  of  MgF2.  The  measurement  of  He-like  and  H-like  Si  and  Mg  lines,  as  a  function 
of  MgF2  thickness,  enabled  the  determination  of  the  energy  penetration  depth.  It  was  found 
that  this  depth  falls  in  the  range  of  250-300  nm  for  a  laser  intensity  of  —  3  X  1016  W/ 
cm2.  Based  on  numerical  simulations,  it  is  estimated  that  solid  density  material  to  a  depth  of 
-  250  nm  is  heated  to  an  electron  temperature  of  -  500  eV. 


Recently  developed  subpicosecond  high  intensity  lasers 
are  enabling  the  generation  of  high-energy  density  plasmas 
which  can  be  used  both  to  test  our  understanding  of  the 
dynamics  of  highly  excited  solids  and  for  the  development 
of  ultrafast  x-ray  sources.1"3  The  use  of  subpicosecond  laser 
pulses  to  produce  such  dense  plasmas  significantly  simpli¬ 
fies  the  analysis  of  the  laser-plasma  interaction,  since  hy¬ 
drodynamic  expansion  of  the  plasma  during  the  laser  pulse 
is  relatively  small.  It  is  crucial  to  understand  the  electron 
thermal  conduction  process  under  these  conditions,  since  it 
governs  both  the  peak  temperature  achieved  during  the 
laser  pulse  and  the  subsequent  cooling  rate,  thus  control¬ 
ling  the  intensity  and  duration  of  the  x-ray  pulse  produced. 

The  theoretical  uncertainties  involving  electron  heat 
transport  arise  mainly  due  to  the  presence  of  a  steep  ther¬ 
mal  gradient  (  -eV/A).  The  classical  treatment  of  thermal 
conduction  is  a  local  analysis/  wherein  the  scale  length  of 
the  temperature  gradient  is  assumed  to  be  much  larger 
than  the  electron  mean  free  path.  For  subpicosecond  laser 
plasmas,  this  condition  does  not  necessarily  hold,  so  a  non¬ 
local  treatment,5  in  which  hot  electrons  stream  into  cooler 
regions,  may  become  necessary.  In  the  absence  of  a  sophis¬ 
ticated  nonlocal  description,  hydrodynamic  simulation 
codes,  when  treating  gradients  with  a  scale  length  as  short 
as  a  mean  free  path,  simply  preserve  causality  by  limiting 
the  heat  flux  to  be  less  than  fnvT,  where  n  is  the  electron 
density,  u  the  electron  thermal  velocity,  and  T  the  electron 
temperature.  In  this  expression  /  is  the  "flux  limiter,"  a 
free  parameter,6  usually  with  a  magnitude  of  -0.1 

In  this  letter  we  report  on  the  experimental  determi¬ 
nation  of  the  heat  penetration  depth  of  the  laser-plasma 
interaction.  This  is  accomplished  by  monitoring  the  abso¬ 
lute  intensities  of  the  characteristic  line  radiation  emitted 
from  a  senes  of  ultraviolet  transmitting  targets  consisting 


*  Also  associated  with  F  M  Technologies,  Fairfax.  VA  22032 


of  a  fused  silica  ( Si02 )  substrate  coated  with  MgF:  layers 
of  known  thickness.  Measurements  of  the  relative  emis¬ 
sions  radiated  by  highly  ionized  Si  and  Mg  atoms  are  then 
used  to  establish  the  depth  of  energy  penetration.  Since  a 
subpicosecond  KrF*  laser  was  used  to  irradiate  the  targets 
in  this  study,  it  is  important  that  both  constituents  of  the 
layered  targets  are  sufficiently  transparent  to  low  levels  of 
the  ultraviolet  (248  nm)  radiation  that  no  prepulse  plasma 
ic  formed.7 

The  laser  system8  used  in  the  present  study  ( KrF*/248 
nm)  delivered  a  subpicosecond  pulse  (  —600  fs)  having  an 
energy  of  135  ±  15  mJ.  An  accompanying  prepulse  ( ASE ) 
had  a  duration  of  -  20  ns,  an  energy  of  8-30  mJ.  and  was 
focusable  to  a  spot  with  an  area  of  —  8x  10*  2  cm2  The 
effect  of  the  prepulse  has  been  carefully  considered.  Sub¬ 
strates  exposed  to  the  ASE  show  no  evidence,  under  mi¬ 
croscopic  examination,  of  a  tightly  focusable  (  -  10  jam 
diameter)  component.  Only  the  larger  (—8x10"'  cm’) 
pattern  can  be  made  visible  on  suitably  absorbing  materi¬ 
als.  Also,  theoretical  estimates  of  the  maximum  focusable 
intensity  arising  from  the  ASE  are  found  to  be  well  below 
damage  thresholds  for  the  target  materials  used.  All  the 
available  experimental  evidence7  and  analysis  indicate  that 
no  prepulse  plasma  was  formed  in  these  experiments 

The  subpicosecond  radiation  was  focused  onto  a  flat 
rotation  target  using  an  // 10  CaF2  lens.  The  peak  laser 
intensity  achievable  on  target  arising  from  the  subpicosec 
ond  pulse  was  estimated  to  be  — 1017  W/cm2.  This  esti¬ 
mated  intensity  is  based  on  direct  imaging  of  the  ultrav  idet 
focal  zone  with  a  microscope  and  CCD  camera,  a  measure 
ment  which  indicated  a  focal  spot  diameter  of  -  10  (im  I  n 
order  to  perform  this  measurement,  the  seed  beam  was 
focused  onto  a  surrogate  target  composed  of  a  250-um- 
thick  planar  pyrex  substrate.  This  target,  which  complete!  v 
absorbs  the  ultraviolet  radiation  in  a  scale  length  less  than 
a  micrometer,  produces  fluorescence  from  the  irradiated 
region  which  is  subsequently  imaged  with  an  optical  ~ 
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croscope  having  a  resolution  of  —0.5  /am. 

X-ray  pinhole  camera  measurements4  were  used  to  ex¬ 
amine  the  focal  region  at  full  beam  power.  For  these  mea¬ 
surements  a  pinhole  camera  used  DEF  film  and  a  filter 
composed  of  two  1000  A  layers  at  A1  on  a  2-/im-thick 
polycarbonate  substrate.  The  camera  had  a  5  jum  pinhole 
and  yielded  a  magnification  of  8.2 ±0.2  on  the  recording 
film.  The  x-ray  emission  was  measured  by  viewing  at  an 
angle  of  45°  with  respect  to  the  target  normal.  The  target 
used  for  this  examination  of  the  focal  spot  was  SiO-  coated 
with  250  nm  of  MgF:.  The  pinhole  camera  recorded  a 
radiating  plasma  with  a  diameter  in  the  range  of  8-9  pm 
for  radiation  in  the  kilovolt  spectral  region.  Furthermore, 
subsequent  microscopic  examination  of  the  SiO;  target  sur¬ 
face  revealed  6-7-^m-diam  craters.  These  findings  are  in 
reasonable  correspondence  with  the  measurement  of  the 
focusability  of  the  laser  seed  beam  described  above. 

In  these  experiments,  the  emissions  of  H-like  ( \s-lp) 
La  and  He-like  (lr-ls3p)  He/5  lines,  from  both  Si  and 
Mg,  were  monitored  as  a  function  of  both  incident  laser 
intensity  and  target  MgF:  layer  thickness.  Using  the  imag¬ 
ing  technique  described  above  to  monitor  the  focal  dimen¬ 
sions.  several  positions  of  the  focusing  lens  were  selected  in 
order  to  vary  the  incident  peak  laser  intensity.  The  x-ray 
emission  from  the  targets  was  collected  by  flat  KAP 
(2d  =  26.6  A)  and  ADP  (2^=10.64  A)  crystal  spec¬ 
trometers  using  DEF  film.  The  spectrometers  contained 
filters  composed  of  two  1000  A  layers  of  Al  on  a  2-pm- 
thick  polycarbonate  substrate  for  shielding  the  DEF  film 
from  unwanted  radiation.  Each  spectrum  was  produced  by 
25  laser  shots  and  the  absolute  line  intensities  were  de¬ 
duced  from  published  data  on  the  spectrograph-crystal  ef¬ 
ficiency,10  th*  filter  transmission,  and  the  film  calibration." 

The  experiments  discussed  below  were  modeled  using 
the  LASNEX  simulation  code.  This  analysis  involves  the 
solution  of  a  system  of  nonlocal  thermodynamic  equilib¬ 
rium  rate  equations,  which  incorporates  a  hydrogenic 
atomic  model  for  each  of  the  elements  (Mg.Si.O.F),  along 
with  flux-limited  thermal  conduction  (  /  =  0. 1 )  and  La- 
grangian  hydrodynamics.12  In  particular,  the  process  of 
thermal  conduction  was  treated  in  the  manner  described  by 
Lee  and  More.13 

The  measured  absolute  line  intensities  of  the  Si  and 
Mg  La  and  He/3  lines  from  bare  (uncoated)  SiO:  and 
MgF:  targets,  as  a  function  of  incident  laser  intensity,  are 
shown  in  Fig  1  along  with  the  corresponding  LASNEX 
simulations.  These  measurements  were  performed  in  order 
to  test  how  well  the  LASNEX  code  represented  the  inter¬ 
action.  Assuming  that  30%  of  the  incident  laser  energy  is 
absorbed  by  the  target,  the  LASNEX  calculation  agrees 
well  with  the  experimental  data  up  to  an  irradiance  of 
-2-3x  1014  W/cm2.  A  deviation  between  the  theoretical 
results  and  the  experimental  data,  however,  is  indicated  at 
the  highest  laser  intensity  of  -  1017  W/cm2.  For  compar¬ 
ison,  an  absorption  of  10%  was  assumed  for  laser  intensi¬ 
ties  above  -  1.5  X  1014  W/cm2.  Support  for  the  use  of  this 
value  comes  from  the  recent  work  at  high  intensity  by 
Mumane  et  al?  These  calculations  produce  a  reasonable 
match  of  the  data  corresponding  to  the  maximum  intensity 


FIG.  1  Absolute  line  intensities  of  Si  [La  C  and  He#  A)  and  Mg  •  La 
-  and  He#  x  )  radiation  measured  as  a  function  of  incident  laser  in¬ 
tensity  using  bare  SiO,  and  MgF;  targets.  The  curves  represent  LASNEX 
calculations  assuming  an  energy  absorption  of  30%  [ —  (solid!  Si  La. 

(closed  dash-dot)  Si  He#.  -  ■  -  (open  dash-dot)  Mg  La. - 

I  dash)  Mg  He#)  and.  beginning  at  an  intensity  of  I  5x  10'”  W/cm-.  an 
energy  absorption  of  10%.  shown  by  the  corresponding  extensions  with 
dotted  lines. 

(  -  1017  W/cm2)  used  in  these  studies.  Therefore,  for  the 
range  of  intensity  examined,  an  absorption  of  laser  energy 
in  the  10-30%  range  is  indicated.  A  variation  of  this  mag¬ 
nitude  could  arise  from  an  absorption  dynamically  depen¬ 
dent  on  the  incident  laser  intensity.  Further  studies  over  a 
larger  range  of  intensity  will  be  needed  to  more  fully  un¬ 
derstand  the  nature  of  this  dependence.  Nevertheless,  the 
LASNEX  simulation  is  found  to  give  a  reasonable  repre¬ 
sentation  of  the  experimental  results  for  the  conditions  of 
this  work. 

The  energy  penetration  depth  was  determined  by  mon¬ 
itoring  the  change  of  the  absolute  intensity  of  the  La  line 
emitted  by  H-like  Si  as  a  function  of  MgF2  coating  thick¬ 
ness.  Figure  2  presents  the  absolute  line  intensities  of  both 
Si  and  Mg  La  and  He/3  lines  as  a  function  of  MgF:  coating 
thickness.  The  measurements  were  carried  out  for  a  peak 
incident  laser  intensity  of  —  3x  1014  W/cm2.  an  intensity 
for  which  the  LASNEX  simulation  seems  adequate  As 
seen,  the  Si  La  line  intensity  decreases  rather  sharply,  by  a 
factor  of  approximately  e  from  the  bare  target  interaction, 
for  a  coating  thickness  in  the  range  of  250-300  nm.  indi¬ 
cating  that  the  energy  penetrates  to  a  depth  of  this  scale 
Qualitatively  similar  behavior  is  exhibited  by  the  Si  He/3 
transition.  Moreover,  both  the  Mg  La  and  He/3  lines  show 
the  expected  corresponding  increase  in  the  same  range  of 
thickness. 

The  results  stemming  from  the  LASNEX  simulations 
for  these  four  lines  are  also  shown  in  Fig.  2.  For  the  Si  La 
and  He/3  transitions,  it  is  seen  that  the  slope  representing 
the  rate  of  decrease  of  the  line  intensities  is  roughly  within 
a  factor  of  2  of  the  experimentally  observed  behavior  in  the 
region  where  most  of  the  energy  is  deposited  (MgF-  thick¬ 
ness  <  200  nm).  A  larger  deviation  occurs  in  the  region  in 
which  the  signal  is  relatively  small  (thickness  >  200  nm  t 
This  difference  at  the  larger  depths  indicates  that  some 
corrections  in  the  energy  transport  including  the  possibility 
of  nonlocal  heat  conduction,  are  indicated,  but  that  these 
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FIG.  2.  Absolute  line  intensities  of  Si  and  Mg  La  and  He/?  radiation  as  a 
function  of  MgF>  thickness  on  SiO;  substrates.  The  measured  data  are 
represented  as  (b )  Si  La-,  (A)  Si  He/?;  (  -»  )  Mg  La;  (  x  )  Mg  He0 
Solid  curves  are  drawn  through  the  data  to  guide  the  eye.  The  data  points 
indicated  at  1000  nm  of  MgF,  are  for  bare  MgF:  targets.  The  correspond¬ 
ing  LASNEX  calculations  assuming  i0%  absorption  and  an  incident 
laser  intensity  of  -  3  <  I0‘“  W/cttr,  are  also  shown  with  designations 

given  by  -  -  (dash-dot)  St  La. - (dash)  Si  He/?,  ...  (small  dot) 

Mg  La.  ■  ■  ■  ( large  dot )  Mg  He/? 

effects  are  relatively  small  and  do  not  involve  a  substantial 
fraction  of  the  energy.  For  the  Mg  La  and  He0  lines,  the 
LASNEX  simulations  are  seen  to  agree  quite  well  with  the 
measured  data  over  the  full  range  of  conditions  studied. 

While  the  data  presented  in  Fig.  2  indicate  that  the 
heat  wave  penetrates  to  a  depth  of  -  250  nm,  these  data  do 
not  directly  reveal  the  electron  temperature  to  which  the 
material  is  heated.  In  order  to  estimate  the  peak  tempera¬ 
ture  achieved  during  the  excitation,  the  electron  tempera¬ 
ture  was  calculated  as  a  function  of  the  incident  laser  in¬ 
tensity  by  the  LASNEX  code.  The  results  of  these 
calculations  are  shown  in  Fig.  3.  For  the  laser  intensity 
used  in  the  determination  of  the  energy  penetration  depth 
-  3  X  1016  W/cm;,  the  results  of  this  analysis  indicate  that 
the  electron  temperature  reaches  a  value  of  —500  eV  to  a 
depth  of  —250  nm  inside  the  solid  target  material. 

Spectral  studies  of  layered  SiCL/MgFj  targets  reveal 
that  an  energy  penetration  depth  of  250-300  nm  is  pro¬ 
duced  with  subpicosecond  excitation  at  248  nm  with  an 
intensity  of  -  3  x  1016  W/cm2  This  value  is  considerably 
in  excess  of  the  skin  depth  associated  with  the  plasma  being 
produced.  LASNEX  modeling  of  the  interaction  was  found 
to  be  in  good  agreement  over  the  range  of  conditions  stud¬ 
ied.  An  estimate  of  the  electron  temperature  of  the  heated 
material,  as  indicated  by  the  LASNEX  simulation,  gives  a 
value  of  -  500  eV  at  a  depth  of  -  250  nm  inside  the  solid 
material. 


FIG.  3  LASNEX  calculations  of  electron  temperature  as  a  function  of 
peak  laser  irradiance.  Calculations  are  for  a  peak  electron  density  of 
6  -  10"  cm  "  at  the  time  of  peak  laser  irradiance.  with  the  assumption  of 
30^  absorption  for  all  irradiances. 
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Studies  of  the  interaction  of  condensed  matter  with  short  pulse  (  ~  600  fs)  high  intensity 
(  -  101  W/cm:)  ultraviolet  (248  nm)  radiation  show  that  intense  spatially  compact 
sources  for  x-ray  emission  in  the  kilovolt  range  (0.5-1. Ox  10'5  W/cm2)  can  be  generated 
from  solid  targets  at  close  to  the  maximum  volume  specific  rate  allowed. 


Solid  materials  offer  the  possibility  of  combining  a  high 
particle  density  with  a  very  high  level  of  electronic  excita¬ 
tion,  provided  that  the  deposition  of  energy  in  the  solid  can 
be  made  to  occur  sufficiently  rapidly.  Ultraviolet  high- 
brightness  subpicosecond  light  sources1*7  enable  such  rapid 
deposition  of  energy  to  be  achieved. 

The  laser  system5  used  in  the  present  study  (KrFV248 
nm)  delivered  a  subpicosecond  pulse  (  -  600  fs)  having  an 
energy  of  135  ±  15  mJ.  An  accompanying  prepulse  had  a 
duration  of  —  20  ns,  an  energy  of  8-30  mJ,  and  was  focus- 
able  to  a  spot  with  an  area  of  -  8  X  10  ~ 2  cm2.  The  subpi¬ 
cosecond  radiation  was  focused  onto  a  flat  rotating  BaF2 
target  using  an  f/ 10  CaF2  lens.  The  peak  laser  intensity 
achievable  on  target  arising  from  the  subpicosecond  pulse 
was  estimated  to  be  -  1017  W/cm2.  This  estimated  inten¬ 
sity  is  based  on  direct  imaging  of  the  ultraviolet  focal  zone 
with  a  microscope  and  a  charge  coupled  device  (CCD) 
camera,  a  measurement  which  indicated  a  focal  spot  diam¬ 
eter  of  6-7  ^m. 

Independent  x-ray  pinhole  camera  measurements  were 
also  used  to  examine  the  focal  region.  For  these  measure¬ 
ments  a  pinhole  camera  using  Kodak  DEF  film  and  a  filter 
of  2000-A-thick  A1  on  a  2-^m-thick  polycarbonate  sub¬ 
strate  was  used.  The  camera  had  a  5  fsm  pinhole  and 
yielded  a  magnification  of  8.2  ±0.2  on  the  recording  film. 
Viewing  at  an  angle  of  45*  with  respect  to  the  target  nor¬ 
mal.  the  x-ray  emission  from  two  types  of  solid  targets  was 
measured.  The  targets  were  Si02  coated  with  250  nm  of 
MgF:  and  bare  BaF:.  The  former  indicated  a  radiating 
plasma  with  a  diameter  in  the  range  of  8-9  fim  for  radia¬ 
tion  in  the  kilovolt  spectral  region.  Microscopic  examina¬ 
tion  of  the  Si02  target  surface  revealed  6-7  ^m  diameter 
craters,  a  finding  in  reasonable  correspondence  with  the 
measurement  of  the  focusability  of  the  laser  beam  de¬ 
scribed  above.  The  measurements  with  the  BaF2  targets 
indicated  a  radiating  plasma  with  a  diameter  of  -20  ±4 

and  corresponding  craters  with  diameters  in  the  8-10 
(urn  range  were  observed  in  the  target  substrate. 

For  spectral  studies  of  BaF2  ,  the  x-ray  em:  .  on  from 
the  targets  was  collected  by  a  flat  KAP  (2 d  =  26.6  A) 
crystal  spectrometer  which  produced  a  recording  on  DEF 
film.  The  absolute  line  intensities  were  deduced  from  ex¬ 
tant  data  on  the  crystal  efficiency,5  the  filter  transmission. 


1  Also  of  FM  Technologies.  Fairfax.  V  A  22032 


and  the  DEF  film  calibration.9  To  insure  the  reproducibil¬ 
ity  of  the  data,  the  spectra  were  collected  several  times  over 
a  period  of  six  months.  The  film  was  developed  according 
to  established  procedures  for  the  calibration  of  the  film  and 
the  optical  density  was  digitized  at  the  facility  at  the  Naval 
Research  Laboratory.  Each  spectrum  was  produced  by  an 
exposure  accumulated  over  25  laser  shots. 

An  experimental  concern  in  studies  of  the  interaction 
of  intense  radiation  with  solid  matter  is  the  effect  of  any 
low  intensity  prepulse  on  the  true  state  of  the  target  mate¬ 
rial.  If  a  sufficient  intensity  is  present  in  the  prepulse,  a 
plasma  can  be  formed  which  will  shield  the  solid  target 
from  the  high  intensity  subpicosecond  energy.  Experimen¬ 
tal  studies,  however,  have  demonstrated  that  the  effect  of 
the  prepulse  can  be  eliminated  at  248  nm  by  using  a  target 
composed  of  ultraviolet  transmitting  materials.10  For  ex¬ 
ample,  MgF2  and  BaF2  exhibit  essentially  no  absorption  " 
up  to  a  relatively  high  intensity  ( =*  1011  W/cm2).  There¬ 
fore,  if  the  prepulse  intensity  is  appreciably  less  than  10n 
W/cm2,  negligible  absorption  is  present  and  no  interfering 
plasma  is  generated.  In  this  case,  the  threshold  intensity  for 
nonlinear  absorption  is  sufficiently  high  to  prevent  any  sig¬ 
nificant  influence  by  the  prepuise  in  these  studies,  since  the 
maximum  prepulse  intensity  occurring  was  only  —  2.5 
x  107  W/cm2. 

The  spectrum  observed  from  BaF2  targets  exhibited 
very  strong  emission  in  the  9-13  A  range,  as  shown  in  Fig. 
1 .  The  radiation  in  the  broad  band  centered  at  -12  A 
arises  from  Af-shell  transitions  of  barium  ions  and  the  ob¬ 
served  structure  can  be  put  in  reasonable  correspondence 
with  the  previously  measured12  3d"— 3d"~  1  4/  sequence  in 
that  material.  ”'ith  consideration  of  the  response  of  the 
DEF  film,  the  crystal  response,  and  the  filter  calibrations, 
these  measurements  indicated  that  —  2-4  mJ  of  x  rays  were 
produced  in  this  spectral  region  on  each  shot,  the  greater 
production  corresponding  to  the  higher  range  of  the  inci¬ 
dent  laser  energy.  The  isotropy  of  the  emission  was  estab¬ 
lished  experimentally  by  recording  the  spectra  at  two 
widely  different  angles  (45*  and  80*)  with  respect  to  the 
target  normal.  In  both  cases,  the  properties  of  the  recorded 
spectra  were  similar.  Although  the  duration  of  the  emis¬ 
sion  was  not  directly  measured,  an  pper  bound  of  1-2  ps 
is  consistent  with  ( 1 )  the  size  of  the  radiating  region  mea¬ 
sured  with  the  pinhole  camera  for  a  velocity  of  expansion 
of  a  few  Angstroms  per  femtosecond.  (2)  with  LASNEX 
calculations  for  our  experimental  conditions,1'  (3)  with 
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FIG  I  Spectrum  of  x-ray  emission  observed  from  BaF.  at  a  laser  inten¬ 
sity  of  -  10  1  'V/cm; 

the  time  resolved  studies  of  laser-produced  plasmas  re¬ 
ported  by  Mumame  et  a/.,14  and  (4)  recent  measurements 
of  Cobble  et  al.  with  solid  aluminum  targets.15 16  The  latter 
explicitly  indicates  rapid  cooling  that  would  cause  the 
emission  of  kilovolt  radiation  to  be  of  short  duration.  It 
follows,  assuming  isotropic  emission  from  the  radiating 
zone  determined  by  the  pinhole  camera,  that  for  BaF:,  an 
x-ray  intensity  in  the  range  of  5x  1014  </t<1015  W/cm2 
was  produced  on  each  laser  shot. 

An  upper  bound  on  the  x-ray  intensity  achievable  from 
this  type  of  source  can  be  estimated  by  considering  a  slab  of 
unit  area  with  density  p  and  thickness  5  which,  upon  ex¬ 
citation  with  intense  subpicosecond  radiation,  produces  ra¬ 
diation  in  the  kilovolt  spectral  range.  Since  the  time  of 
excitation  is  short,  in  first  approximation  we  will  ignore  the 
hydrodynamical  expansion,  a  process  which  will  occur 
roughly  at  a  rate  of  =1-10  A/fs.  An  upper  limit  on  the 
x-ray  intensity  Ixm  can  then  be  written  as 

[x„=,fuoxp6/2rx  (1) 

in  which  fuox  represents  a  kilovolt  quantum  and  rx  is  the 
effective  radiative  rate  of  the  excited  material.  For  the 
present  discussion,  we  will  take  p  as  BaF,  solid  density 
(  —  1.7  x  10'2  cm  "’)  and  use  for  r,  a  typical  radiative  rate 
for  kilovolt  transitions  in  atoms  or  ions.  The  latter  param¬ 
eter  will  be  taken  in  the  range1’  of  2  X  10  “  ,4<r<  s;  10 ' 1 5  s. 

It  remains  to  estimate  the  thickness  of  the  radiating 
zone  denoted  by  5.  Using  a  senes  of  MgF2/SiO,  targets 
coated  with  a  range  of  thicknesses  of  MgF2,  the  measured18 
depth  of  energy  penetration,  for  plasmas  generated  at  an 
incident  laser  intensity  of  —  3x  I016  W/cm2,  was  found  to 
be  in  the  range  of  250-300  nm.  The  depth  8  was  deduced 
from  the  observed  relative  strengths  of  the  Mg  and  Si 
H-like  and  He-like  transitions  as  a  function  of  the  thick¬ 
ness  of  the  MgF-  layer  Taking  this  energy  penetration 
depth  for  8.  the  approximate  upper  bound  represented  by 
Eq.  ( 1 )  can  be  evaluated  From  the  given  ranges  on  r,  and 
<5,  and  assuming  that  the  radiating  material  is  thin  to  the  x 
rays,  we  find  3  x  1014  y  /tm  y  2  x  10i?  W/cm"  The  avail¬ 
able  information  on  opacities19  indicates  that  the  assump¬ 
tion  of  a  thin  plasma  is  consistent  with  the  range  of  8  given 

Comparison  of  this  upper  limit  on  /[m  with  the  exper¬ 
imental  value  /,  indicates  that  the  observed  intensity  falls 
within  the  estimated  range  for  maximum  possible  emission 


Although  such  a  source  is  incoherent,  a  radiation  rate  of 
this  magnitude,  coupled  with  the  relatively  narrow  spec¬ 
trum  over  which  it  occurs,  constitutes  a  source  of  ex¬ 
tremely  high  spectral  brightness  that  can  produce  unusu¬ 
ally  high  excitation  rates  in  absorbing  material  to  which  it 
can  be  coupled.  We  note  that  since  the  simple  model  used 
to  estimate  the  upper  radiative  bound  is  far  from  an  equi¬ 
librium  picture  and  since  the  experimental  result  is  quite 
close  to  that  limit,  the  implication  is  that  the  plasma  state 
may  deviate  significantly  from  equilibrium  for  this  type  of 
intense  and  rapid  excitation. 

In  conclusion  it  is  found  that  the  use  of  high  brightness 
ultraviolet  subpicosecond  light  source  for  the  study  of  the 
interaction  of  condensed  matter  with  strong  electromag¬ 
netic  fields  is  leading  to  the  development  of  ( 1 )  the  ability 
to  couple  a  strong  laser  field  directly  to  solid  target  mate¬ 
rial,  without  the  production  of  a  prepulse  plasma,  and  (2) 
the  generation  of  powerful  x-ray  sources  associated  with 
spatially  compact  high  density  plasmas. 
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Plasma  produced  from  ultraviolet-transmitting  solid  targets  undergoing  intense  (>10‘*  W/cm2)  subpicosecond 
(—600  fs)  ultraviolet  (248  nm)  irradiation  have  been  studied  under  conditions  for  which  no  interfering  prepulse 
plasma  is  generated.  Time  and  spatially  integrated  measurements  of  the  x-ray  emission  for  H-like  and  He-like 
Mg  and  Si  were  found  to  be  consistent  with  lasnex  calculations  that  model  the  laser- target  interaction. 


Subpicosecond  high-peak-brightness  lasers  capable 
of  delivering  focal  intensities  of  greater  than 
10 16  W/cm2  permit  the  generation  of  high- temperature 
plasmas  from  solid  targets.  These  plasmas  can  pro¬ 
duce  short  pulses  of  x  rays,  of  the  order  of  a  picosec¬ 
ond  or  less,  owing  to  the  rapid  cooling  of  the  plasma 
on  the  solid  target  surface. 

Numerous  experiments  involving  plasma  produc¬ 
tion  with  picosecond  and  subpicosecond  lasers  at 
different  wavelengths  and  intensities  have  been  re¬ 
ported.1'6  An  important  experimental  concern  in 
subpicosecond  studies  involving  the  coupling  of  laser 
energy  to  solid  targets  is  the  potential  influence  on 
the  target  of  any  relatively  low-intensity,  but  long¬ 
time-scale  prepulse  that  may  be  present,  since  such 
prepulses  are  generally  associated  with  amplified 
spontaneous  emission  (ASE)  in  high-power  laser  sys¬ 
tems.  Moreover,  since  the  ablation  threshold  for 
most  solid  materials  is  less  than  -108  W/cm2,  even  a 
relatively  low  level  of  ASE  originating  from  the  laser 
amplifier  may  be  sufficient  to  ablate  and  ionize  the 
target  material,  creating  a  long-scale-length  plasma 
and  thus  altering  the  dynamics  of  interaction  of  the 
intense  subpicosecond  pulse  with  the  solid  target. 
Consequently,  the  relatively  low-density  plasma  that 
may  be  generated  by  the  prepulse  can  severely  alter 
the  interaction  of  the  subpicosecond  energy  with  the 
solid  matter. 

In  this  Letter  we  report  experimental  findings 
showing  that  the  effect  of  the  prepulse  can  be  elimi¬ 
nated  when  a  KrF*  laser  system  is  used  at  248  nm 
to  irradiate  optically  polished  ultraviolet-transmit¬ 
ting  targets  such  as  CaF2,  MgF2,  BaF2,  and  Si02.  In 


addition,  we  report  on  the  results  of  spectrally  re¬ 
solved  x-ray  emission  from  CaF2,  MgF2,  and  Si02 
targets  irradiated  at  the  peak  intensity  of  =3  x 
101®  W/cm2.  Based  on  the  properties  of  the  observed 
line  emissions,  the  characteristics  of  the  plasmas 
generated  in  these  studies  were  compared  with 
lasnex7  simulations. 

The  KrF*  laser  system  used  for  these  experiments 
has  been  described  elsewhere.8  It  operated  with  a 
pulse  duration  of  =600  fs,  had  a  subpicosecond  pulse 
energy  of  -90  mJ  accompanied  with  =20  mJ  of  ASE 
that  had  a  duration  of  —20  ns,  and  was  focused  onto 
a  rotating  planar  target  with  an  fj  10  CaF2  lens 
The  accompanying  prepulse  (ASE)  was  focusable  to 
a  spot  with  an  area  of  —8  x  10~2  cm2  owing  to  the 
relatively  large  divergence. 

The  effect  of  the  prepulse  has  been  carefully  con¬ 
sidered.  Substrates  exposed  to  the  ASE  show  r.o 
evidence,  under  microscopic  examination,  of  a 
tightly  focusable  (~10-/im-diameter)  component 
Only  the  larger  (-8  x  10~2  cm2)  pattern  can  be 
made  visible  on  suitably  absorbing  materials.  These 
tests  were  performed  with  the  beam  blocked  at  the 
input  to  the  amplifier  chain.  Also,  theoretical  esti- 
mates  of  the  maximum  focusable  intensity  that 
arises  from  the  ASE  are  found  to  be  well  below  dam 
age  thresholds9,10  for  the  target  materials  used 
Since  the  mechanisms  of  absorption  in  these  materi 
als  at  248  nm  are  a  multiphoton  process,  the  mm: 
mum  intensity  required  to  induce  significant 
absorption910  is  in  the  range  of  1010-10n  W/cm:  \ . 
the  available  experimental  evidence  and  analysis  r. 
dicate  that  no  prepulse  plasma  was  formed  in  :!>•»• 
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Fig.  1.  X-ray  spectrum  from  a  Si02  target  showing  tran¬ 
sitions  from  highly  ionized  Si  ions.  The  main  features 
originating  from  He-like  Si  xili  and  H-like  Si  Xiv  ions  are 
identified.  Sats.  indicates  the  position  of  satellite  lines 
identified  with  Si  xii  ions.  The  horizontal  bars  indicated 
on  the  spectrum  corresponding  to  the  He-a,  L -a,  He-/3, 
He-y,  and  He-d  transitions  designate  the  lasnex  calcu¬ 
lated  intensities  of  these  lines.  These  calculated  results 
are  shown  normalized  to  the  He-a  feature. 

experiments.  Therefore  these  materials  and  the 
plasmas  produced  therefrom  should  not  be  affected 
by  the  laser  prepulse  present  in  these  studies. 

The  expected  absence  of  a  prepulse  plasma  was 
experimentally  verified.  This  was  achieved  by  us¬ 
ing  only  the  ASE  to  irradiate  comparatively  a  MgF2 
target  and  a  multilayer  target  (100-nm  Ni  on  a 
5-mm-thick  MgF2  substrate).  When  the  ASE  was 
focused  onto  the  bare  MgF2  target,  no  visible  break¬ 
down  was  observed,  and  no  surface  damage  was  evi¬ 
dent.  When  the  multilayer  Ni  target,  with  the  Ni 
side  facing  the  laser,  was  irradiated  with  the  ASE, 
breakdown  was  observed  as  well  as  surface  damage 
on  the  Ni.  After  reversal  of  the  target  so  that  the 
uncoated  MgF>  surface  faced  the  laser,  a  plasma  was 
again  observed  on  the  Ni  side  of  the  target,  with  no 
observable  plasma  formed  on  or  damage  to  the  front 
MgF2  surface.  In  this  latter  configuration  the  en¬ 
ergy  was  freely  transmitted  to  the  rear  surface. 
Therefore  it  is  concluded  that  the  use  of  properly  se¬ 
lected  ultraviolet-transmitting  materials  was  effec¬ 
tive  in  eliminating  the  unwanted  influence  of  the 
nanosecond  prepulse  present  in  this  research. 

The  x-ray  emission  from  all  targets  was  collected 
by  flat  potassium  hydrogen  phthalate  (2d  = 
2.66  nm)  and  pentaerythritol  (2d  =  0.87  nm)  crystal 
spectrometers  and  recorded  on  DEF  film.  The 
spectrometers  contained  filters  composed  of  two 
100-nm  layers  of  A1  on  a  2-um-thick  polycarbonate 
substrate  for  shielding  the  DEF  film  from  unwanted 
radiation.  Approximately  25  shots  were  accumu¬ 
lated  for  each  spectrum.  Portions  of  the  spectra  ob¬ 
tained  from  a  Si02  target  are  shown  in  Fig.  1.  Also 
shown  in  Fig.  1  are  the  calculated  line  intensities  for 
the  He-a,  L-a,  He-0,  He-y,  and  He-6  transitions 
from  ionized  Si.  The  ratio  of  line  strengths,  nor¬ 
malized  to  the  He-a  feature,  are  seen  to  agree 
within  approximately  20%.  For  MgF2  and  Si02 
targets,  both  He-like  and  H-like  emissions  were 
identified  from  Mg  and  Si  ions.  For  the  CaF2 


target,  only  the  ls2-ls2p  line  emitted  by  He-like 
Ca  ions  was  found.  The  absolute  line  intensities 
were  deduced  from  known  data  for  the  film  calibra¬ 
tion,11  the  filter,  and  the  response  of  the  crystal 
spectrometer.12 

The  lasnex  code  was  used  to  simulate  the 
experiment  at  3  x  1016  W/cm2,  assuming  a  30%  ab¬ 
sorption  by  the  target  material,  a  value  that  is  com¬ 
parable  with  the  absorption  measured  in  Ref  4 
The  treatment  incorporates  flux-limited  thermal 
conduction  (with  the  thermal  flux  limited  to  its 
free-streaming  value  times  a  coefficient  f,  taken  to 
be  0.1),  hydrodynamic  expansion,  and  a  nonlocal 
thermodynamic  equilibrium  solution7  to  the  atomic 
rate  equations  in  the  context  of  a  hydrogenic  atomic 
model.  The  choice  of  f  =  0.1  is  based  on  previous 
findings  of  its  suitability  for  Fokker-Planck  heat 
transport  simulations.13  Furthermore,  the  results 
are  not  sensitive  to  the  magnitude  of  the  flux  lim¬ 
iter  (/“);  the  use  of  f  =  0.2  in  the  calculations 
changes  the  outcome  by  less  than  20%.  Figure  2 
exhibits  the  calculated  development  of  both  the  elec¬ 
tron  temperature  Tt  and  the  electron  density  n,  for  a 
Mg  plasma  produced  by  a  KrF*  laser  intensity  of 
3  x  1016  W/cm2  at  the  time  of  peak  laser  irradiation 
[Fig.  2(a)]  and  1  ps  after  the  peak  laser  intensity 
[Fig.  2(b)].  As  shown,  a  high-density  hot  plasma  is 
formed  initially  at  the  surface.  After  rapid  expan¬ 
sion,  a  hot  subcritical  corona  is  produced.  Thermal 
conduction  inward  creates  a  fairly  hot,  solid-density 
region  of  nearly  fully  ionized  Mg  and  F.  Based  on 
the  simulated  density  and  temperature  profile,  the 
absolute  line  intensities  were  calculated  for  ls2-ls3p. 
Is2-ls2p  He-like,  and  ls-2p  H-like  transitions  of 
Mg  and  Si  ions.  These  calculations  were  compared 
with  the  experimental  values  obtained  from  the 
DEF  film  and  found  to  be  in  agreement  to  better 
than  a  factor  of  2.  Furthermore,  for  MgF2  targets 
subjected  to  an  intensity  of  =®3  x  1016  W/cm2,  the 
simulation  indicates  that  a  plasma  with  an  electron 
temperature  of  the  order  of  500  eV  at  a  solid  density 
is  produced. 

In  conclusion,  when  a  subpicosecond  KrF*  laser 
system  is  used  to  irradiate  a  high-damage-threshold 
ultraviolet-transmitting  target,  the  formation  of  a 
prepulse  plasma  is  prevented.  This  allows  the  high- 
intensity  laser  energy  to  create  solid-density  hot  ma¬ 
terial.  It  should  be  noted  that  this  technique 
substantially  relaxes  the  requirement  for  the  con¬ 
trast  ratio  of  the  ASE  and  pulse  intensities  for  high- 
power  short-pulse  laser  systems. 


Fig.  2.  lasnex  simulation  indicating  the  electron  terr.p-- 
ature  and  density  established  (a)  at  the  time  of  peak  ;  n  t  -  - 
sitv  of  irradiation  and  (b)  1  ps  after  the  peak  vaiue 
irradiation. 
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ABSTRACT 


Studies  of  the  interaction  of  matter  with  high  intensity  radiation  io;‘ 
W/cm3)  are  leading  to  the  observation  of  coulomb  explosions  from  highly  excited 
molecules  and  intense  x-ray  emission  from  solid  targets. 

I.  INTRODUCTION 

Advances  in  femtosecond  lasers  are  extending  the  exploration  of  multiphoton 
interactions  well  into  the  strong-field  regime,  the  physical  situation  in  which  the 
external  field  is  greater  than  an  atomic  unit  (e/a^).  The  performance  projected  for 
ultraviolet  rare  gas  halogen  technology  is  currently  being  realized1-10  and  new 
near-infrared  solid  state  systems,  such  as  Ti;AI20,.  are  under  vigorous  development 
Both  technologies11  should  reach  a  field  strength  of  -  iGO  e,a^  with  instruments 
that  produce  an  output  energy  of  -  U. 

There  is  a  growing  body  of  evidence  that  fundamentally  new  forms  of  hignly 
excited  matter  can  be  produced  with  strong-field  interactions.  Specifically,  there 
is  considerable  interest  in  (1)  mechanisms  leading  to  high  levels  of  electronic 
excitation,  particularly  those  involving  ionization  and  the  excitation  of  inner-shell 
states,  and  (2)  the  properties  of  matter  which  govern  the  generation  of  intense  short 
wavelength  radiation,1 

II.  EXPERIMENTAL  CONSIDERATIONS 

The  KrF*  (248  nm)  source  used  for  the  studies  discussed  herein  has  been 
described  elsewhere.®1®  This  instrument  currently  has  an  output  pulse  width  of  = 

*  Stab#  published  in  the  Proceedings  of  the  International  Conference  on 
Mulciphoton  Processes  (ICOMP  V,  Saciav,  Iy9l). 


85 


600  fs.  a  typical  pulse  energy  in  the  300-400  mj  range,  and  focusablllty  nominally 
within  a  factor  of  two  of  the  diffraction  limit.3  For  example,  with  a  simple  f/10 

CaF,  lens,  peak  intensities  of  -  3  x  io1*  W/cm1  are  produced  even  though  such  a 
lens  exhibits  appreciable  spherical  aberration.  The  application  of  more  sophisticated 
focusing  systems  can  c'early  generate  substantially  greater  field  strengths.3,9 

A  recent  modification  has  incorporated  a  Ti:AljO,  amplifier  into  the  system 
for  amplification  of  the  radiation  at  745  nm.  This  change  has  resulted  in  approximate¬ 
ly  a  40-fold  increase  n  the  248  nm  output  energy  from  the  mixing  crystal  to  a 

value  of  -  40  yj.  This  enables  the  generation  of  a  significantly  more  intense  seed 
beam  to  drive  the  power  amplifier  and  a  reduction  of  the  intensity  of  the  amplified 
spontaneous  emission. 

A  more  detailed  characterization  of  the  pulse  shape12  for  this  system  Is 

becoming  available  with  the  use  of  a  technique  Involving  a  third  order  intensity 

correlation13.  One  of  the  properties  is  a  clearly  observed  asymmetry  In  rhe 
profile  of  the  pulses,  in  our  case,  a  waveform  rising  more  rapidly  than  It  falls.  The 
current  seed  beam  generates  pulses  with  a  width  as  short  as  -  260  fs. 

The  ability  to  concentrate  the  pulse  energy  in  a  small  focal  zone  is  another 
crucial  property  necessary  for  the  generation  of  high  field  strengths.  With  an  f/7 
focusing  system,  focal  zones  with  a  diameter  of  -  6.6  ym  have  been  measured,  a 
dimension  which  corresponds  to  a  value  of  -  t  5— fold  the  diffraction  limit. 

III.  STUDIES  OF  MOLECULAR  ION  PRODUCTION 

Studies  of  collision-free  Ionization  have  been  a  valuable  source  of  Information 
on  multiphoton  processes.14-17  Although  most  studies  have  been  conducted  with 
atoms,  recent  work  has  begun  to  examine  molecules,  including  certain  polyatomic 
systems  Indeed,  studies  of  molecular  ionization  have  produced  valuable  new  Insights 
into  the  dynamics  of  multiphoton  processes.  Energetic  coulomb  explosions 

have  been  studied  in  which  (1)  the  fragments  have  kinetic  energies  of  several 
tens  of  eV  per  particle  (2)  evidence  for  molecular  inner-shell  excitation  has  been 
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observed.26  and  (3)  molecular  atomic  site-dependent  energy  deposition  has  been 

detected.27 

A.  Na 

It  is  oossible  to  measure  the  kinetic  energy  distributions  of  the  products 
arising  from  a  molecular  coulomb  explosion  with  an  appropriately  arranged  time— 
of— flight  apparatus  22  N^‘  Ionic  fragments  produced  with  a  peak  Intensity  of  - 

3  *  TO  W/cm3  are  observed  with  charges  in  the  range  1  «  q  £  5.  The  fragment 

distributions  generally  fall  off  at  sufficiently  high  ion  energy  and  commonly  exhibit 
observable  structure.  If  we  take  the  energy  at  which  the  distribution  falls  to  one 
half  its  maximum  value  as  an  approximate  measure  of  the  characteristic  kinetic 
energy  of  a  fragment,  we  can  summarize  the  findings  on  Na  as  shown  In  Fig.  (i). 

The  coulomb  energy  associated  with  the  production  of  N*‘  fragments  at  - 
38  eV  implies,  given  the  equilibrium  Internuciear28  separation  of  N,  as  r„  =  1.1  A, 

the  generation  of  a  nascent  ionized  molecule  of  the  form  N ,9*  with  q  >  5.  This 


Fig.  (1):  Summary  of  fragment  energies  observed  In  Na. 
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Indicates  that  during  the  course  of  the  Interaction  the  Ionization  must  proceed  quite 
rapidly  to  at  least  this  stage.  This  finding  Is  in  general  concurrence  with  conclusions 
on  the  rate  of  energy  transfer  stemming  from  earlier  studies.22  The  minimum  energy 
required  to  produce  N,s*.  assuming  an  initial  internuclear  separation  of  -  r0.  Is  - 
224  eV. 

8.  Comparison  of  CO  and  CO, 

The  comparison  of  the  carbon  and  oxygen  fragments  observed  from  CO  and 
CO,  are  providing  further  information  on  strong-field  coupling  to  molecules.  In 
considering  the  data  shown  below  for  C“*  fragments  in  Fig.  (2),  it  should  be 
understood  that  the  shapes  of  the  time— of— flight  Ion  signals  are  determined  by 
both  the  instrumental  resconse  and  the  kinetic  energies  of  the  fragments.  The 
individual  Ion  signals  are  actually  composed  of  two  components  which  become 
resolved  at  sufficiently  high  particle  energy.  One  component  arises  from  ions 
Initially  formed  with  velocities  parallel  to  the  tlme-of-flight  direction  while  the 


C4*  from  C02  C **"  from  CO 


3  :!1J  T  l.J  it  UTO 


Fig.  (2):  Ion  tlme-of-fMght  spectra  for  C'*:  (a)  Ions  produced  from  CO,  appear 

at  low  energy;  (b)  Ions  produced  from  CO  show  two  well  resolved  peaks 
with  a  maximum  kinetic  energy  -  40  eV. 
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other  arises  from  Ions  whose  initial  velocity  is  anti-parallel  to  that  vector.  The 

greater  the  separation  of  these  two  peaks,  the  greater  is  the  energy  associated 

with  the  kinetic  energy  of  the  fragments.  The  two  peaks  are  of  comparable 

height  and  width  for  appropriate  settings  of  the  extraction  and  acceleration 

20  22 

voltages  used  in  the  time— of— flight  apparatus.  ’ 

The  time— of— fl ight  data  for  C“  observed  from  CO,  are  shown  in  Fig.  (2a). 

In  this  case,  only  one  peak  Is  seen,  a  clear  sign  that  the  particles  are  slow. 

This  conforms  to  the  expectation  that  the  carbon  atom  in  the  center  of  the  molecule 

would  develop,  in  the  process  of  fragmentation,  a  relatively  small  kinetic  energy  as 

25  29 

a  result  of  the  opposing  forces  from  the  oxygen  atoms  on  either  side.  ' 

The  corresponding  data  for  C**  fragments  produced  from  CO  are  illustrated 
in  Fig.  (2b).  These  data  clearly  exhibit  two  well  resolved  peaks  Indicating  that  a 
substantial  kinetic  energy  is  present.  The  maximum  kinetic  energy  observed  in  this 
case  is  approximately  40  eV.  Basically.  In  the  comparison  of  CO,  and  CO.  the 
carbon  Ions  are  produced  with  a  low  kinetic  energy  from  the  former  and  a  high 
kinetic  energy  from  the  latter. 

The  characteristics  of  the  oxygen  Ions  produced  represent  a  rather  different 
pattern.  The  situation  can  be  illustrated  by  consideration  of  the  observed  kinetic 
energy  spectra  for  O3*  which  are  shown  In  Fig.  (3).  In  this  case  both  molecules 
produce  energetic  ions,  but  those  arising  from  CO,  are  significantly  more  energetic. 
Comparison  of  Fig.  (2b)  and  Fig.  (3b)  indicates  that  the  maximum  energies  of  the 
C“*  and  O1*  ions  are  roughly  comparable  when  produced  from  CO.  This  Is  not 
surprising,  since  the  molecule  is  nearly  symmetric  and  the  two  charge  states 
differ  modestly. 

Comparison  of  these  results  for  the  C**  and  O’*  Ions  with  extant  data  on 
ion-molecule  collisions  is  informative.  The  collisional  process  that  was  studied2^ 
Involved 

Ar13*  X  — —  Ar13*  +  fragments  (C'\  Of*,  e*).  0) 
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Energy  Spectrum  of  0*‘  from  CO,  Energy  Spectrum  of  O**  from  CO 

lit-  J.1U 


kinetic  energy  (eV)  kinetic  energy  (eV) 

Fig.  (3):  Kinetic  energy  spectra  for  O’*:  (a)  Ions  produced  from  CO,  show  a 

maximum  energy  -  70  eV;  (b)  Ions  produced  from  CO  exhibit  a  maximum 
energy  of  -  30  eV. 


with  X  denoting  CO  and  CO,.  The  collision  energy  corresponded  to  ~  1  MeV/au 
so  that  the  time  of  interaction  is  sufficiently  small  that  negligible  nuclear  motion 
will  occur  in  the  molecular  frame.  The  resulting  coulomb  explosion  commences 
from  an  intact  ionized  molecule  in  a  spatial  configuration  determined  by  the  neutral 
target  molecule. 

The  collisional  data^9  of  Matsuo  et  al.  for  carbon  and  oxygen  Ions  produced 

from  CO  and  CO,  are  shown  in  Fig.  (4).  Also  indicated  for  comparison  with  their 

data  are  the  maximum  kinetic  energies  of  the  O'**  and  O1*  Ions  observed  In  the 

multiphoton  studies.  In  all  cases,  the  multiphoton  results  correspond  rather 

closely  to  the  values  produced  by  the  Ion-molecule  interaction.  This  similarity  In 

outcome  leads  to  the  conclusion  that  the  nature  of  the  moleculai  ion  which 

undergoes  the  coulomb  explosion  Is  essentially  the  same  In  both  cases.  Specifically. 

we  conclude  that  the  strong-field  interaction  has  many  features  In  common  with 

ion-atom  and  ion-molecule  collisions.  The  nature  of  this  analogy  has  been  explored 

30  31 

previously  In  other  work. 
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IV.  X-RAY  GENERATION  FROM  SOLID  TARGETS 


Solid  materials  offer  the  possibility  of  combining  a  high  particle  density  with 

a  very  high  level  of  electronic  excitation,  if  the  deposition  of  energy  in  the  solid 

can  be  made  to  occur  sufficiently  rapidly.  High  brightness  subpicosecond  light 
2- 10 

sources  now  enable  such  rapid  excitation  to  be  achieved. 

Consider  a  slab  of  unit  area  with  density  p  and  thickness  5  which,  upon 
excitation  with  intense  subpicosecond  radiation,  produces  radiation  in  the  kilovolt 
spectral  range.  In  order  to  set  the  scale,  we  want  to  obtain  a  simple  upper 
bound  on  the  x-ray  intensity.  Since  the  time  of  excitation  is  short,  as  a  first 
approximation,  we  will  ignore  the  hydrodynamical  expansion,  a  process  which  will 
occur  roughly  at  a  rate  of  -  1-10  A/fs.  An  upper  limit  on  the  x-ray  intensity 
Ixm  can  then  be  simply  written  as 


(2) 


in  which  ho>x  represents  a  kilovolt  quantum  and  tx  is  the  effective  radiative  rate 

of  the  excited  material.  For  the  present  discussion,  we  will  set  p  as  solid  density 

(-  5  x  10”  cm'1)  and  take  as  the  lower  limit  for  tx  a  typical  radiative  rate  for 

kilovolt  transitions.  For  the  latter  will  use  the  range  2  x  I0*lw  ^  xx  <  10'13  s 

as  a  rough  measure  If  we  take  for  6  the  recently  measured  value33  of  the  depth 

of  energy  penetration  of  -  250  nm  for  MgF,/Si02  targets  irradiated  with  248  nm 

at  an  intensity  of  -  3  x  101*  W/cm3,  Eq.  (2)  can  be  evaluated.  With  these  parameters 

we  find  10IS  <  Ixm  <  5.  x  1015  W/cm3.  This  simple  upper  limit  on  Ixm  can  now 

34 

be  compared  to  the  experimental  result. 

An  experimental  concern  in  studies  of  the  interaction  of  intense  radiation 
with  solid  matter  is  the  influence  of  any  low  intensity  prepulse  on  the  true  state 

of  the  irradiated  system  It  has  been  experimentally  found33  that  the  interference 
of  the  prepulse  can  be  entirely  eliminated  at  248  nm  by  the  use  of  a  target 
composed  of  ultraviolet  transmitting  materials.  For  example.  MgF,  and  BaF. 
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exhibit  essentially  no  absorption  up  to  relatively  high  intensities  and  studies  have 

demonstrated  that  these  values  are  sufficiently  high  to  prevent  any  significant 

influence  oy  the  prepulse.  This  enables  the  intense  subpicosecond  energy  to 

33  34 

couple  directly  to  virgin  solid  matter. 

The  KrF*  system  used  in  this  experiment  delivered  a  subpicosecond  pulse 

energy  of  -  150  mj.  The  prepulse  had  a  duration  of  -  20  ns,  an  energy  of  ~  10 

mj,  and  is  known  to  be  very  poorly  focusable.8  The  radiation  was  focused  onto  a 

flat  rotating  8aF3  target  with  a  f/10  CaFa  lens.  The  peak  laser  intensity  achievable 

on  target  was  estimated  to  be  -  3  x  101*  W/cm2,  based  on  previous  threshold 

ionization  intensity  measurements  at  low  density35  as  well  as  recent  direct  imaging 

of  the  focal  zone.  The  x-ray  emission  from  the  target  was  collected  by  flat  KAP 

(2d  =  26. 6^)  and  PET  (2d  =  8.7&)  crystal  spectrometers  and  recorded  on  DEF  film. 

For  nigher  resolution,  a  convex  KAP  crystal  spectrometer  was  used.  Approximately 

10  shots  were  accumulated  for  each  spectrum. 

The  spectrum  observed  from  BaF,  exhibited  strong  emission  in  the  £-13  A 

range.  In  addition  to  the  spectrum,  calibration  of  the  film  revealed  that  -  6  mj 

of  x-rays  were  produced  on  each  shot.  Subsequent  pinhole  camera  measurements 

of  the  x-ray  emission  from  solid  targets  (MgFa/SiOa)  indicated  a  radiating  spot 

size  with  a  diameter  in  the  range  of  8-9  pm.  This  value  is  quite  close  to  the 

34 

independently  measured  focal  diameter  of  the  optical  system.  Finally,  calculations 
estimate  that  the  radiation  occurs  over  a  time  of  -  1  ps,  a  value  that  is  consistent 
with  the  x-ray  pinhole  camera  measurements.  Therefore,  assuming  isotropic 
emission  from  a  volume  of  these  dimensions,  the  experimentally  radiated  intensity 
Ir  was  approximately  Ir  -  2-3  x  10*’  W/cm2.  We  see  immediately  that  this  value 
falls  in  the  range  given  by  the  estimate  the  upper  bound  Ixm.  The  measured 
performance,  therefore,  appears  close  to  the  maximum  possible.  Although  such  a 
source  is  incoherent,  a  radiation  rate  of  this  magnitude,  coupied  with  the  relatively 
narrow  spectrum  over  which  it  occurs,  constitutes  a  source  of  extremely  high 
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1  23456  1  23456 

CHARGE  STATE  i 

Fig.  (4):  The  data  of  Ref.  (29)  on  carbon  and  oxygen  Ions  produced  In  Ar13* 

ion-molecule  collisions  In  the  reaction  Ar13*  +  X  — — *  Ar13*  fragments 
(C'*,  oi*.  e‘)  at  a  collision  energy  of  -  1  MeV/au.  The  multiphoton 
data  for  C"*  and  O1*  are  also  shown  for  comparison  for  CO,  (A)  and 
CO  (o).  Rather  close  agreement  of  the  collislonal  and  multiphoton  data 
is  seen.  Figure  reproduced  with  permission. 

spectral  brightness  that  can  produce  unusually  high  excitation  rates  in  absorbing 
material. 

IV.  SUMMARY  AND  CONCLUSIONS 

The  use  of  recently  developed  high  brightness  subpicosecond  light  sources 
for  the  study  of  the  Interaction  of  matter  with  strong  electromagnetic  fields  Is  leading 
to  the  observation  of  new  phenomena  and  new  classes  of  excited  states  of  matter. 
Two  important  consequences  of  these  developments  are  the  production  of  new  forms 
of  Ionized  molecular  matter  and  the  generation  of  powerful  x-ray  sources  associated 
with  high  density  plasma  environments. 
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Appendix  E: 


‘Stable  Channeled  Propagation  of  Intense  Radiation  in  Plasma  Arising 
from  Relativistic  and  Charge-Displacement  Mechanisms* 
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ABSTRACT 

Dynamical  calculations  of  propagating  intense  axisymmetric  circularly  polarized 
laser  pulses  in  coilisionless  plasmas,  incorporating  relativistic  and  charge-displacement 
mechanisms,  show  that  their  combined  action  can  lead  to  stabilized  confined 
modes  of  propagation.  The  effect  of  the  charge-displacement  is  large  and  cavitation 
of  the  electron  distribution  is  a  general  feature.  An  example  shows  that  ~  50%  of 
the  incident  power  (~  4  TW)  can  be  trapped  in  a  channel  which,  for  a  wavelength 
of  248  nm,  an  electron  density  of  -  75  x  lOao  cm'1,  and  an  initial  beam  radius 
of  ~  3  pm,  develops  propagating  intensities  of  -  10SI  W/cma. 
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The  dynamics  of  propagation  of  intense  subpicosecond  laser  pulses  in  plasmas 
in  the  strongly  relativistic  regime  is  currently  undergoing  vigorous  theoretical 
analysis1-10  Of  primary  significance  is  the  possible  formation  of  quasistable  self- 
channeled  modes  of  propagation.  In  this  Letter  we  examine  theoretically  the 
consequences  of  the  combined  action  of  relativistic  and  charge-displacement  mechanisms 
with  particular  attention  to  the  development  and  character  of  highly  stable  confined 
modes  of  propagation.  It  is  assumed  that  the  pulse  duration  t  considered  in  this 
treatment  satisfies  the  inequality  Tj  >  >  r  >  >  Te,  with  tj  and  te  designating  the 
response  times  of  the  ions  and  electrons,  respectively,  so  that  ion  displacement 
can  be  neglected.  For  radiation  with  a  wavelength  of  248  nm,  the  range  of  intensity 
of  principal  interest  is  I  >  1019  W/cm3,  since  both  relativistic  and  charge-displacement 
mechanisms  are  very  significant  in  that  region.  Additional  experimental  studies 
concerning  this  matter  are  presented  in  Ref.  (11). 

The  basic  formalism  and  procedures  describing  calculations  have  been  described 
previously.6  Assuming  collisionless  and  lossless  propagation,  the  slowly  varying 
complex  field  amplitude  E(t,z,r)  is  governed  by  the  nonlinear  Schrodinger  equation 
(3.  6-8] 

/  1  3  3  \  i 

( - -  —  ) E  -t-  -  [4E  t  k3(1  -  rlNe)E]  =  0,  (1) 

\Cl  3t  3z )  2k  x 

Ne  =  max(0,f(r)  ♦  kp3AxY).  Y  «  (1  *  I/Ir)l/a  (2) 

Here  =  33/3r3  +  r'l3/3r  (r  is  the  transverse  coordinate),  Ci  is  the  group 

velocity  in  the  plasma,  k  is  the  wave  number,  k3  =  kg  -  kp,  k0  =  u/c,  kp  =  Up_0/c, 
to  is  the  laser  angular  frequency,  Up  0  is  the  plasma  frequency,  c  is  the  vacuum 
speed  of  light.  Ne  denotes  the  electron  density  N  normalized  by  the  peak  value 
of  the  initial  (unperturbed)  electron  density  distribution  N0;  f(r)  denotes  the 
initial  transverse  distribution  of  N0.  I  =  (c/4ir)  |  E  | 3  is  the  intensity  of  laser 
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beam,  Ir  =  me20  wV/4ire2  is  the  parameter  called  the  relativistic  intensity.  The 

term  kp2  AXT  in  (2)  describes  the  charge-displacement  process. 

The  ponderomotive  force  is  capable  of  completely  expelling  electrons  from  certain 
spatial  regions,  a  condition  described  as  "electronic  cavitation".3  The  calculations 

show  that  the  charge-displacement  plays  a  strong  role  in  determining  the  spatial 
character  of  the  propagating  energy  and,  particularly,  in  stabilizing  the  confined 
high-intensity  modes. 

The  essential  finding  of  these  calculations  is  insight  into  the  formation  and 
stabilization  of  spatially  confined  modes  of  propagation  arising  from  the  cooperative 
action  of  the  relativistic  and  charge-displacement  processes.  The  numerical  simulations 
were  performed  with  model  pulses  having  Gaussian  or  hyper-Gaussian  incident 
transverse  intensity  distributions  interacting  with  initially  homogeneous  plasmas  or 
plasma  columns  having  a  hyper-Gaussian  radial  profiles. 

The  results  of  computations  presented  in  this  work  are  given  in  the  frame  (q.z.r ), 
q  =  t  -(z/cj,  connected  with  the  wave  front  of  the  beam.  The  graphical  data 

represented  below  illustrate  the  propagation  of  the  radiation  along  the  z-axis  for 
q  =  constant.  We  present  the  results  of  simulations  of  dynamics  of  initially 
Gaussian  beams  for  the  following  values  of  the  parameters:  X  =  248  nm,  Ir  =  1.34 
x  102°  W/cm2,  peak  value  of  incident  transverse  intensity  distribution  I0  =  3  x 

101S  W/cm2,  initial  beam  radius  r0  =  3  pm,  N0  =  7.5  x  1020  cm'1,  and  f(r)  a  i 
(initially  uniform  plasma).  In  this  case,  the  ratio  of  the  beam  power  P0  to  the 
critical  power  (Pcr)  of  relativistic  and  charge-displacement  self-focusing3 


cr 


(m2  ac5/e2)(  f  g2(p)pdp  j(u/u)p  0) 


is  22  252.  In  the  above  expression,  g0(p)  is  the  so-called  Townes  mode.12 
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The  solutions  for  the  spatial  profiles  of  the  laser  intensity  [Fig.  (i)]  and  the 
corresponding  electron  density  [Fig.  (2)]  exhibit  several  important  characteristics. 
Among  them  are  (1)  the  generation  of  intense  focal  regions,  (2)  the  development 
of  a  stable  confined  channel  of  high-intensity  propagation,  (3)  the  trapping  of  a 
substantial  fraction  of  the  incident  power  in  the  channel,  and  (4)  strong  cavitation 
of  the  electron  density. 

On  the  basis  of  tne  dynamical  picture,  the  combined  effect  of  the  relativistic 
nonlinearity  and  the  charge-displacement  can  be,  in  reasonable  approximation, 
summarized  in  the  following  simple  way.  The  relativistic  effect  leads  to  the 
initial  concentration  of  the  radiation  and  the  resulting  displacement  of  electronic 
charge  reinforces  this  tendency  and  stabilizes  the  confinement.  The  cooperative 
nature  of  this  action  appears  to  lead  to  highly  stable  conditions  of  propagation. 

The  behavior  of  the  computed  solutions  for  large  z  can  be  compared  with 
the  z-independent  eigenmodes8  of  the  nonlinear  Schrodinger  equation  (1).  The 
main  result  of  this  work  is  that  the  field  amplitude  tends  asymptotically  to  the 
lowest  eigenmode  Vs(r)  of  this  equation.  We  call  this  phenomenon  relativistic 
and  charge-displacement  self-channeling.  Vs(r)  is  a  real-valued  positive  monotonic 
function  of  r,  carries  a  real  dimensionless  index  s,  vanishes  as  r  — ■  ®,  and  is 
related  to  the  electric  field  amplitude  by  the  statement3'8 

E(r.z)  =  Vs(r)  exD[iz(s  -  1)kp(2k)'1].  (3) 

Computations  have  shown,  for  the  case  studied,  that  the  field  distribution 
for  large  z  tends  to  the  amplitude  Vs(r)  for  the  index  s  =  0.554.  The  appropriate 
intensity  distribution  is  then  given  by  Is(r)  a  Vs3(r).  Figure  3  exhibits  the  normalized 
asymptotic  field  amplitude  (Is(r)/I0)1/3  and  the  corresponding  asymptotic  electron 
density  distribution  Ns(r)/N0.  In  the  asymptotic  distribution  of  intensity,  46%  of 
the  initial  power  of  the  beam  is  retained. 
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Several  other  cases  have  also  been  considered.  In  particular,  calculations 
examining  the  dynamics  of  propagation  have  been  studied  for  hyper-Gaussian  initial 
intensity  distributions  for  both  homogeneous  plasmas  [f  (r)  *  i]  and  simulated 
plasma  columns  [f(r)  =  exp  (-{r/r0)*j.  The  behavior  was  found  to  be  qualitatively 
the  same  as  that  shown  in  Fig.  (1)  -  Fig.  (3)  with  asymptotically  trapped  fractions 
of  the  initial  laser  power  ranging  from  34%  to  77%.  Generally,  studies  of  the 
relativistic/charge-displacement  self-channeling  thresholds  showed  that  the  only 
condition  required  for  initially  flat  waveforms  to  propagate  in  this  regime  is  P0  > 
Pcr.  Self-channeled  propagation  also  develops  in  the  cases  involving  initially 
focused  and  defocused  beams.  We  note  that  the  relativistic  and  charge-displacement 
self-channeling,  characterized  by  the  asymptotic  evolution  to  the  distributions 
given  by  the  lowest  stationary  eigenmodes,  described  above,  differs  essentially 
from  the  previously  known  oscillatory  regime  arising  from  relativistic  and  charge- 
displacement  mechanisms3  and  the  stabilization  seen  in  purely  relativistic  propagation 
in  plasma  columns.6 

In  summary,  the  dynamical  behavior  of  channeled  propagation  arising  from 
both  relativistic  and  charge-displacement  mechanisms  has  been  studied  theoretically 
Two  principal  findings  have  been  established.  (1)  These  two  processes  can 
cooperatively  reinforce  one  another  and  lead  to  stable  confined  high-intensity 
modes  of  propagation  in  plasmas  which  are  capable  of  trapping  a  substantial 
fraction  (-  50%)  of  the  incident  power.  (2)  The  spatial  profiles  of  both  the 
intensities  and  electron  density  distributions  tend  asymptotically  to  those  derived 
as  lowest  eigenmodes  of  a  z-independent  analysis. 
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FIGURE  CAPTIONS 


Fig.  1:  The  normalized  intensity  distribution  I(r,z)/I0.  I0  =  3  x  101*  W/cm5  r0 

=  3  pm,  X  =  248  nm,  and  N0  =  7.5  x  IQ20  cm'3. 

Fig.  2:  The  normalized  electron  density  distribution  N(r,z)/N0  corresponding  to 

the  data  shown  in  Fig.  (1).  N„  =  7.5  x  IQ20  cm'3. 

Fig.  3:  Radial  dependence  of  the  asymptotic  solutions  for  the  normalized  amplitude 

[Ig(r)/Ia]l/2  and  the  normalized  electron  density  Ns(r)/N0  for  the  index 
s  =  0.554. 
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Figure  1 
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Figure  2 
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Appendix  F: 


"Observation  of  Relativistic  and  Charge-Displacement  Self-Channeling 
of  Intense  Subpicosecond  Ultraviolet  (248  nm)  Radiation  in  Plasmas" 
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Experimental  studies  examining  a  new  relativistic  regime  of  high-intensity  short-pulse  propagation  in 
plasmas  have  been  performed  which  present  evidence  for  the  formation  of  a  stable  mode  of  spatially 
confined  (channeled)  propagation.  For  an  electron  density  of  —  l.35x|0:i  cm"'  and  a  power  of 
—  3x  10"  W.  the  results  indicate  a  channel  radius  <  I  pm  and  a  peak  intensity  —10”  W/cm:.  Com¬ 
parison  of  these  findings  with  a  dynamical  theory  yields  agreement  for  both  the  longitudinal  structure 
and  the  radial  extent  of  the  propagation  observed. 


P ACS  numbers  52.-IO.Db.  42.2J.-p.  42.65.Jx 

A  fundamentally  new  regime  of  electromagnetic  propa¬ 
gation  is  expected  to  arise  in  plasmas  for  short-pulse  radi¬ 
ation  at  sufficiently  high  intensity.  Calculations  of  the 
propagation  in  plasmas,  incorporating  both  relativistic 
[1.2]  and  charge-displacement  mechanisms  [3—7],  indi¬ 
cate  that  the  combined  action  of  these  processes  can  lead 
to  a  new  stable  form  of  spatially  confined  (channeled) 
propagation.  This  Letter  (I)  reports  the  results  of  the 
first  experimental  study  probing  the  physical  regime 
relevant  to  the  observation  of  relativistic  and  charge- 
displacement  self-channeling  and  (2)  presents  the  initial 
comparison  of  these  experimental  findings  with  matching 
theoretical  calculations  performed  with  the  computation¬ 
al  procedures  described  in  Ref.  [7], 

The  experimental  arrangement  used  in  these  studies  is 
illustrated  in  Fig.  1(a).  The  source  of  radiation  was  a 
subpicosecond  KrF*  (A.  “248  nm)  laser  that  has  been  de¬ 
scribed  elsewhere  [8],  It  delivered  a  linearly  polarized 
power  of  —  3x|0"  W  (  —  150  mJ,  pulse  duration  —500 
fs)  in  a  beam  with  a  diameter  of  —42  mm.  When  this 
radiation  was  focused  into  the  chamber  with  lens  L  l 
(// 7).  a  focal  radius  r0—3. 5  pm  was  measured,  giving  a 
maximum  intensity  /o— 8.6x|017  W/cm2.  The  medium 
was  provided  by  filling  the  chamber  statically  with  gas 
[He,  Ne.  Ar,  Kr.  Xe,  N;,  CO;,  or  a  mixture  of  Xe  (4%) 
and  N;  (96c7)l  up  to  a  maximum  density  of  —  l.89x  |0:° 
cm  ~  \ 

The  diffracted  248-nm  radiation  was  measured  as  a 
function  of  the  angle  (0)  with  respect  to  the  direction  of 
the  incident  radiation.  The  incident  laser  beam  was 
blocked  by  a  metal  disk  on  the  output  window  of  the 
chamber  and  lens  L  2  imaged  the  region  near  the  focal 
zone  on  a  fluorescent  screen  5.  The  diaphragm  D  in  front 
of  lens  L  2  restricted  the  collection  of  the  diffracted  light 
to  a  solid  angle  of  —5°  while  simultaneously  increasing 
the  depth  of  field.  The  angle  between  the  axis  of  the  lens 
LI  and  the  axis  of  the  incoming  laser  radiation  could  be 
readily  varied  up  to  a  maximum  angle  of  0—15°  Two 
flat  mirrors  coated  for  high  reflection  (—  99%)  at  248 
nm.  both  having  a  spectral  bandwidth  of  —  10  nm.  served 


in  reflection  as  spectral  filters  (F)  for  the  diffracted  laser 
radiation  so  that  only  the  scattered  248-nm  radiation 
could  illuminate  the  screen.  An  attenuator  A  was  em¬ 
ployed  to  adjust  the  intensity  on  the  screen  and  the  im¬ 
ages  formed  were  recorded  through  the  visible  fluores- 


Axial  Distance  z  (nm) 

FIG  I  (a)  Experimental  apparatus  used  in  studies  of  propa¬ 
gation.  See  text  for  description,  (b)  Data  concerning  the  pat¬ 
tern  of  propagation  observed  with  a  single  pulse  in  N;  ai  j  den¬ 
sity  of  —  1 .35  x  I020  cm  The  maximum  intensity  is  half  the 
detector  (CCD)  saturation.  The  radiation  is  incident  from  ihe 
left.  Inset:  Photographic  data  with  a  vertical  spatial  resolution 
of  —10  pm.  The  graph  illustrates  the  one-dimensional  axial 
profile  taken  along  the  direction  of  propagation  (:)  of  ihe  pho¬ 
tographic  data  (inset).  The  spacing  of  the  maxima.  0  =  300 
±  20  uni.  is  indicated. 
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cence  produced  with  a  microscope  and  a  charge-coupled- 
device  (CCD)  camera. 

The  characteristic  behavior  observed  is  well  illustrated 
by  the  data  recorded  with  Nj.  The  measured  result, 
shown  in  Fig.  1(b),  corresponds  to  a  density  pn,~I  35 
*  I030  cm  -3.  To  the  left,  in  the  photographic  inset,  a  rel¬ 
atively  large  cone  of  light  Rayleigh  scattered  from  the 
plasma  is  visible  at  all  angles  as  the  energy  propagates  to¬ 
ward  the  focal  point  of  the  lens,  while,  in  the  region  to  the 
right  of  the  conical  apex,  a  narrow  filament  developed. 
The  diameter  of  this  filament  is  not  greater  than  10  pm. 
the  measured  spatial  resolution  of  the  imaging  system. 
The  distribution  of  intensity  observed  along  the  filament 
exhibited  several  bright  features  attributed  to  diffraction 
because  they  could  not  be  seen  for  8>  20°.  Since  the 
axis  of  the  imaging  lens  corresponded  to  an  angle 
8  ■■7. 5°,  the  scale  along  the  abscissa  of  the  photographic 
data  is  reduced  by  almost  eightfold,  giving  the  maximum 
length  of  the  filament  as  —  I  mm.  The  graph  in  Fig.  I  (b) 
represents  a  one-dimensional  axial  profile,  taken  along 
the  direction  of  propagation  (r).  of  the  observed  intensity 
pattern  (inset).  Three  peaks  (a,f),  y)  are  visible  with  a  a 
spatial  separation  of  $"200  ±  20  pm.  The  normal  Ray¬ 
leigh  range  for  the  focal  geometry  used  was  —200  pm. 

The  diameter  of  the  filamentary  channel  is  an  impor¬ 
tant  dynamical  variable,  which  we  estimated  by  measur¬ 
ing  the  maximum  angular  deviation  of  the  diffracted 
light.  The  experimental  value  <t>  of  this  diffracted  cone 
was  «— 20°,  a  magnitude  indicating  a  radius  /•,— 0.9  pm 
though  the  relation  $™1.22A./ry  Filaments  of  this  gen¬ 
eral  nature  were  observed  at  densities  above  — 1.35  x  I0:o 
cm  "  '  in  N;,  Ne,  Ar.  Kr.  COi.  and  a  mixture  of  Xc  (4%) 
and  N:  (96%).  but  not  in  He  and  Xc.  two  materials  dis¬ 
cussed  further  below 

Two  mechanisms  exist  that  could  modify  the  refractive 
index  of  the  medium  and  lead  to  the  observed  behavior. 
They  are  (I )  the  Kerr  effect  stemming  from  the  ions  and 
(2)  the  relativistic  and  charge-displacement  process  [7], 
Since  the  pulse  duration  is  very  short  (  —  500  fs),  the 
motion  of  the  ions  is  negligible  (3).  and  no  contribution 
can  arise  from  expulsion  of  the  plasma  from  the  high- 
intensity  zone.  An  implication  of  the  estimate  of  the 
channel  radius  (r#— 0.9  pm)  is  that  the  observed  propa¬ 
gation  is  associated  with  intensities  in  the  1 0 1>‘—  1 0 
W/cm:  range.  Under  these  conditions,  available  experi¬ 
mental  evidence  (9.10)  on  multiphoton  ionization  indi¬ 
cates  that  He  should  be  fully  ionized  and  the  C.  N.  and  O 
atoms  constituting  the  molecular  materials  would  retain, 
at  most,  only  1 5  electrons. 

Consider  explicitly  the  case  of  N;,  which  has  estimated 
(9.10)  threshold  intensities  for  the  production  of  N(*. 
N**.  and  N7  +  of  l.6x|0'‘,  6.4x|0'\  and  l.3x|0,'> 
W/cm:.  respectively.  Hence,  the  volume  of  the  channel 
would  be  largely  ionized  to  N,+.  with  certain  localized 
high-intensity  regions  contributing  some  N6+.  Two 
consequences  of  this  pattern  of  ionization  follow,  namely, 
(i)  the  Kerr  effect  arising  from  the  ions  is  small,  since  the 


polarizabilities  of  the  remaining  1 5  electrons  are  low,  and 

(11)  the  electron  density  in,)  initially  produced  in  the  fo¬ 
cal  region  is  nearly  uniform.  Therefore,  nr  =*  1 .35  x  |0:i 
cm  “ 3  for  the  data  on  Ns  shown  in  Fig.  I  (b). 

A  critical  power  Pc,  for  self-channeling,  arising  from 
the  relativistic  and  charge-displacement  mechanism,  can 
be  defined  [II]  as 

/\r  -  (m-c -/e  ’ )  gs  (r)r  dr(co/<or  0) : 

=  !.62x|Olo(ncr/n,)  W.  (|) 

where  m,.,  e.  and  c  have  their  customary  identifications,  cu 
is  the  laser  angular  frequency.  (or  o  is  the  plasma  frequen¬ 
cy  for  the  uniform  unperturbed  plasma  with  electron  den¬ 
sity  ne.  ncr  is  the  critical  electron  density  (for  X  ”248  nm. 
nLr“  l.82x  10"  cm-3),  and  go(r)  is  the  Townes  mode 

(12) . 

The  critical  powers  associated  with  the  experimental 
conditions,  for  He  and  N;  at  a  medium  density 
p-l.35x  |020  cm-3,  are  l.08x|0,;!  and  2.l9x|o"  W, 
respectively.  Therefore,  since  the  incident  power  was 
P  =  3x  10"  W.  no  filament  was  expected  in  He.  a  predic¬ 
tion  conforming  with  the  observation  of  none.  Moreover, 
the  diffracted  cone  of  radiation  was  also  absent  with  He 
In  contrast.  P/Pcr—\.37  for  N2.  a  condition  that  held 
generally  (PI Pa  >  I )  for  all  materials  which  exhibited 
evidence  for  channel  formation.  We  note,  however,  that 
some  contribution  from  the  Kerr  effect  may  be  present, 
even  for  the  light  materials  (Ne.  N;.  and  CO;),  in  the 
early  stage  of  channel  formation  prior  to  the  development 
of  a  substantial  level  of  ionization,  and  that  the  heavier 
gases  (Ar.  Kr,  and  Xe)  may  involve  a  more  significant 
influence  from  the  Kerr  process.  A  specific  estimate  of 
the  nonlinear  index  change  arising  from  both  N>+  and 
N6+  at  an  intensity  of  — 10'9  W/cm:  indicates  that  their 
contribution  is  less  than  10  -  3  that  of  the  free  electrons, 
hence  the  ionic  contribution  can  be  neglected  in  N:  for 
the  conditions  studied. 

A  direct  comparison  will  now  be  made  between  the 
theoretical  analysis,  fully  described  in  Ref.  (7).  and  the 
experimental  findings  for  Ns.  This  comparison  can  be  ac¬ 
complished  for  both  the  longitudinal  intensity  profile  and 
the  radial  extent  of  the  channel.  Figure  2(a)  illustrates 
the  intensity  profile  l(r,:)/l 0  calculated  with  physical  pa¬ 
rameters  corresponding  to  those  of  the  experiment  ( i.e  . 
F  — 3x|0"  W,  r0“3.5  ^m.  n,  *l.35x  I0:|  cm-',  and 
PI  Pi :r“l  37).  Importantly,  all  of  these  parameters  are 
based  on  independent  measurements  of  (I )  the  laser  pulse 
involving  determinations  of  the  energy  and  power  1  PI 
(2)  the  focal  radius  (r0)  of  the  incident  radiation,  and  <  3 1 
the  characteristics  of  the  multiphoton  ionization  (10)  gen¬ 
erating  the  electron  density  (n, ).  Therefore,  this  compar¬ 
ison  does  not  involve  a  fit  with  a  free  parameter  The 
normalized  electron  density  calculated  is  presented  in 
Fig  2(b).  from  which  it  is  seen  that  electronic  cavitation 
occurs  only  near  the  positions  of  the  maxima  in  the  inten¬ 
sity  profile  [Fig.  2(a)).  Curve  A  in  Fig.  2(c)  represent' 
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Radial  Distance  r(um) 


the  one-dimensional  axial  intensity  profile  HO.:)/ to  cor¬ 
responding  to  the  calculated  distribution  shown  in  Fig. 
2(a).  The  spacing  (6)  of  the  maxima  is  seen  to  be 
S—  185  pm,  a  value  in  close  agreement  with  the  experi¬ 
mental  figure  (6  — 200±20  am)  illustrated  in  Fig.  Mb) 
Furthermore,  analysis  has  shown  that  the  spacing  o  is 
quite  sensitive  to  the  power  P o  and  electron  density  n,.. 
particularly  in  the  region  close  to  the  threshold  (see  Eq. 
(1)1.  With  respect  to  the  results  illustrated  in  Fig.  2(c). 
an  increase  in  nr  by  less  than  5%  causes  a  reduction  in 
the  spacing  5  by  approximately  25%.  Therefore,  substan¬ 
tial  changes  in  the  physical  parameters  would  grossly 
alter  the  comparison  of  the  experimental  and  theoretical 
results. 

Theoretical  studies  [3-7,11]  indicate  that  the  charge 
displacement  plays  a  very  important  dynamical  role.  In 
order  to  test  this  hypothesis,  identical  calculations  were 
made  for  N>  for  the  purely  relativistic  case  [2]  which  ex¬ 
plicitly  neglects  the  charge-displacement  term,  namely, 
elimination  of  the  term  (o2/tu/io/’ij)Aj.(l +///,)l/:  in  Eq. 
(24)  of  Ref.  [7],  Significantly,  the  resulting  axial  profile 
[curve  B  in  Fig.  2(c)]  exhibits  only  a  single  relatively 
weak  maximum,  for  0  <  z  <  600  am.  an  outcome  sharply 
at  variance  with  both  the  full  theoretical  analysis  and  the 
experimental  observation.  Although  the  expected  charge 
displacement  is  highly  localized  [Fig.  2(b)] .  this  compar¬ 
ison  reveals  the  strong  influence  it  has  on  the  propagation 
[3.7,11],  At  a  greater  incident  power  (P/P«— 10).  a 
continuous  channel  in  the  electron  distribution  is  expect¬ 
ed  to  develop  [7], 

The  measurements  indicated  an  approximate  value  of 
r#— 0.9  am  for  the  radial  extent  of  the  channel,  a  result 
that  can  be  compared  with  the  corresponding  theoretical 
figure.  Figure  2(d)  illustrates  five  radial  intensity  profiles 
l(r.:,)l to  of  the  distribution  pictured  in  Fig.  2(a).  Since 
the  measurement  of  this  angularly  scattered  radiation  did 
not  correspond  to  a  known  longitudinal  position,  this 
comparison  can  only  be  qualitative,  but  the  radial  distri¬ 
butions  shown  indicate  that  the  expected  value  lies  in  the 
interval  0.5  1.0  am.  a  range  that  comfortably  in¬ 

cludes  the  experimental  value  r* 

The  results  observed  with  Xe  deserve  additional  discus¬ 
sion.  since  those  experiments  did  not  give  evidence  for  the 
formation  of  a  channel.  In  significant  contrast  to  the  case 
involving  N;,  the  electron  density  n,  produced  by  the 
multiphoton  ionization  [10]  in  Xe  is  expected  to  be  ver> 


FIG.  2.  Calculations  for  N;  with  /»- 3x  10"  W.  m -?  5  am. 
«..-l  35*  I0:l  cm*',  and  /0-8.6xlO'7  W/cm:  <a)  Normal¬ 
ized  intensity  /(/-..-  >//o.  (b)  Normalized  electron  densiiy 

V(r,c )//Vo  for  N;  with  <V0-nr.  (c)  Normalized  one-dimen¬ 
sional  axial  intensity  profiles  /(0.;)//o.  Curve  A.  full  theory  lor 
data  in  panel  (a).  5-185  am.  Curve  B.  calculation  with 
charge-displacement  term  neglected,  (d)  Normalized  radial  in¬ 
tensity  profiles  Hr.:,)/lo  corresponding  to  panel  (a)  Longitu¬ 
dinal  positions  :<  — 172  a m.  ::  —  245  am.  c i  —358  am.  -  j  “  44 1 
am.  and  —  559  am  and  r#— 0.9  am. 
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nonuniform  spatially.  For  intensities  spanning  10I6-I0U 
W/cm:,  the  corresponding  density  n,  would  vary  by  over 
a  factor  of  2.  Since  this  nonuniformity  would  tend  to 
reduce  the  refractive  index  locally  in  the  central  high- 
intensity  region,  a  significant  defocusing  action  is  expect¬ 
ed  which  could  suppress  the  channel  formation. 

Finally,  we  note  (I)  that  the  intensity  distribution  is 
not  expected  to  depend  strongly  on  the  state  of  polariza¬ 
tion  (13.14]  and  (2)  that  losses  to  the  plasma  may  be 
significant,  particularly  at  electron  densities  close  to  nu. 

The  first  experiments  examining  a  new  relativistic  re¬ 
gime  of  high-intensity  pulse  propagation  in  plasmas  have 
been  performed  and  the  findings  indicate  the  formation  of 
a  channeled  mode  of  propagation  over  a  length  consider¬ 
ably  greater  than  the  Rayleigh  range.  Specific  compar¬ 
isons  of  the  experimental  observations  with  a  dynamical 
theory,  which  explicitly  includes  both  the  influence  of  the 
relativistic  mass  shift  and  the  displacement  of  the  elec¬ 
tronic  component  of  the  plasma,  produce  excellent  agree¬ 
ment  for  both  the  longitudinal  structure  of  the  intensity 
profile  and  the  radial  extent  of  the  channel.  W'hile  the 
present  channel  contains  several  foci,  a  continuous  chan¬ 
nel  is  predicted  to  develop  at  higher  power.  Finally,  the 
intrinsically  very  high  concentration  of  power  associated 
with  this  mechanism  of  channeled  propagation  provides 
an  efficient  and  general  method  for  the  production  of  con¬ 
ditions  necessary  for  x-ray  amplification  (15]. 
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A  simple  derivation  in  the  Coulomb  gauge  of  the  nonlinear  Schrodinger  equation  describing  propaga¬ 
tion  of  powerful  ultrashort  circularly  polarized  laser  pulses  in  underdense  cold  inhomogeneous  plasmas 
is  presented.  Numerical  solutions  are  given  for  the  two-dimensional  axisymmetric  case  for  both  initially 
homogeneous  plasmas  and  static  preformed  plasma  columns.  These  solutions  account  for  (i)  diffraction, 
iii)  refraction  arising  from  variations  in  the  refractive  index  due  to  the  spatial  profile  of  the  electron  den¬ 
sity  distribution,  mil  the  relativistic  electronic  mass  shift,  and  (iv)  the  charge  displacement  resulting  from 
the  transverse  ponderomotive  force.  The  most  important  spatial  modes  of  propagation  corresponding  to 
1)  purely  relativistic  focusing  and  >2)  the  combined  action  of  both  the  relativistic  and  charge- 
displacement  mechanisms  are  described.  The  latter  leads  to  the  formation  of  stable  confined  modes  of 
propagation  having  paraxially  localized  regions  of  high  intensity  and  corresponding  paraxially  situated 
cavitated  channels  in  the  electron  density.  It  is  further  demonstrated  that  the  dynamical  solutions  of  the 
propagation  tend  asymptotically  to  the  lowest  eigenmodes  of  the  governing  nonlinear  Schrodinger  equa¬ 
tion.  Finally,  the  calculations  illustrate  the  dynamics  of  the  propagation  and  show  that  the  relativistic 
mechanism  promotes  the  initial  concentration  of  the  radiative  energy  and  that  the  subsequent  charge 
displacement  stabilizes  this  confinement  and  produces  waveguidelike  channels. 

PACS  number(s):  52.40.Nk,  42.65.Jx,  52.35.Mw,  52.40.Db 


I,  INTRODUCTION 

The  interaction  of  reiativistically  intense  subpi¬ 
cosecond  laser  pulses  with  gaseous  media  has  been  an 
area  of  vigorous  research  for  the  past  several  years.  For 
ultraviolet  wavelengths  on  the  order  of  200-300  nm.  the 
intensity  region  of  interest,  in  which  relativistic  effects  be¬ 
come  important,  lies  above  —  1018  W/cmJ.  The  propaga¬ 
tion  of  radiation  in  such  media,  for  intensities  greater 
than  -  1016  W/cm;,  naturally  causes  strong  nonlinear 
ionization  in  ail  matter.  Hence,  the  pulse  itself,  even  in 
regions  where  the  intensity  is  relatively  low  compared  to 
the  peak  value,  removes  many  electrons  [1,2]  from  the 
atomic  or  molecular  constituents  creating  a  plasma 
column  in  which  the  main  high-intensity  component  of 
the  pulse  propagates.  Therefore,  in  a  reasonable  first  ap¬ 
proximation,  the  investigation  of  the  resulting  propaga¬ 
tion  can  be  divided  into  two  separate  and  distinct  areas. 
They  are  tl)  the  atomic  and  plasma  physics  occurring  in 

45 


the  field  of  the  intense  electromagnetic  wave  leading  to 
the  ionization  and  (2)  the  subsequent  nonlinear  propaga¬ 
tion  of  the  radiation  in  the  plasma  that  is  generated.  The 
work  described  below  concerns  the  latter  issue. 

To  our  knowledge,  Akhiezer  and  Polovin  published  the 
first  treatment  of  high-intensity  electromagnetic  waves  in 
a  cold  plasma  [3],  They  derived  the  equations  describing 
the  propagation  as  a  function  of  the  single  canonical  ar¬ 
gument  (a >t—kz)  appropriate  for  one-dimensional  wave 
motion,  reduced  the  problem  to  Lagrangian  form  with 
two  integrals  of  motion,  and  presented  either  exact  or  ap 
proximate  solutions  corresponding  to  particular  vim 
plified  special  cases. 

Several  subsequent  treatments  have  been  devoted  * 
the  acceleration  of  charged  particles  in  either  relativist 
beat  waves  or  in  the  tail  of  a  single  relativistic  puNe 
Specifically,  Noble  [4]  applied  the  equations  derived  -• 
Akhiezer  and  Polovin  to  the  study  of  single-  and  douK 
wave  propagation.  In  other  work  [5]  beat  waves  ;n 
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plasmas  are  described  kineticaUy.  The  excitation  of  par¬ 
ticles  in  cold  plasmas  occurring  in  the  tail  of  a  single  in¬ 
tense  laser  pulse  has  also  been  investigated  [6]. 

Importantly,  this  previous  work  [6],  which  considered 
the  case  of  linear  polarization,  established  that  harmonics 
of  the  fundamental  frequency  are  not  produced  at  a 
significant  power,  if  the  plasma  is  underdense  and  the 
group  velocity  of  light  is  close  to  the  speed  of  light. 

Utilizing  the  quasistatic  approximation  for  the  fluid 
equations  of  cold  underdense  plasmas,  Sprangle,  Esarey, 
and  Ting  [7]  derived  the  first  nonlinear,  fully  self- 
consistent  set  of  equations  describing  the  propagation  of 
relativistic  laser  pulses  into  plasmas  for  a  one¬ 
dimensional  geometry.  This  formulation  was  used  to  ob¬ 
tain  insight  into  relativistic  self-focusing,  including  the 
self-consistent  electron-density  profiles,  wake-field  gen¬ 
eration,  optical  guiding,  and  second-harmonic  genera¬ 
tion.  However,  since  this  treatment  was  one-dimensional, 
the  multidimensional  case  remained  unsolved. 

An  important  finding  of  Sprangle.  Esarey,  and  Ting  [7] 
concerned  the  diffractional  erosion  of  the  leading  edge  of 
a  pulse  propagating  in  the  plasma.  However,  in  the  work 
described  below,  we  concentrate  on  the  evolution  of  the 
central  portion  of  the  pulse  as  it  propagates  in  the  plasma 
and  neglect  the  associated  erosion  of  the  leading  edge. 

Extensive  literature  exists  on  the  motion  of  electrons  in 
radiation  fields  of  certain  explicitly  given  forms.  The 
solution  for  the  case  of  plane  monochromatic  waves  is 
known  [8],  and  this  result  has  been  considerably  extended 
in  subsequent  work  [9,10].  Of  course,  the  well- 
established  Volkov  solution  for  the  Dirac  equation  also 
exists  [11].  In  particular.  Bardsley,  Penetrante,  and  Mit- 
tleman  [12]  have  numerically  simulated  the  relativistic 
dynamics  of  electrons  in  a  one-electron  picture  that  in¬ 
cludes  the  effects  of  space  charge  and  the  spatial  distribu¬ 
tion  of  the  radiation  field. 

The  first  treatment  of  the  relativistic  self-focusing  in 
plasmas  was  developed  by  Max,  Arons,  and  Langdon 
[13].  In  addition,  the  general  character  of  the  elec¬ 
tromagnetic  propagation  in  plasmas  has  undergone  con¬ 
siderable  analysis.  Schmidt  and  Horton  [14],  Hora  [15], 
and  Sprangle  et  al.  [16,17]  have  evaluated  the  thresholds 
for  relativistic  self-focusing  using  analytic  methods. 
Especially  germane  is  the  work  of  Sun  et  al.  [18],  who, 
for  initially  homogeneous  plasmas,  derived  the  two- 
dimensional  i  r,r  i  nonlinear  Schrodinger  equation  govern¬ 
ing  propagation,  including  consideration  of  the  combined 
effect  of  the  relativistic  nonlinearity  and  charge  displace¬ 
ment.  This  work  presented  the  lowest  eigenmode  of  the 
nonlinear  Schrodinger  equation  and  included  numerical 
evaluation  of  the  threshold  of  relativistic  self-focusing, 
the  value  of  which  has  been  approximately  estimated  in 
earlier  work  [14-17],  Finally,  this  analysis  [18]  presented 
the  lr,z)  dynamics  of  the  propagation  for  the  perturbed 
lowest  eigenmode  of  the  nonlinear  Schrodinger  equation 
for  cases  not  involving  spatial  cavitation  of  the  electron 
density.  We  note  that  Kurki-Suonio,  Morrison,  and  Taji- 
ma  [19]  have  also  developed  the  stationary  analytic  solu¬ 
tions  to  this  equation  for  a  one-dimensional  geometry. 

Additional  related  works  can  be  cited  in  this  context. 
Relativistic  self-focusing  and  beat-wave  phase-velocity 

i  ' 


control  in  plasma  accelerators  are  kindred  subjects  [20], 
Computations  involving  particle  simulations  revealing 
the  initial  process  of  self-focusing  and  subsequent  pon- 
deromotively  driven  electron  motions  have  been  per¬ 
formed  [21].  Other  calculations  analyzing  the  plasma  dy¬ 
namics  and  self-focusing  in  heat-wave  accelerators  [22], 
as  well  as  the  consideration  of  the  nonlinear  focusing  of 
coupled  waves  [23],  have  also  appeared. 

In  prior  publications,  we  have  (1)  investigated  the  gen¬ 
eral  behavior  of  two-dimensional  (r,z)  axisymmetric  rela¬ 
tivistic  self-focusing,  (2)  presented  results  on  the  process 
of  stabilization  of  laser  pulses  in  plasma  columns  [24].  '3i 
described,  with  the  use  of  an  analytical  model  [25],  the 
steady-state  characteristics  in  cavitated  channels  having 
overdense  walls,  (4)  reported  preliminary  results  of  calcu¬ 
lations  that  evaluated  the  combined  action  of  the  relativ¬ 
istic  and  charge-displacement  mechanisms  and  indicated 
the  formation  of  stable  confined  cavitated  modes  of  prop¬ 
agation  [26],  and  (5)  presented  the  experimental  evidence 
of  relativistic  and  charge-displacement  self-channeling  of 
intense  subpicosecond  ultraviolet  radiation  in  plasmas 
[27],  including  specific  comparison  with  the  results  of 
this  computational  model  [27], 

The  current  work  presents  a  full  description  of  the 
theory  of  nonlinear  propagation  of  intense  axisymmetric 
ultrashort  laser  pulses  in  cold  underdense  plasmas.  In 
this  analysis  we  use  the  term  ultrashort  to  indicate  that 
the  duration  of  the  pulses  r  satisfies  the  inequality 
rl»r»re,  with  r,  and  re  designating  the  response 
times  of  the  ions  and  electrons,  respectively.  We  consider 
both  homogeneous  plasmas  and  preformed  plasma 
columns.  In  this  study,  no  attempt  is  made  to  establish 
consistency  between  the  local  ionization  state  and  the  lo¬ 
cal  laser  intensity,  the  issue  outlined  in  Sec.  I.  For 
sufficiently-low-Z  materials,  such  as  hydrogen  (HD.  this 
calculation  would  be  unimportant,  since  full  imaximal1 
ionization  would  be  achieved  even  in  rather  low  intensity 
(~1015-1016  W/cm:)  regions  [28].  Since  the  main  re¬ 
gime  of  interest  for  this  work  involves  intensities  greater 
than  — 1018  W/cm:,  we  believe  that  the  plasma  condi¬ 
tions  we  have  chosen  for  analysis  are  sufficiently  close  to 
the  true  self-consistent  state  to  be  adequate  for  the  in¬ 
tended  scope  of  this  study. 

The  calculations  discussed  below  have  been  performed 
with  the  specific  goal  of  exploring  the  characteristic  dy¬ 
namics  and  stability  of  the  propagation,  including  partic¬ 
ularly,  the  interplay  of  the  relativistic  and  charge- 
displacement  processes.  Therefore,  in  order  to  illustrate 
this  behavior,  numerical  results  are  presented  that  por¬ 
tray  the  propagation  and  self-focusing  action  as  a  conse¬ 
quence  of  il)  the  purely  relativistic  nonlinearity  and  2' 
the  combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms.  Although  the  relativistic 
influence  cannot  be  truly  physically  separated  from  the 
charge-displacement  process,  it  has  been  examined  sepa¬ 
rately  because  this  comparison  provides  insight  on  the 
dynamics  of  the  focusing  action,  specifically  the  process 
by  which  the  charge-displacement  and  the  confined  prop¬ 
agation  in  the  electronically  cavitated  channel  develop 
Since  both  the  relativistic  effect  and  the  charge  displace- 
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ment  tend  to  locally  produce  a  reduction  in  the  plasma 
frequency  {cop )  that  is  more  significant  in  the  high- 
intensity  regions,  both  of  these  effects  perturb  the  wave 
fronts  in  a  manner  that  encourages  convergence  of  the 
wave  and  the  formation  of  localized  high-intensity  zones. 
Furthermore,  for  both  the  purely  relativistic  and  the 
combined  cases,  we  have  calculated  the  stationary  eigen- 
modes  of  the  nonlinear  Schrodinger  equation  and  show 
their  relationship  to  the  modes  of  propagation  dynami¬ 
cally  developed  in  the  plasmas.  In  particular,  it  is  shown 
that  the  charge  displacement,  especially  that  resulting  in 
dynamical  cavitation,  has  a  very  strong  effect  on  the  spa¬ 
tial  character  of  the  propagation  [26] . 

The  calculations  indicate  that  the  efficiency  of 
confinement  of  the  propagating  energy  is  potentially 
high,  namely,  that  a  large  fraction  of  the  incident  power 
can  be  trapped  in  a  channeled  mode  by  the  combined  ac¬ 
tion  of  the  relativistic  and  charge-displacement  processes. 
Furthermore,  in  the  asymptotic  regime,  it  is  found  that 
the  channel  is  characterized  spatially  by  an  intensity 
profile  and  electron-density  distribution  corresponding  to 
the  lowest  z-independent  eigenstate  of  the  nonlinear 
Schrodinger  equation.  Overall,  the  principal  conclusion 
of  this  work  is  that,  under  appropriate  conditions,  a  new 
dynamical  mode  of  stable  highly  confined  propagation 
naturally  evolves  for  the  propagation  of  sufficiently  short 
'*  »'»V  pulses  of  coherent  radiation  in  plasmas.  By 
new  dynamical  modes,  we  mean  the  self-channeling  of  the 
radiation  through  the  formation  of  stabilized  electroni¬ 
cally  cavitated  paraxial  modes,  which  result  from  the 
combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms.  Interestingly,  for  ultraviolet 
wavelengths  in  the  i200-300)-nm  range,  the  power  densi¬ 
ties  naturally  associated  with  these  environments  can  ap¬ 
proach  thermonuclear  values. 

In  Sec.  II  the  underlying  physical  concepts  are  present¬ 
ed  and  the  governing  nonlinear  Schrodinger  equation  is 
derived.  Section  III  presents  the  analysis  of  the  eigen- 
modes  of  this  equation.  Representative  results  of  numeri¬ 
cal  simulations  of  the  two-dimensional  axisymmetric 
propagation  are  given  in  Sec.  IV.  Section  IV  also  con¬ 
tains  a  comparison  with  the  purely  relativistic  case.  Fi¬ 
nally.  the  conclusions  are  summarized  in  Sec.  V. 

II.  GENERAL  CONSIDERATIONS 
Physical  model 

Several  physical  phenomena  [1-7,14-19.24-26,29] 
play  a  role  in  the  nonlinear  dynamics  governing  the  prop¬ 
agation  of  intense  coherent  radiation  under  the  condi¬ 
tions  being  examined  in  this  study.  They  are  the  follow¬ 
ing. 

'a)  The  creation  of  a  plasma  column  by  ionization  in 
the  temporally  early  region  of  the  laser  pulse. 

(bt  The  influence  of  the  spatial  variation  of  the  refrac¬ 
tive  index  arising  from  the  nonlinear  response  of  the 
dielectric  properties  of  the  medium.  Two  mechanisms 
are  related  to  the  electronic  component;  specifically,  the 
relativistic  shift  in  electron  mass  and  the  ponderomotive- 
ly  driven  electron  motion  which  tends  to  displace  elec¬ 


trons  from  the  high-intensity  zone.  For  sufficiently  short 
pulses,  only  the  electrons  are  expelled  from  the  laser 
beam  and  the  more  massive  ions,  due  to  their  substantial¬ 
ly  greater  inertia,  are  regarded  as  motionless  [17,18,25], 
A  third  mechanism  is  the  nonlinear  response  arising  from 
the  induced  dipoles  of  the  ions,  but  this  is  generally  small 
and  negligible  [25]. 

ic)  Defocusing  mechanisms,  caused  by  diffraction  from 
the  finite  aperture  and  refraction  by  the  transverse  mho- 
mogeneities  in  the  electron  density. 

id)  Dissipation  of  laser-beam  energy  by  (i)  motion  of 
the  electrons,  lii)  ionization  of  the  gas  atoms,  (iii)  genera¬ 
tion  of  harmonic  radiation,  tiv)  production  of  inverse 
bremsstrahlung,  (v)  Compton  scattering,  and  (vi)  other 
amplitudes  of  nonlinear  scattering. 

The  present  work  incorporates  four  phenomena:  1 1 ) 
the  nonlinear  response  of  the  refractive  index  of  the  plas¬ 
ma  due  to  the  relativistic  increase  in  the  mass  of  the  free 
electrons,  1 2)  the  refractive  index  variation  due  to  the  per¬ 
turbation  of  the  electron  density  by  the  ponderomotive 
force,  (3)  the  diffraction  caused  by  the  finite  aperture  of 
the  propagating  energy,  and  (4)  the  refractioi^generaied 
by  the  transversely  inhomogeneous  plasma  density  associ¬ 
ated  with  the  formation  of  a  plasma  column.  In  these 
calculations,  preformed  static  plasma  columns  were  used 
in  order  to  approximate  the  radial  distribution  of  ioniza¬ 
tion  that  is  expected,  if  the  incident  laser  pulses  were  pro¬ 
ducing  the  ionization  on  their  temporally  leading  edge. 
Finally,  the  calculations  were  performed  for  a  length  of 
propagation  that  is  much  shorter  than  the  characteristic 
length  for  dissipation  of  the  energy  of  the  pulse  through 
ionization  or  other  modes  of  energy  loss. 

It  should  be  noted  that  the  relativistic  intensities 
characteristic  of  the  phenomena  examined  in  this  work 
can  be  currently  obtained  experimentally  [30-32],  In 
particular,  the  experimental  parameters  characteristic  of 
the  ranges  that  would  apply  to  the  study  of  these  phe¬ 
nomena  are  presented  in  Table  I.  With  the  wavelength 
and  range  of  electron  densities  shown  in  Table  I,  the  plas¬ 
ma  is  always  underdense,  namely,  i  o>p /w ):  «  1 ,  where 
u>  =  2 trc  /k  is  the  angular  frequency  of  the  laser  radiation 
and  wp=t4 rre'~\e/me)]/2  is  the  customary  plasma  fre¬ 
quency. 

B.  The  propagation  equation 

Consider  the  propagation  of  an  intense  ultrashort  laser 
pulse  in  a  plasma  with  an  initially  radially  inhomogene¬ 
ous  electron  density,  described  by  the  function  /in,  so 


TABLE  I.  Experimental  parameters  characteristic  of  the 
ranges  that  would  apply  for  the  study  of  relativistic  and  charge- 
displacement  self-focusing  of  laser  pulses  in  plasmas. _ 


Peak  intensity 

/=10"-102°  W/cm: 

Pulse  length 

100-1000  fs 

Initial  focal  spot  radius 

r0  ^  1  -3  pm 

Wavelength 

k  =  0.248  jim  iKrF*  laser  light 

Target  gas  densities 

1016- 1020  cm’1 

Initial  unperturbed 

.V,  ==  10]T- 1021  cm'1 

electron  density 
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that  S{,Q)  —  Staf(r),  max/(r)=l.  We  denote  the  vector 
and  the  scalar  potentials  of  the  electromagnetic  fields  as 
A  and  <&,  respectively,  and  the  corresponding  electric 
field  as  E.  Let  the  momentum  of  the  electrons,  the 
current  density,  and  the  charge  density  be  denoted  as  Pe, 
j,  and  p,  respectively.  We  assume,  for  the  short  timescale 
of  interest,  that  the  ions  are  inertially  frozen  in  space 


[17,18,24,25].  Then. 

[2  A  =  c  'V——  (4rr/c)j  ,  (1) 

3r 

V26  —  —4 ?rp  ,  ( 2) 

<VA)  =  0,  (3) 

3  —  _  1 3  A  . 

—  -Hv-V)  p,  =  —  e  ~c  — - V<b 

3 1  e  *  3/ 

+  c''[v(XlVX  A)]  ,  (4) 

j  =  —  eSt ve ,  p=e(.\e°'-.\')  ,  (5) 

v,  =p,/m,,  me=me0y, 

(6) 


Ym[  1+  p, i2/( m#i0c  )2]l/2  . 


both  space  and  time  is  much  greater  than  both  the  plas¬ 
ma  (2irc/<up0)  and  electromagnetic  (3.)  wavelengths. 
This  leads  to  the  validity  of  the  inequalities 


'  !’  ilr  «!fall*VI 


where  k1=kQ—kp,  k0=ca/c,  and  kp=<jj)Q/c.  We  use 
the  notation  d>po  =  4ire2iVt0/m(0  specifically  to  denote 
the  unperturbed  plasma  frequency.  With  the  assump¬ 
tions  and  approximations  stated  above,  Eqs.  (1M8)  be¬ 
come 


□  A=*;.v,r-‘A  , 

■10) 

V-6  =  kl{Se-f(r))  , 

ill) 

Vf  <b~y  )  =  0  , 

<12) 

/=(  1+  p,  2)1/2,  pf=  A  . 

1 13) 

The  term  c  _l3V<4/3r  is  omitted  in  Eq.  (10),  since  y  and 
6,  according  to  Eqs.  (9),  (12),  and  (13),  do  not  have  a 
high-frequency  dependence.  Equations  (11)  and  (12)  re¬ 
sult  in  the  expression  for  the  electron  density, 

:V,=max[0,/(r)  +  Jfc-2V2r]  .  (14) 


In  the  set  of  statements  above,  Eqs.  (1)  and  (2),  are  the 
Maxwell  equations,  Eq.  (3)  is  the  Coulomb  gauge  condi¬ 
tion,  Eq.  (4)  is  the  equation  of  motion  of  the  electrons, 
Eqs.  (5)  are  the  definition  of  the  current  density  and  the 
charge  density  (3  C  is  the  initial  charge  density,  while 
is  the  dynamical  charge  density),  and  Eqs.  (6)  represent 
the  relativistic  relation  between  the  velocity  v,  and  the 
momentum  of  the  electrons.  In  Eqs.  (6),  me  Q  is  the  rest 
electron  mass,  Ci  =  V2  — c  -23:/3f2. 

It  is  convenient  to  normalize  the  values  in  the  equa¬ 
tions  presented  above  as  follows: 

\  =  <e/me0c:)  A,  $  =  ( e/mf  0c2)i, 

£  =  te/me0c:  lE  ,  P^=pf/me0c,  (7) 

V,*v,  /c.  .V,  =.V,  /\eQ  , 


The  logical  function  max(0,  )  provides  for  the  physical¬ 
ly  obvious  and  necessary  condition  Ne  >  0.  The  analo¬ 
gous  expression  for  the  electron  density  has  been  previ¬ 
ously  derived  by  Sun  et  al.  [18]  for  the  case  /( r)  =  1.  It 
should  be  noted  that  Eq.  (12)  states  the  condition  for  the 
balance  of  the  ponderomotive  and  the  electrostatic  forces 
for  the  relativistic  case.  Through  combining  Eqs.  '9>. 
1 10),  and  (14),  we  establish  the  equation  for  the  slowly 
varying  complex  amplitude  of  the  vector  potential 
a[r,z,t)  as 

J_  JL  +  _3_ 
vg  3r  3 z 

-+■  —  ( V2a  +kh  l  —  y  ~  lmax[0,/(rH-/lc~2V2}'  ]  \a  ) 

Ik  *■  y  p 

=  0  .  15 


a 


with  the  understanding  that,  henceforth,  the  tilde  sign 
will  be  suppressed.  Using  the  relations  <p,  V)p, 
=  V  p,  2/2-p,  Xi  V  Xp, )  and  Vy  =  V  p,i2/2y,  Eq.  (4) 
becomes 

-^-ip,  -  A)— v^XfVxlp,  —  A)]  =  V(<$  —  y)  .  (8) 

3f 

In  the  limit  r»2ir /wp,  the  expression  pe  —  A  is  approxi¬ 
mately  valid.  As  we  shall  see  below,  this  condition  means 
that  the  electron  response  can  be  regarded  as  adiabatic. 
Furthermore,  assuming  the  vector  potential  to  be  circu¬ 
larly  polarized,  we  write 

A(r,z.r)  =  4 1  (et  'a  ( r.z.f  )exp(i< cot  —  Icz)]  +  c.c.  j  . 

i9) 

Consistent  with  the  statements  made  above,  we  em¬ 
phasize  the  use  of  the  assumption  that  the  pulse  length  in 


In  accord  with  the  previous  assumptions  used  in  deriving 
Eq.  (15),  we  have  neglected  the  second  z  and  t  derivatives 
in  this  expression.  In  Eq.  (15),  ug  =cel0/Z  is  the  group  ve¬ 
locity  of  light  in  unperturbed  plasma,  c0=  1  —  ( 
is  its  corresponding  dielectric  constant,  ar.q 
y—(  1+  ai2)1'2.  It  should  also  be  noted  that  the  electro 
field  vector  E  and  the  vector  potential  A  are  relatej 
through  the  approximate  relation 


E(r,z,r)  =  —  ik0  A(r,z,t)  . 

The  calculation  of  the  propagating  wave  form  for 
central  high-intensity  region  of  the  pulse  is  accompli': 
by  considering  the  solutions  of  Eqs.  (15)  along  its  char :. 
terisucs.  Changing  the  variable  t  to  q  =t  —  z/vg,  Eq 
becomes 
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~a  +  ~(V?a  +  kp\  1  -y  'max[0,/(r)  +  Ic~2V2y]|a) 

«  2*  ' 

=0,  (16) 

for  which  we  seek  solutions  of  the  form  <j  =  const.  We  re¬ 
call  that  certain  solutions  of  Eq.  (16)  for  the  special  case 
fir  i=l  have  already  been  presented  [18,19,33]. 

The  results  of  the  computations  presented  in  this  work 
are  given  with  respect  to  the  variables  i  q,r,z). 
Specifically,  the  data  presented  in  the  figures  below  illus¬ 
trate  the  propagation  of  the  radiation  along  the  z  axis  for 
<7=const.  Equation  1 16)  describes  the  two-dimensional 
ir.z)  dynamics  of  propagation  for  coherent  circularly  po¬ 
larized  radiation  in  plasmas  having  an  inhomogeneous 
transverse  i  r-dependent)  distribution  of  electron  density. 
The  basic  physical  phenomena  explicitly  embodied  in 
these  equations  have  been  outlined  in  Sec.  II  A. 

III.  STATIONARY  SELF-LOCALIZED  MODES 
OF  PROPAGATION 

Equation  M6),  in  the  case  of  homogeneous  plasmas 
[/( r)  =  1  ],  has  axisymmetnc  partial  solutions  of  the  form 

a  (r,z>  =  Uf  „t  k^nexpft  [k^/2k)is  —  1  )z]  .  (17) 

where  s  is  a  real-valued  dimensionless  parameter  and  the 
real-valued  function  U%  „  obeys  the  ordinary  differential 
equation  [18] 

.  '18) 

In  this  discussion  the  dimensionless  argument  p  —  kpr  is 
used  and  V2  in  Eq.  1 18)  denotes  the  Laplacian  with  the 
derivatives  designated  with  respect  to  this  variable  (p>. 

The  nonlinear  term  in  Eq.  1 18)  is 

.V. ,  =maxi0. 1  n  i  ,  '19) 

-ti-t’,:,.' : . 

The  natural  boundary  conditions  for  Eq.  (18)  are 
di’in  i0) 

- — - =0.  (.'.„(  *)=0.  '20) 

dp 

The  first  condition  assures  axial  symmetry  of  the  solu¬ 
tions.  while  the  second  is  necessary  for  the  property  of 
finite  energy  for  the  solutions.  We  designate  these  solu¬ 
tions  as  "eigenmodes"  and  their  significance  is  explored 
below. 

Consider  initially,  however,  the  eigenmodes  corre¬ 
sponding  to  the  purely  relativistic  case,  functions  which 
are  obtained  from  Eqs.  '18)  and  (19)  by  neglecting  the 
term  V:y  Since,  we  are  explicitly  eliminating  the  pon- 
deromotive  potential  and  its  influence  on  the  motion  of 
the  electrons,  in  this  situation  the  electron  density  will 
not  be  self-consistent.  In  this  restricted  case,  Eq.  1 18)  can 
be  rewritten  in  the  form 

d_  '  ,  | 

dp  [  ;  j  dp 


where  V(U,n,s)=(s/2)Uj;n-(l  +  U*n)l/2+\. 

The  second  boundary  condition,  stated  in  Eq.  (20),  en¬ 
sures  finiteness  of  energy  only  in  the  situation  where  the 
dynamical  system  defined  by  the  last  two  expressions  has 
three  rest  points,  one  of  them  being  zero,  that  is  for  the 
case  0<r  <1.  Exactly  as  in  the  classic  theory  of  cubic 
media  [34,35]  for  values  of  the  parameter  s  belonging  to 
this  interval,  the  localized  eigenmodes  constitute  a  count¬ 
able  set  that  is  ordered  by  the  number  of  zeros- n  in  each 
of  its  members  for  finite  values  of  p.  Figure  1  illustrates 
the  first  four  eigenmodes  for  s  =0. 95. 

Consider  now  the  eigenmode  corresponding  to  the 
unrestricted  situation  that  includes  both  the  relativistic 
and  charge-displacement  mechanisms,  a  case  which  is 
found  to  be  fully  analogous.  The  eigenmodes  of  the  equa¬ 
tion  describing  the  combined  relativistic  and  charge- 
displacement  processes  for  the  interval  0  <  s  <  1  also 
make  a  countable  set  that  can  be  ordered  in  the  same  way 
as  the  one  associated  with  the  purely  relativistic  case. 
The  zeroth  (lowest)  eigenmode  £/t0(p)  with  s  =0.95. 
along  with  the  corresponding  electron  density  .V^p1, 
given  by  the  second  formula  in  Eq.  (19),  are  depicted  in 
Fig.  2.  These  two  functions  have  been  developed  in  ear¬ 
lier  work  [18].  The  first  and  second  relativistic  and 
charge-displacement  eigenmodes  for  the  value  of  the  pa¬ 
rameter  s  =0.95  and  the  corresponding  electron-density 
eigenmodes  <VS  ,(p)  and  \s  2lp>  are  depicted  in  Figs.  3  and 
4.  respectively.  These  higher  eigenmodes  exhibit  the  im¬ 
portant  feature  that  cavitation  may  occur  in  the  electron¬ 
ic  component  of  plasma  [see  Fig.  (4)],  even  if  it  does  not 
occur  in  the  case  of  the  lower  eigenmodes  [see  Figs.  2  and 

3]- 

A  prior  study  [18]  treated  the  (r,z)  dynamics  ot  the 
evolution  of  the  perturbed  lowest  eigenmodes  for  the  case 
in  which  no  cavitation  of  the  electronic  component  of 
plasma  occurs.  The  lowest  eigenmodes  were  perturbed 
by  multiplying  them  by  constants  in  the  neighborhood  of 
unity.  In  this  case,  the  known  regime  involving  oscillato¬ 
ry  propagation  was  observed.  However,  the  instability  of 
the  numerical  method  applied  in  that  work  [18]  prevent¬ 
ed  the  performance  of  similar  computations  for  cases  in¬ 
volving  cavitation.  Our  calculations  show  that  the  analo¬ 
gously  perturbed  lowest  eigenmodes  Fs0(p)  also  propa¬ 
gate  in  the  same  oscillatory  fashion  when  electronic  cavi¬ 
tation  occurs,  the  corresponding  density  of  which  is  given 
by  o'P1-  I{  1S  also  found  that  the  same  behavior  devei- 


FIG  1.  Stationary  axially  symmetric  eigenmode' 
s  =0.95  corresponding  to  purely  relativistic  self-focusmc 


y<  v 


,  :  -  -p 


dU ,, 

dp 
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FIG.  2.  The  zeroth  'lowest)  axially  symmetric  stationary 
eigenmode  with  s  =0.95  corresponding  to  relativistic  and 
charge-displacement  self-focusing:  the  normalized  held  ampli¬ 
tude  distribution  £/,,0(p)  and  the  normalized  electron  density 
distnbution  .V,  0(p). 


ops  for  the  perturbed  first  eigenmodes,  which  have  one  or 
two  cavitated  channels  in  the  electronic  density  described 
by  -Vt  i (p).  The  results  of  these  simulations  can  be  re¬ 
garded  as  evidence  of  the  (r,z)  stability  of  axially  sym¬ 
metric  lowest  and  first  eigenmodes  against  small  pertur¬ 
bations  in  the  amplitude.  It  should  be  noted,  however, 
that  these  higher  eigenmodes  are  presumably  unstable 
against  small  azimuthal  perturbations,  the  nature  of 
which  violates  the  assumed  axial  symmetry  of  the  distri¬ 
butions. 

The  important  consequence  of  the  character  of  the 
solutions  discussed  above  is  the  fact  that  the  power 

P,~2ir  fju20(p)pdp 

contained  in  the  intensity  distribution  corresponding  to 
the  lowest  eigenmode  U!  0(p),  unlike  the  case  of  a  cubic 
medium,  depends  on  the  parameter  s.  Namely,  this 
power  decreases  as  s  increases. 

The  infimum  of  Ps  by  s,  in  the  internal  0<s  <1,  is 
called  the  critical  power  (Pcr)  of  the  relativistic  and 
charge-displacement  self-focusing.  This  power  equals  the 
critical  power  of  the  purely  relativistic  self-focusing 
[18,36]  and 


/\-r 


=  inf  P,  =  lim  P,  . 

0  <  S  <  1  5  —  1  “0 


y,,(c) 


i 


FIG.  3.  The  first  axially  symmetric  stationary  eigenmode 
with  j=0.95  corresponding  to  relativistic  and  charge- 
displacement  propagation:  the  normalized  amplitude  distribu¬ 
tion  t/,,(p)  and  the  normalized  electron-density  distribution 


FIG.  4.  The  second  axially  symmetric  stationary  eigenmode 
with  s  =0.95  corresponding  to  relativistic  and  charge- 
displacement  propagation:  the  normalized  amplitude  distnbu¬ 
tion  £/,,(p)  and  the  normalized  electron-density  distnbution 


The  explicit  evaluation  of  this  critical  value  is  presented 
below.  As  UsQ(p)~* 0,  Nj  0(p)— *1  for  s  — 1— 0;  we  have 
[18]  therefore, 

r~o\o ~  1  =  Ns0(  l  +  U,20)-l/2-l 

=  ( 1  +  U20  /2 )  ~ 1  —  1  =r  -Us:0/  2. 

j  — 1-0  '  s  — 1-0 


Furthermore, 


tf,o  <P>  *  W. 

j  — 1-0 

where  U0(p)  is  the  positive,  monotonically  decreasing 
(reaching  zero  for  no  finite  value  of  p)  solution  to  the 
boundary-value  problem  [18,36] 

V2LU0-€U0  +  {Ul=0  , 

dU0 

— —  (0>=0,  U0(  *  )=0  , 

dp 

with  €—  1  —  s.  A  change  of  the  variables,  which  is  a  stan¬ 
dard  procedure  for  treating  the  case  of  the  cubic  non¬ 
linearity  [37],  shows  that  [36] 

U0(p)  =  (2e)w2g0(ewlp)  , 

where  g0  is  the  customary  Townes  mode,  i.e.,  the  posi¬ 
tive,  monotonically  vanishing  solution  to  the  following 
boundary-value  problem: 

V2.go~go+go=0  ’ 
dg0 

— ~  ( 0 )  =  0,  g0(  so  1=0  . 
dp 

From  the  definition  of  Pct  and  the  relation  between  C/0» p  1 
and  g0(p),  it  follows  that 


P 


cr 


=  inf  P. 

0  <  s  <  1 


=  lim  P, 

s  — 1-0 


=  2n  f  * U\(p)pdp 
J  o 

=4w  f  * gl(p)p dp  =  2P c  , 

y  0 

where  Pirc  =  2rr  J*0xgo*P  WP  is  the  critical  power  of  the 
Kerr  self-focusing  in  cubic  media  [34,37], 

Using  the  resulting  relation  between  the  normalized 
values  of  the  critical  power  of  the  relativistic  and  charge- 
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displacement  self-focusing  PCT  and  the  critical  power  of 
the  Kerr  self-focusing  Pcr  c  we  have  the  elementary  result 
that 

P  =  ip 
*  cx  *"*  cr.c 

Finally,  by  calculating  g,yp )  numerically,  we  find  the  nor¬ 
malized  value  of  P  T  with  high  precision  [36]  to  be 

P„  =2 P,r  c  =4rr  f  ^g.ypip  dp  ^23.4018  . 

Hence,  we  establish  the  expression  for  the  critical  power 
[36]  as 

P;  r.1  =  lm-0c5/e:)  fJg'o{p)pdp(co/wpQ)2 

^  1.6198X  [O^lco/cOpQp  W  . 

The  constant  factor  involved  in  this  statement  improves 
on  that  given  previously  [18]. 

The  relation  between  Lrnip)  and  g0(p )  enables  the  de¬ 
velopment  of  the  asymptotic  expressions  (in  the  case 
0)  for  both  the  peak  value  of  the  amplitude  of  the  sta¬ 
tionary  eigenmode  and  for  the  radius  of  the  eigenmode  by 
using  the  appropriate  characteristics  of  the  Townes 
mode.  The  results  are 

=s< 2e)‘  :g0l0) 

and 

r,  -  e  ‘ • 

fore— 0. 

Special  computations  were  performed  in  order  to 
determine  how  close  the  values  for  C'!0tO)  and 
C'0(0)  =  <  2e)1,  :g0(0)  are  in  the  case  where  e«l.  The 
calcuiat’ons  show  that  U.  0(0)  —  C/0( 0 ) I  / f/0< 0 )  is 
2.026X  10~?  and  2.096X  lO"0  for  e=  10~5  and  10"6,  re¬ 
spectively. 

IV.  THE  TWO-DIMENSIONAL  CASE 

It  is  convenient  to  treat  Eq.  (16)  numerically  using  the 
normalized  coordinates 

r,  =r/r.-t  .  -z/{2krl)  , 

u(  r] 1  =£»,,'  'a  (  r,z)  , 

where  t,  is  a  characteristic  radius  of  the  initial  intensity 
profile  and  a0  =  max  a1  r, 0)  for  the  selected  value  of  q. 
For  simplicity,  we  put  /,(r,  )=/(r)  and  omit  below  the 
subscript  l  of  r,  and  c,  for  brevity.  Thus,  the  mathemat¬ 
ical  statements  can  be  expressed  in  the  following  set  of 
equations: 


4^+-iV2u-iF(/,<r),  u  :)u  =0,  z  >0 
az 

i21l 

u  (r,0)  =  u0(r),  max  u0.  =  l 

r 

(22) 

—  (0,z)=0,  ui  x,z)  =  0  . 

‘23) 

dr 

The  nonlinear  term  F  is  the  real-valued  operator 


*Vi.£)=a,{l-(l+a2£r1/2 

x  max[0,  f \(r)+a  j"1  V2(  1  a ) 1  :  '  , 

24 1 

and  the  dimensionless  parameters  a{,  az  are  defined  as 

a,=(r0A:p)2,  a2=IQ/Ir,  I0  =  m20(i)2cial/<4ire: '  , 

'25) 

where  I0  is  the  peak  intensity  at  the  entrance  plane 
<z  =  0)  of  the  medium.  The  parameter 
I,  =  m}0arcl /{4irel)  is  known  as  the  relativistic  intensity 

[15]- 

The  ratio  of  the  power  of  the  beam  P0  and  the  critical 
power  of  the  relativistic  and  charge-displacement  self- 
focusing  Pcr,  defined  in  Sec.  Ill,  is  an  important  parame¬ 
ter  characterizing  Eqs.  (21)  — (24).  In  the  notations  of  the 
present  section,  the  value  of  P0/PCT  can  be  expressed  in 
the  following  way  [36]: 

P0/Pcr=(aia-,/B)  f  *  u0{r)\2r  dr  , 

J  o 

with  dimensionless  constant  B  given  by  [36] 

B—2 fQXgl(p)pdp^3.1245\  . 

When  ,u~  initial  transverse-intensity  distribution  is 
Gaussian,  namely,  iu0(r)|2  =  exp(  —  r2),  we  have 

P0/Pzl=a\a-,/(2B)  . 

We  note  that  in  several  other  studies,  the  critical  power 
of  the  -elativistic  self-focusing  is  defined  alternatively,  ba¬ 
sically  with  B  =4.0.  In  the  present  discussion,  the  ex¬ 
pression  for  Pq / PCT  involves  a  value  of  B  <4.0.  so  that 
the  corresponding  ratio  Pq/P^  is  raised.  Therefore,  a 
pulse  having  a  Gaussian  initial  transverse-intensity  distri¬ 
bution  and  a  flat  initial  wave  front  with  the  parameters 
a,  =248.6192  and  a2  =0.031,  undergoes  self-focusing.  Ir 
this  case,  the  ratio  P0/Pcr,  as  defined  above,  evaluates  to 
1.0347,  namely,  P0>Pcr.  However,  using  the  value 
B  =4.0  would  give  P0  < /\r,  a  statement  contradicting 
the  results  of  the  computations. 

A.  The  initial  conditions 

In  this  section,  we  examine  the  self-focusing  of 
coherent  radiation  for  pulses  having  Gaussian  or  hyper- 
Gaussian  transverse  and  longitudinal  intensity  distribu¬ 
tions  [24]  of  the  form 

[[!^I0lr,t)  =  Imtxp{-(t/r)S'-(r/r0)'':]  , 

AT,  >2,  <V:  -  2  2b 

with  r  and  f  being  dimensional.  We  assume  [24]  that 
Im  =  /,»/*  =  1016  Wcm"2  with  I*  designating  the  ap¬ 
proximate  value  of  the  threshold  for  rapid  nonlinear  ion¬ 
ization  [1,2].  The  spatial  amplitude  distribution  of  t he 
cident  radiation,  defined  by  Eq.  (26)  for  the  case  ot  a  o 
incident  phase  front,  is  of  the  form 
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u0(r)=exp(  — r Vz/2)  ,  iV2> 2.  (27) 

The  dimensionless  parameter  a2  in  Eq.  (25),  correspond¬ 
ing  to  the  incident  pulse  intensity  on  the  axis  (r=  0)  of 
the  beam  /0(r),  is 

a2  =  /0i t )//,  =  < Im  //,  )exp[  — ( t/r)'S'  ]  .  (28) 

The  transverse  profile  of  the  plasma  column,  created  by 
the  temporally  leading  edge  of  the  pulse,  is  simulated  by 
the  hyper-Gaussian  function 

/(r)  =  exp[  —  (r/r.  )'  5  ],  ,\}>  2,  (29) 

where  r  is  the  radial  coordinate.  The  aperture  of  the 
plasma  column  r0  can  be  estimated  [24]  with  the  use  of 
the  relation 

/olr.,f0)=/ouo,exP(_(r./ro)'V;]=:f*  ■  * 30) 

For  example,  in  the  case  of  the  Gaussian  transverse  inten¬ 
sity  distribution  (.V,  =2),  the  aperture  of  the  plasma 
column  for  /*  =  1016  Wcm";,  /,^0.45X  1020  Wcm~:, 
and  /0(f0)  =  (0. 1 )/,,  gives  r.  =  2.47r0.  In  this  situation, 
the  homogeneous  plasma  approximation  /(r)=  1  is  valid 
[24].  In  contrast,  it  follows  from  Eq.  (30)  for  plateaulike 
incident-transverse-intensity  distributions  that  the  aper¬ 
ture  r9  of  the  simulated  plasma  column  tends  to  the 
value  of  the  beam  aperture  r0.  Therefore,  in  the  example 
given  above  for  .V,  =8,  we  have  r.  =  1.25r0.  Thus,  de- 
focusing  of  the  beam,  which  is  significant  because  of  the 
near  coincidence  of  the  apertures  of  the  laser  beam  and 
the  plasma  column,  must  be  taken  into  account  when  the 
evolution  of  a  beam  with  a  plateaulike  incident- 
transverse-intensity  distribution  is  studied  [24]. 


where 

4(r,g)  =  f*F(r,i))di} 
v  0 

=  a1|£-(2/a,)/l(r)[(l+a2£)l/2-l]|  .  (33) 

Figure  5  illustrates  the  nature  of  the  calculated  solu¬ 
tions  corresponding  to  this  purely  relativistic  case.  The 
parameters  of  the  incident  radiation  and  plasma,  for  the 
examples  presented  in  Fig.  5,  are  A.  =  0.248  pm, 
/,  =0.45  X  1020  W  cm-2,  /0  =  }/,  =  3X1019  Wcm':, 
r0  =  3  pm,  and  .Vf  0=7.5X  1020  cm-3.  The  correspond¬ 
ing  values  of  the  associated  dimensionless  parameters  are 
a  |  ^  2.486  192  X  102  and  a,  =  y. 

Figure  5(a)  presents  the  result  for  purely  relativistic 
propagation  in  a  homogeneous  plasma  along  the  r  axis 
for  a  pulse  having  an  incident  Gaussian  transverse  inten¬ 
sity  distribution  and  a  flat  wave  front  [#2=2  in  Eq.  (27)] 
for  the  value  of  q  defined  in  Sec.  II B,  corresponding  to 
Figure  5(b)  presents  the  analogous  graph  for  a 
plateau  like  incident  transverse  intensity  distribution 
[.V,  =  8  in  Eq.  (27)].  We  conclude  from  these  results  that 
the  solution  is  critically  dependent  upon  the  initial  condi¬ 
tion  represented  by  the  detailed  character  of  the  incident 
transverse  intensity  distribution. 

The  ratio  of  the  beam  power  ( P0 )  to  the  critical  power 
( P.r )  of  the  relativistic  self-focusing  for  the  given  values 
of  the  parameters  a,  and  a2  yields  F0/Pcr  =  22.252  in  the 
case  of  the  Gaussian  initial  transverse  intensity  distribu¬ 
tion  [,V2=2  in  Eq.  (27),  Fig.  5(a)]  and  P0/PCT  —  20. 168  in 
the  latter  example  of  the  plateaulike  [.V2  =  8  in  Eq.  (27), 
Fig.  5(b)]  initial  transverse  intensity  distribution. 

Note  that  for  the  values  of  the  parameters  a,. a 2  and 
incident  wave  forms  studied,  we  have  from  Eq.  (32) 
P2  <  0.  In  this  situation,  the  following  inequality  is  valid: 


B.  Relativistic  self-focusing 


Consider  the  two-dimensional  ( r,z)  solutions  of  the  sys¬ 
tem  of  equations  embodied  in  Eqs.  (21),  (22),  and  (23)  for 
the  purely  relativistic  nonlinear  term,  the  form  of  which 
can  be  obtained  from  Eq.  (24)  by  disregarding  the  term 
involving  V2.  In  this  case,  we  find 

F(/, .*>=<!, [1  -/,<r)(l+<j;sri  ;]  . 


The  relativistic  self-focusing  mechanism  prevails  over  the 
charge-displacement  mechanism  outside  of  the  focal  spot 
under  the  conditions  a(  »I,  a2  —  1-  (Discussion  of  the 
conditions  for  the  prevailing  of  the  relativistic  self- 
focusing  can  also  be  found  in  Ref.  [33]).  The  situation 
represented  by  Eqs.  '21),  22),  and  (23),  with  the  above 
nonlinear  term,  describes  self-focusing  with  a  nondissipa- 
tive  saturation  of  the  nonlinearity.  The  properties  of  the 
solutions  for  this  case  depend  essentially  on  the  values  of 
two  conserved  integrals  given  specifically  by 

P,  =  f  u  (  r)  :r  dr  (31) 

J  o 

and 


du 

dr 


—  <t>(r,  u  " ) 


r  dr 


(32) 


max  u  (r,z)\z  >  (4/a,a2  )iP2!  /P,  , 

r 


(34) 


namely,  with  respect  to  the  radial  coordinate  r,  the  max¬ 
imum  beam  intensity  has  a  positive  lower  bound  indepen¬ 
dent  of  z.  This  conclusion  can  be  established  in  the  same 
way  that  Zakharov,  Sobolev,  and  Synakh  [37]  demon¬ 
strated  the  analogous  result  in  their  earlier  work  on  self- 
focusing.  Therefore,  a  powerful  relativistic  beam,  which 
self-focuses  in  a  homogeneous  nonabsorbing  plasma 
p=0,  Ref.  [24]),  results  in  a  field  distribution  represent¬ 
ing  a  pulsing  waveguide  when  P:  <0.  Figures  5ia)  and 
5(b)  explicitly  illustrate  the  formation  of  such  a  regime. 
These  pulsing  waveguides  consist  of  alternations  of  ring 
structures  and  focal  spots  on  the  axis  of  the  beam.  The 
power  confined  in  these  complex  modes  represents  ap¬ 
proximately  50%  and  90%  of  the  total  incident  power  for 
the  cases  depicted  in  Figs.  5(a)  and  5(b),  respectively.  It 
should  be  noted  that  oscillating  periodic  solutions  of  this 
type  have  also  been  presented  by  other  workers  [17], 

In  the  case  P2  >0,  an  inequality  comparable  to  Eq.  < 34) 
cannot  be  obtained.  Following  the  method  of  Zakharov. 
Sobolev,  and  Synakh  [37],  we  simulate  the  case  P:  iO  by 
considering  a  beam  initially  focused  or  defocused  by  a 
lens  at  the  entrance  plane  (z  =0)  of  the  medium.  Let  the 
Reusing  (or  defocusing)  length  of  this  lens  be  R  =kr\R , 
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FIG  5.  Purely  relativistic  propagation  with  /0  =  3xlO'‘ 
W  cm;.  r  =  3  urn. =0. 248  jim.  and  V,  j -7.5  x  lO:o  cm  .  al 
The  formation  of  the  pulsing  waveguide  regime  in  the  case  of 
the  relativistic  self-focusing  of  a  pulse  with  a  flat  incident  wave 
front:  Gaussian  initial  transverse-intensity  distribution  [.V.  =2 
in  Eq.  22!],  homogeneous  plasma.  (bi  The  formation  of  the 
pulsing  waveguide  regime  in  the  case  of  the  relativistic  self- 
focusing  of  a  pulse  with  a  flat  incident  wave  front;  hyper- 
Gaussian  initial  transverse-intensity  distribution  [.V,=8  in  Eq. 
2’i],  homogeneous  plasma.  ■  c )  The  single-focus  regime  in  the 
case  of  the  relativistic  self-focusing  of  an  initiallv  focused  pulse 
R.  R.  .,2  ,n  Eq.  ■  35']:  Gaussian  incident-transverse-intensity 
distribution  V.  =2  in  Eq.  ‘ 35 l],  homogeneous  plasma,  d)  The 
formation  of  the  quasistabilized  regime  in  the  case  of  the  rela¬ 
tivistic  self-focusing  of  a  pulse  with  a  flat  incident  wave  front  in 
a  plasma  column  f.V,  =8.  r,  =r  in  Eq.  1 29)];  hyper-Gaussian 
initial  transverse-intense  distribution  [.V.  =8  in  Eq.  1 27)]. 


(Rf  dimensionless).  The  condition  Rf> 0  signifies  that 
the  pulse  is  initially  focused  and  Rf  <  0  indicates  that  it  is 
initially  defocused.  Then,  the  corresponding  initial  con¬ 
dition  can  be  written  as 

u0l  r)  =  exp[  —  r  1  .'2  +  irz/(2Rf  )],  .V,  >  2  .  5  5 

Furthermore,  let  P2=  0  for  Rr]  =  Rf  0.  We  note  that  the 
case  P,  >  0  corresponds  physically  to  a  high  degree  of  the 
initial  focusing  or  defocusing:  Rl  0>  Rr  >0.  Moreover, 
when  R  fQ  <  Rr  <*,/>,<  0. 

Figures  5(a)  and  5(c)  display  the  calculated  intensity 
distributions  for  the  case  of  relativistic  self-focusing  in 
homogeneous  plasma  of  beams  with  Gaussian  mcident- 
transverse-intensity  distributions  and  values  of  R ,  given 
by  R ,  =  x  (flat  wave  front)  and  Rf%0/2,  respectively.  A 
more  detailed  study  of  the  dynamics  of  the  transition 
from  the  pulsing  waveguide  regime  [see  Fig.  5 ( a > j ,  to  the 
single-focus  regime  [see  Fig.  5(c)  representing  the  propa¬ 
gation  of  an  initially  sufficiently  sharply  focused  beam  ' 
can  be  found  in  Ref.  [38].  Importantly,  the  computations 
show  that  the  value  P:  =0  is  not  the  threshold  separating 
these  two  regimes  of  relativistic  self-focusing.  This  tran¬ 
sition  occurs  as  the  first  focus  gains  power  and  is  shifted 
closer  to  the  entrance  of  the  medium  iz  =0>,  while  the 
remaining  foci  are  shifted  in  the  opposite  direction  and. 
in  the  limit  of  large  displacement,  become  diffused.  Note 
that  single-focus  regimes  of  propagation  have  also  been 
observed  by  Sprangle,  Tang,  and  Esarey  [17]. 

Defocused  beams  can  evolve  in  a  different  fashion  Ini¬ 
tially  sufficiently  sharply  defocused  pulses  R 
Rf  %Rr 0/2,  P:>  0)  propagating  in  the  relativistic  re¬ 
gime  monotonically  diffuse  on  the  radial  periphery  anj 
do  not  exhibit  the  phenomenon  of  self-focusing. 

The  detailed  spatial  character  of  the  plasma  column 
can  have  a  strong  influence  on  the  evolution  of  the  propa¬ 
gation.  Figure  5(d)  illustrates  the  relativistic  propagat:  n 
in  a  plasma  column  along  the  z  axis  for  a  pulse  corre¬ 
sponding  to  .V,=8  in  Eq.  (27).  The  transverse  profile 
the  plasma  column  is  given  by  fir)  as  defined  by  Eq  I  • 
with  ,V3  =  8  and  r,  =r0.  The  comparison  between  F 
5(b)  and  5(d)  demonstrates  that  the  defocusing  <t 
hyper-Gaussian  beam  in  a  plasma  column,  with  an  apvr 
ture  close  to  the  radius  of  the  beam,  fundamentally  ait.  ■  - 
the  spatial  dynamics  of  the  propagating  energy  [24  1 

defocusing  causes  a  fraction  of  the  beam  to  spread  aw.r. 
from  the  column,  but  the  remaining  energy  of  the  be  i- 
resolves  into  a  state  that  balances  the  relativistic  • 
focusing,  defocusing,  and  diffraction.  The  power  trarr  : 
in  this  quasistabilized  state  so  formed  is  approximu 
25%  of  the  incident  power. 

Previous  analytical  estimates  [38]  have  shown  that 
relativistic  self-focusing  length  is  minimal  when  the 
cident  pulse  intensity  on  the  axis  of  the  beam  sati^fic- 
condition  /0(t)  =  2/r.  Moreover,  specific  compui.r 
also  show  that  this  inference  remains  valid  for  the  . 
the  relativistic  self-focusing  of  pulses  with  fla; 
phase  fronts  and  Gaussian  incident-transverse- 
distnbutions  in  homogeneous  plasmas. 

The  locus  of  the  first  focus  for  the  case  of  the  - 
tic  self-focusing  has  been  presented  previously  3  • 
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that  study,  the  beam  initially  had  both  a  flat  phase  front 
and  a  Gaussian  transverse-intensity  distribution.  The  pa¬ 
rameters  were  the  same  as  those  applying  to  Fig.  5(a),  and 
the  duration  of  the  pulse  corresponded  to  r=0. 5X  10-13 
s.  The  minimal  z  of  the  locus  is  reached  when  I0(t)  =  2Ir. 
If  the  maximum  intensity  on  the  axis  of  the  beam  is  such 
that  >21,,  then  the  trajectory  of  the  first  focus  exe¬ 
cutes  a  path  that  has  three  reversal  points.  Two  of  them, 
corresponding  to  an  identical  value  of  z,  are  due  to  the 
condition  /0(r>  =  2/,.  The  central  one,  corresponding  to 
a  greater  value  of  z,  arises  when  IQ(t)=I„.  In  principle, 
this  behavior  makes  it  possible  to  distinguish  the  process 
of  relativistic  self-focusing  from  that  associated  with  the 
Kerr  nonlinearity,  since  the  latter  results  in  a  locus  of  the 
first  focus  having  only  a  single  reversal  point. 

C.  Relativistic  and  charge-displacement  self-channeling 

This  section  is  devoted  to  the  description  of  the  self¬ 
channeling  occurring  when  both  the  relativistic  and 
charge-displacement  mechanisms  are  included  in  the  in¬ 
teraction.  This  situation  is  described  by  Eqs.  (21)-i24) 
and  includes  the  important  nonlinear  term  involving 
appearing  in  Eq.  (24).  The  numerical  results  show  that 
sufficiently  intense  short  duration  l  r,  »r»r(,)  axisym- 
metric  pulses  readily  undergo  self-channeling  in  plasmas 
over  a  rather  broad  range  of  conditions.  Moreover,  it  is 
found  that  a  large  fraction  of  the  total  incident  power  of 
the  beam  can  be  trapped  in  a  stabilized  mode  confined  to 
the  axis  of  propagation. 

A  specific  example  is  informative  in  representing  the 
general  behavior  exhibited  by  the  propagation  in  the  re¬ 
gime  in  which  the  influence  of  relativity  and  charge  dis¬ 
placement  are  both  significant.  These  results  are  exhibit¬ 
ed  in  Fig.  6.  In  the  case  of  propagation  of  an  initially 
Gaussian  transverse  wave  form  incident  in  a  homogene¬ 
ous  plasma,  with  the  parameters  a.  =0.248  jtm. 
/„  =  -/,  =  3  X  10|J  Wcm':,  r0  =  3/um,  and  ,V,0  =  7.5 
X  10:o  cm  ~ 3  (a,  =248.6192,  a,  =  -=-,  F0/7>cr  =  22.252), 
the  numerical  computations  show  that,  as  soon  as  the 
first  focus  on  the  axis  of  propagation  is  formed,  electronic 
cavitation  occurs.  Specifically,  this  leads  to  complete  ex¬ 
pulsion  of  the  electronic  component  of  plasma  from  the 
paraxial  domain  [Fig.  6(b)].  This  process  results  in  a 
quasistabilized  cavitated  channel  in  the  electron  distribu¬ 
tion  which  extends  along  the  entire  axis  of  propagation 
past  the  location  of  the  initial  focus  [26].  We  note  that 
some  of  the  lowest  stationary  solutions  corresponding  to 
the  relativistic  and  charge-displacement  problem,  includ¬ 
ing  cavitation,  were  developed  by  Sun  et  al.  [18].  The 
analysis  shows  that  the  first  focus  involves  about  45%  of 
the  total  incident  power  of  the  propagating  energy.  A 
fraction  of  the  remaining  power  is  dissipated  through 
diffraction  on  the  periphery,  while  another  component  is 
temporarily  involved  in  the  formation  of  a  pulsing  ring- 
shaped  structure  [Fig.  6(a)].  Subsequent  energy  exchange 
between  the  ringlike  feature  and  the  paraxial  zone  is  ob¬ 
served,  and  as  a  consequence  of  this  interaction,  a  certain 
part  of  the  energy  of  the  pulsing  ring  is  diffracted  away 
while  the  remaining  power  joins  the  paraxial  domain. 

In  finer  detail,  the  following  aspects  of  the  evolution  of 


the  radiative  energy  are  also  revealed  by  the  calculations. 
As  described  above,  after  the  formation  of  the  first  focus, 
considerable  power  is  transferred  from  the  paraxial  focal 
zone  to  the  ring-shaped  feature,  which,  at  this  stage  in 
the  evolution  of  the  pulse,  contains  approximately  68% 
of  its  total  initial  power.  This  intense  ring,  which  spreads 
away  from  the  paraxial  domain,  produces  a  correspond¬ 
ing  ring-shaped  cavity  in  the  electronic  distribution  [see 
Fig.  6(b)].  The  refraction  resulting  from  this  strongly 
perturbed  electron-density  profile,  together  with  the  rela¬ 
tivistic  self-focusing  mechanism,  causes  the  wave  to  re¬ 
turn  energy  to  the  core  of  the  beam.  Thus,  the  charge 
displacement  produces  a  potent  additional  self-focusing 
action,  which  leads  to  the  formation  of  a  confined  paraxi¬ 
al  mode  of  high  intensity  stabilized  along  the  axis  of 
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FIG.  6.  The  self-channeling  of  a  pulse  with  a  Gaussian  initial 
transverse-intensity  distribution  (.V,=2  in  Eq.  (27)]  and  a  rial 
incident  wave  front  in  initially  homogeneous  plasma  in  the  ca^e 
of  the  relativistic  and  charge-displacement  propagation  with 
/,  =  3X  lO'”  W/cnr.  r„  =  3  fim,  A.=0.248  jim.  arc 
V,  =.V  ,  =  7.5X  lCF0  cm'1,  (a)  The  distribution  of  the  norma 
ized  intensity,  (bi  The  distnbution  of  the  normalized  eiectr 
density,  >ci  Radial  dependence  of  the  asymptotic  solutions  ’ 
the  normalized  amplitude  [ /, ( r ) //0 ] 1  and  the  normalized  e:e. 
iron  density  ,Vf(ri/.V0  for  s  =0.554. 
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propagation.  This  phenomenon  of  confined  propagation 
is  designated  as  self-channeling.  It  should  be  noted  that 
the  calculations  of  the  electrostatic  energy  associated 
with  the  charge  displacement,  which  is  given  by  the  ex¬ 
pression  =  1 J0  E2r  dr,  with  the  electrostatic  field 

E.  defined  by  the  equation  VE.  =  —4 rrp  and  r.  designat¬ 
ing  the  radial  extent  of  the  channel,  show  that  this  energy 
can  be  relatively  small.  Specifically,  for  the  conditions 
represented  in  Fig.  cbi  'c  =  45.4  ^m],  the  electrostatic 
energy  ( IF)  accounts  for  only  0.18^  of  the  totai  energy 
of  the  laser  radiation  per  unit  length. 

The  essential  finding  of  these  calculations  is  that  the 
combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms  produces  a  strong  tendency  for 
the  generation  of  spatially  highly  confined  modes  of  prop¬ 
agation  which  are  stabilized  along  the  axis  of  propaga¬ 
tion.  Furthermore,  the  study  of  a  range  of  other  cases  in¬ 
dicates  that  these  modes  are  exceptionally  stable  and  that 
a  considerable  fraction  of  the  incident  power  can  be 
confined  in  them.  The  result  is  the  controlled  generation 
of  a  very  high  peak  intensity  in  these  channeled  modes 
with  values  reaching  —  10:i  W/cm;  for  the  range  of  con¬ 
ditions  studied. 

Important  characteristics  of  the  asymptotic  behavior 
of  these  confined  modes  have  also  been  established.  It 
has  been  shown  that  the  distribution  of  the  amplitude 
u  i  r,z  1,  for  large  values  of  z,  tends  asymptotically  to  the 
lowest  eigenmode  of  the  nonlinear  Schrodinger  equation 
see  Sec.  III).  Specifically,  for  the  example  discussed 
above,  the  computations  have  demonstrated  that  the 
asymptotic  radial  amplitude  distribution  corresponds  to 
the  lowest  eigenmode  with  s  ==0.554.  In  this  case,  the 
asymptotic  intensity  distribution  /t(  ri  =  Cr,:0(  n  contains 
46c‘c  of  the  total  incident  power. 

The  normalized  asymptotic  field  amplitude 
[/,.)r 2  and  the  corresponding  normalized  asymp¬ 
totic  plasma  electron  density  .V((n/.V0  are  depicted  in 
Fig.  6(c).  It  should  be  noted  that  the  profiles  of  the  inten¬ 
sity  /<r,z>  and  the  electron  density  .V(r,z),  obtained  as 
the  results  of  the  dynamical  calculations  of  the  propaga¬ 
tion  for  z=400  jim.  differ  from  /,(/•)  and  .V,i r)  for 
5  =0.  554  by  much  less  than  1 CT.  In  addition,  we  observe 
that  the  energy  of  charge  displacement  for  rc  =2.5  /im  in 
the  asymptotic  state  accounts  for  a  fraction  of 
4.45X10  4  of  the  total  energy  of  the  beam  per  unit 
length.  We  observe  that  this  tendency  for  a  solution  of  a 
nonlinear  Schrodinger  equation  involving  a  saturable 
nonlinearity  to  converge  to  the  lowest  stationary  solution 
was  originally  discovered  by  Zakharov,  Sobolev,  and 
Synakh  [37], 

For  the  range  of  parameters  studied,  the  calculations 
clearly  show  that  the  charge  displacement  has  a  very 
strong  influence  on  the  character  of  the  propagation  after 
the  first  focus  is  formed.  The  pulsing  intensity  structure, 
consisting  of  alternating  foci  on  the  axis  of  propagation 
and  peripheral  focal  rings,  which  is  the  usual  behavior  for 
the  purely  relativistic  seif-focusing  [see  Figs.  5<a)  and 
5'bi],  is  converted  into  a  stabilized  and  uniform  channel. 
Collaterally,  a  stabilized  cavitated  channel  in  the  electron 
density  is  also  formed.  We  remark  that  the  periodic  in¬ 


tensity  structure  characteristic  of  the  relativistic  self- 
focusing  occurs  for  great  values  of  z  and,  therefore,  may 
be  regarded  as  the  corresponding  asymptotic  solution  of 
the  purely  relativistic  case.  The  charge  displacement 
leads  to  the  formation  of  the  asymptotic  amplitude  distri¬ 
bution  represented  by  the  corresponding  lowest  eigen¬ 
mode  U,  d(  r)  instead  of  this  pulsing  structure. 

Computations  have  been  also  performed  for  the  propa¬ 
gation  of  incident  plateaulike  wave  forms  with  flat  in¬ 
cident  phase  fronts  in  both  homogeneous  plasmas  and 
plasma  columns.  Additional  calculations,  have  also  ex¬ 
amined  the  behavior  of  focused  Gaussian  and  plateaulike 
incident  wave  forms  as  well  as  defocused  Gaussian  in¬ 
cident  wave  forms  in  homogeneous  plasmas  see  Figs. 
7-11).  It  is  found  that  the  main  features  described 
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FIG.  7.  The  self-channeling  of  a  pulse  with  a  hyper -Gaussian 
incident-transverse-intensity  distribution  [ .V,  =  8  in  Eq  2* 
and  a  flat  initial  wave  front  in  initially  homogeneous  plasma  ■  r 
the  case  of  a  relativistic  and  charge-displacement  propagaii' 
with  f0  =  3X10i<>  W/cm;,  rQ  =  3  jam,  >,  =  0.248  qm.  rJ 
V,  =  .V. .,  =  7. 5 x  10;o  cm'1.  ia)  The  distribution  of  the  norma, 
ized  intensity,  ib)  The  distribution  of  the  normalized  elect  - 
density,  c)  Radial  dependence  of  the  asymptotic  solutions 
the  normalized  amplitude  [/,( r) //0 ]*  :  and  the  normalized  -ce. 
iron  density  .V.lrl/.V,  for  s  =0.515. 
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above,  namely  (1)  self-channeling,  (2)  stabilization  of  the 
mode  of  propagation,  (3)  the  confinement  of  a  substantial 
fraction  of  the  incident  power,  and  (4)  the  formation  of 
paraxial  cavitated  channels  in  the  electron  distnbution 
are  common  aspects  of  the  dynamics  over  a  wide  range  of 
conditions. 

The  character  of  these  findings  is  now  illustrated  with 
five  representative  examples.  For  the  propagation  of  a 
beam  with  an  incident  plateaulike  transverse-intensity 
profile,  given  by  Eq.  (27)  with  .V,  =  8,  which  has  a  flat  ini¬ 
tial  wave  front  incident  on  an  initially  homogeneous  plas¬ 
ma  (the  values  of  the  parameters  have  been  adopted  from 
above  and  P0//Vr  =  20. 168),  the  asymptotic  radial  profile 
of  the  intensity  distribution  contains  approximately  77% 
of  the  incident  power  [Fig.  7(a)].  The  corresponding 
asymptotic  amplitude  is  the  lowest  eigenmode  of  the  non¬ 
linear  Schrodinger  equation  with  s  ^0.515.  The  propa¬ 
gation  of  the  same  wave  form  in  a  plasma  column,  with 
the  initial  electron  distribution  defined  by  Eq.  (29)  with 
.V3=8  and  r,  =  r0,  results  in  a  quasistabilized  intensity 
distribution  containing  34%  of  the  incident  power  [Fig. 
8(a)],  The  purely  relativistic  propagation  of  incident  pla¬ 
teaulike  pulses  [  .V ,  =  8  in  Eq.  ( 27)]  in  plasma  columns 
[.V3  =8  and  r.  =r0  in  Eq.  (29)]  also  results  in  the  forma¬ 
tion  of  quasistabilized  regimes  which,  in  this  case,  arise 
dynamically  from  the  defocusing  action  of  the  refraction 


FIG  8.  The  self-channeling  of  a  pulse  with  a  hyper-Gaussian 
incident-transverse-intensity  distribution  [,V.=8  in  Eq.  '27'] 
and  a  fiat  initial  wave  front  in  a  preformed  column-shaped  plas¬ 
ma  [.V,  =8.  r.  =  r  ,  in  Eq.  '29i]  for  the  case  of  relativistic  and 
charge-displacement  propagation  with  /,  =  3X101'*  W/cm,- 
r,  —  3  qm.  A=0.248  jam,  and  .V0  =  .V,  ,  =  7.5X  10‘°  cm  *',  'a) 
The  distnbution  of  the  normalized  intensity.  >b»  The  distribu¬ 
tion  of  the  normalized  electron  density. 


generated  by  the  transverse  profiles  of  the  electron  densi¬ 
ty  [24],  For  this  situation,  approximately  25%  of  the  to¬ 
tal  incident  power  is  confined  [Fig.  5(d)].  The  compar¬ 
ison  of  these  cases  with  the  two  examples  discussed  above 
involving  the  charge-displacement  mechanism  indicates 
that  the  increase  in  the  value  of  the  confined  power  stems 
principally  from  the  substantial  additional  focusing  ac¬ 
tion  arising  from  the  inhomogeneous  electron  distribu¬ 
tion  produced  by  the  ponderomotive  force. 

Overall,  the  computations  reveal  that  the  charge  dis¬ 
placement,  which  generally  results  in  electronic  cavita¬ 
tion,  plays  an  important  role  in  stabilizing  the  mode  of 
confined  propagation  that  develops  dynamically.  This 
stabilization  naturally  occurs  by  refraction  of  the  radia- 
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FIG.  9,  Self-channeling  of  an  initially  focused  [R,  =  R..y  2 
m  Eq  1 3 5 >]  pulse  with  a  Gaussian  incident-transverse-intensity 
distribution  [.V,  =2  in  Eq.  (35)]  in  initially  homogeneous  plasma 
for  the  case  of  relativistic  and  charge-displacement  propagation 
with  /1=3X101’  W/cm:.  r0  =  3  qm,  A.  =  0.248  ^m.  and 
V  =  S',  ,  =  7  5  x  10;o  cm  ~  (a)  The  distnbution  of  the  normal¬ 

ized  intensity,  ibl  The  distribution  of  the  normalized  electron 
density  o  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  [ /,( n  //„  ] 1  :  and  the  normalized  elec¬ 
tron  density  .V,<r »/.V0  for  s  =0.566. 

124 


5842 


A.  B.  BORISOV  et  al. 


45 


tion  into  the  central  paraxial  region.  Indeed,  the 
influence  of  the  charge  displacement  on  the  propagation 
is  so  strong  that  self-channeling  occurs  even  in  the  cases 
of  extremely  focused  and  defocused  incident  wave  forms. 
The  results  of  the  corresponding  calculations  for  ex¬ 
tremely  focused  inctdent  wave  forms,  for  Rf—R.^/2 
Eq.  -35'j,  with  both  Gaussian  \V;=2  in  Eq.  1 3 5 » ]  and 
plateaulike  ‘.V-=8  in  Eq.  '35']  incident-transverse- 
mtensitv  distributions,  are  depicted  in  Figs.  9  and  10.  re¬ 
spectively  'the  parameters  of  beam  and  plasma  being  the 
same  as  the  other  examples  above'.  In  these  two  cases, 
the  asymptotic  transverse  profiles  of  the  amplitude  and 
the  electron  density  are  found  to  correspond  to  the  lowest 
eigenmodes  of  the  nonlinear  Schrodinger  equation  with 
values  of  the  parameter  s  =0.566  'Fig.  9)  and  s  =0.505 
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FIG  10  The  self-channeling  of  an  initially  focused 
R .  =R..,/2  in  Eq.  '35']  pulse  with  a  hyper-Gaussian  tncident- 
transverse-mtensity  distribution  [.V,  =8  in  Eq.  1 35’]  in  initially 
homogeneous  plasma  for  the  case  of  relativistic  and  charge- 
displacement  propagation  with  /.,  =  3x  10'J  W/cm\  r,  =  3  urn. 
>.  =0  248  wm.  and  .V,  =  V  ,  =  7  5  x  10;°  cm  .  ai  The  distribu¬ 
tion  of  the  normalized  intensity  b'  The  distribution  of  the  nor¬ 
malized  electron  density,  d  Radial  dependence  of  the  asymp¬ 
totic  solutions  for  the  normalized  amplitude  [/  in//,]1  -  and 
the  normalized  electron  density  .V.i  n  .V,  for  s  =0. 505. 


(Fig.  10).  Note,  in  strong  contrast  to  the  situation  involv¬ 
ing  charge  displacement,  the  analogously  strongly  fo¬ 
cused  incident  wave  forms  for  the  purely  relativistic  case 
propagate  in  the  single-focus  regime  (Fig.  5(c)],  and  stable 
confinement  does  not  develop. 

Figure  1 1  illustrates  the  relativistic  and  charge- 
displacement  propagation  of  an  initially  strongly  de¬ 
focused  wave  form  [R/=  —R.0/2  in  Eq.  (35)]  having  a 
Gaussian  initial  transverse-intensity  distribution  (.V,  =2 
in  Eq.  '  35)].  After  the  initial  stage  of  the  defocusmg,  this 
pulse  evolves  into  a  paraxial  structure  that  is  analogous 
to  those  described  above.  Furthermore,  in  this  example 
the  asymptotic  transverse  profiles  of  the  amplitude  and 
the  electron  density  are  found  to  correspond  to  the  lowest 
eigenmode  of  the  nonlinear  Schrodinger  equation  with 


(0 


I,  >■  /■  ' 

V,  r  V- 
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FIG.  11.  The  self-channeling  of  an  initially  defocused 
[R  .  -  -  Rr0/2  in  Eq.  i35)]  pulse  with  a  Gaussian  incident- 
transverse-intensity  distribution  [,V3=2  in  Eq.  t35)]  in  initially 
homogeneous  plasma  for  the  case  of  relativistic  and  charge- 
displacement  propagation  with  f0  =  3X1019  W/cm*.  r,  =  ?  urn. 
k  =0.248  /am.  and  .V0  =  .V,0  =  7.5X  10-’°  cm"’,  (a)  The  distribu¬ 
tion  of  the  normalized  intensity,  (b)  The  distribution  of  the  nor¬ 
malized  electron  density,  (c)  Radial  dependence  of  the  asymp¬ 
totic  solutions  for  the  normalized  amplitude  [/,(n//  '  -  ano 
the  normalized  electron  density  .V,(  r)/.V0  for  s  =0. 800 
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j  =0.800. 

For  these  values  of  the  dimensionless  parameters  as¬ 
sumed  in  these  examples  (Figs.  9-11)  (a,  =248.6192, 
a2  =  },  that  is  F0/Fcr  =  22.252  for  Gaussian  incident 
wave  form),  the  magnitude  of  the  quantity  Rf=Rf0/ 2, 
the  factor  determining  the  degree  of  the  initial  focusing 
(defocusing).  is  0.124.  For  A.  =0.248  /im  and  r0  =  3  pm, 
the  focusing  distance  R  =kr^Rf  0/ 2,  used  in  the  exam¬ 
ples  above,  has  a  magnitude  of  28.35  pm.  The  fact  that 
the  self-channeling  of  such  extremely  focused  or  de- 
focused  distributions  occurs  (Figs.  9-11)  means  the  self¬ 
channeling  should  be  expected  to  be  nearly  independent 
of  the  initial  conditions  of  focusing  or  defocusing. 

A  critical  power  P.t  for  self-channeling  arising  from 
the  relativistic  charge-displacement  mechanism  can  be 
defined.  We  now  discuss  this  issue  with  a  model 
developed  in  an  earlier  study.  The  threshold  power  Pct  is 
defined  as  the  power  that  separates  the  asymptotic  behav¬ 
ior,  with  respect  to  large  distance  of  propagation  (z),  into 
two  distinct  classes.  Since  its  definition  rests  on  an 
asymptotic  property,  it  can  only  exist  in  the  limiting  case 
of  vanishing  losses  (/z=0  in  Ref.  [24]).  For  a  power 
P  <  Pcr,  the  asymptotic  transverse-intensity  profile  tends 
to  zero  at  large  z.  In  contrast,  for  P  >  Pcr  the  asymptotic 
profile  of  the  intensity  tends  to  the  lowest  eigenmode  of 
the  governing  nonlinear  Schrodinger  equation.  Addition¬ 


ally,  the  normalized  value  of  this  critical  power  is  exactly 
twice  the  corresponding  normalized  value  of  the  critical 
power  for  the  Kerr  efiFect  self-focusing  found  in  cubic 
media  for  initially  flat  wave  forms  [34,37].  The  expres¬ 
sion  corresponding  to  this  critical  power  [36]  for  the  rela¬ 
tivistic  and  charge-displacement  mechanism  is 

Pcr={m;0ci/e2)  f*g'o(p)pdp(o/cop0)2 

=  1.6198X  10"Ww,,,0)2  W  ,  .36* 

where  g0(p )  is  the  Townes  mode  [34], 

It  has  also  been  shown  that  the  value  of  the  relativistic 
and  charge-displacement  self-channeling  threshold 
power,  in  cases  involving  initially  focused  or  defocused 
beams,  exceeds  Pcr  and  depends  on  the  degree  of  the  ini¬ 
tial  focusing  or  defocusing.  For  a  given  magnitude  of  the 
curvature  of  the  wave  front,  the  value  of  the  threshold 
power  for  initially  defocused  beams  is  greater  than  that 
for  initially  focused  beams.  Finally,  it  has  been  shown 
both  analytically  and  numerically  that  for  self-channeling 
of  an  arbitrary  wave  form  to  occur,  it  is  sufficient  that  the 
Hamiltonian  of  the  purely  relativistic  case,  considered  as 
a  functional  of  the  initial  transverse  amplitude  distribu¬ 
tion,  should  be  negative  [36].  The  precise  statement  for 
this  condition  is 


P:\u0 


u0  2  —  <  2/a; )[(  1  -t-a2!u0:)l/2  -  1  ]) 


r  dr  <  0  . 


We  now  comment  on  the  behavior  of  pulses  having  ini¬ 
tial  amplitude  distributions  C/0(r)  close  to  higher  eigen- 
modes  FI  )f(r>,  n  >  1.  These  higher  modes  generally  are 
associated  with  electronically  cavitated  channels.  The 
self-channeling  in  these  cases  could  result  in  asymptotic 
distributions  of  intensity  and  electron  density  corre¬ 
sponding  to  certain  higher  eigenfunctions 
i  =  F*„tr),  .V.  „i  r),  nil.  In  particular,  this  conjec¬ 
ture  has  established  an  initial  amplitude  distribution 
C'0ir i,  which  is  close  to  the  first  eigenmode  F,j(r)  with 
s  =0. 544.  Direct  numerical  calculations  showed  that  the 
asymptotic  distributions  of  both  intensity  and  electron 
density  that  evolved  in  this  case  correspond  to  the  first 
eigenmode  /5 ,( n  =  V’pr).  ,Vt  l(ri  with  s  essentially  equal 
to  0.544.  It  should  be  additionally  noted,  however,  that 
these  higher  eigenmodes  In  i  1  >  are  quite  possibly  unsta¬ 
ble  against  small  azimuthal  perturbations  that  destroy 
their  axial  symmetry. 

The  principal  result  of  this  section  is  the  finding  that 
the  combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms  can  result  in  self-channeling 
with  the  formation  of  stabilized  paraxial  modes  over  a 
rather  wide  range  of  physical  conditions.  Moreover, 
these  spatially  confined  modes  are  generally  associated 
with  corresponding  cavitated  channels  in  the  electron 
density  Finally,  the  characteristics  of  these  channeled 


modes  have  asymptotic  behavior  that  is  described  by  the 
appropriate  lowest  eigenmodes  of  the  governing  non¬ 
linear  Schrodinger  equation. 

V.  CONCLUSIONS 

A  theoretical  approach  suitable  for  the  numerical  in¬ 
vestigation  of  the  two-dimensional  (r,z)  dynamics  of 
propagation  of  coherent  ultrashort  (r,  »r»rt)  relativ  - 
istic  laser  pulses  in  cold  underdense  plasmas  has  been 
developed.  Four  basic  physical  phenomena  are  included 
within  the  scope  of  this  method.  They  are  (l)  the  non 
linear  dependence  of  the  index  of  refraction  due  to  the 
relativistic  increase  in  the  mass  of  the  electrons.  (21  the 
variation  of  the  index  of  refraction  resulting  from  the 
perturbation  of  the  electron  density  by  the  ponderom. 
tive  force,  (3)  the  diffraction  of  the  radiation,  and  ( 4 >  the 
refraction  caused  by  nascent  transversely  inhomogener  .- 
electron-density  distributions.  The  equations  studied 
this  work  may  be  regarded  as  the  generalization  of  t h.  - 
treated  previously  in  other  studies  [18,19,33]  involving 
itially  inhomogeneous  plasmas.  Further  studies  are  c 
tinuing  with  an  extension  of  this  analysis,  which  incluo 
the  azimuthal  coordinate  in  the  description  of  the  pn  r 
gation. 
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The  main  conclusions  stemming  from  the  calculations 
are  the  following. 

(1)  The  cooperative  effect  of  the  relativistic  and 
charge-displacement  mechanisms  leads  asymptotically  to 
stable  high-intensity  z-independent  modes  of  self¬ 
channeling.  and  a  major  fraction  of  the  incident  power 
can  be  confined  in  these  paraxial  modes.  Stable  cavita¬ 
tion  of  the  electron  density  is  a  general  feature  of  these 
spatially  confined  modes. 

(2)  The  c-independent  modes,  to  which  the  solutions  of 
the  equation  describing  the  relativistic  and  charge- 
displacement  propagation  tend  asymptotically,  are  recog¬ 
nized  as  the  lowest  eigenmodes  of  the  governing  nonlinear 
Schrodinger  equation. 

1 3)  A  separate  study  of  purely  relativistic  propagation 
shows  that  beams  with  flat  incident  phase  fronts  exhibit  a 
pulsing  behavior  in  homogeneous  plasmas  but  can  under¬ 
go  quasistabilization  in  suitably  configured  plasma 
columns.  However,  sufficiently  sharply  initially  focused 


beams  generally  exhibit  the  development  of  only  a  single 
focus.  In  significant  contrast,  the  present  study  shows 
that  with  both  the  relativistic  and  charge-displacement 
mechanisms,  initially  focused  beams  also  generally  lead 
to  confined  modes  of  propagation.  Finally,  it  should  be 
noted  that  the  equation  for  relativistic  self-focusing  can 
be  considered  as  a  general  model  equation  of  propagation 
in  saturable  nonlinear  media. 
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Abstract 


The  spatial  control  of  the  energy  deposited  for  excitation  of  an  x-ray  amplifier  plays  an  important 
role  in  the  fundamental  scaling  relationship  between  the  required  energy,  the  gain  and  the  wavelength 
New  results  concerning  the  ability  to  establish  confined  modes  of  propagation  of  short  pulse 
radiation  of  sufficiently  high  intensity  in  plasmas  lead  to  a  sharply  reduced  need  for  the  total  energy 
deposited,  since  the  concentration  of  deposited  power  can  be  very  efficiently  organized. 


I.  Discussion  of  Research 

Recent  theoretical1'2  and  experimental3  studies  have  led  to  a  fundamental  development  concerning 
the  generation  of  high-brightness  x-ray  sources.  These  results  affect  our  ability  to  controllabiy 
apply  very  high  power  densities  in  materials,  the  basic  issue  for  the  creation  of  bright  and  efficient 
sources  of  radiation  in  the  x-ray  range.  The  main  significance  of  this  work  is  the  establishment 
of  the  scaling  law  concerning  the  energy  requirements  for  x-ray  amplification  in  the  kilovolt  range 
shown  m  Fig.  1.  Importantly,  the  parameters  represented  in  Fig.  1,  and  which  define  the  relationship 
presented,  are  based  on  both  theoretical  and  experimental  information. 

The  critical  governing  issue,  which  determines  the  scaling  relationship  between  the  required 
excitation  energy  (E)  and  the  amplifier  gain  (G)  of  x-ray  lasers,  is  the  spatial  control  of  the 
deposited  energy.  The  information  presented  in  Fig.  2  shows  that  optimizing  the  gain  (G)  per 
unit  energy  (E)  calls  for  the  guided  mode  of  propagation  in  order  to  optimally  control  the  deposition 
of  the  energy.4  Overall,  in  comparison  to  traditional  forms  of  excitation,  for  a  fixed  x-ray  energy 

output  (Ex)  and  wavelength  (Xx),  a  reduction  of  several  orders  of  magnitude  in  the  necessary 
energy  (E)  results,  as  shown  in  Fig.  1,  if  this  form  of  confined  (channeled)  propagation  can  be 
achieved.  Therefore,  if  this  scaling  holds,  a  relatively  small  and  useful  laboratory-scale  technology 
becomes  feasible. 

Recent  experiments,3  which  are  supported  by  carefully  developed  theoretical  analysis,1-2  have 
demonstrated  the  basic  physics  of  a  new  form  propagation  exactly  of  the  type  necessary  for  the 
implementation  of  x-ray  lasers  of  a  fundamentally  new  regime  of  electromagnetic  propagation  is 
expected  to  arise  in  plasmas  for  short-pulse  radiation  at  sufficiently  high  intensity.  Dynamical 
calculations  of  the  propagation  in  plasmas,  incorporating  both  relativistic1'®  and  charge-displacement 
mechanisms,2'®-^  indicate  that  the  combined  action  of  these  processes  can  lead  to  a  new  stable 
form  of  spatially  channeled  propagation.  Specifically,  these  experimental  studies  which  have 
examined  a  new  relativistic  regime  of  high-intensity  short-pulse  propagation  in  plasmas,  present 

evidence  for  the  information  of  such  a  stable  mode  of  spatially  confined  (channeled)  propagation 
For  an  electron  density  of  -  1.35  x  103‘  cm°  and  a  power  of  ~  3  x  1QU  W,  the  results  indicate 

a  i~hannel  radius  <  1  pm  and  a  peak  intensity  -  1018  W/cm3.  Comparison  of  these  findings  with 
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Energy  (E)/Wavelength  (X)  Scaling 


Quantum  Energy  (fiojx)  keV 


Fig  i  Scaling  relationship  between  required  excitation  energy  (E)  and  quantum  energy 
characteristic  of  the  amplifier.  Parameters:  total  gain  exponent  G  =  100,  energy  effic 
nx  =  10'3,  channel  diameter  <5  =  3  pm,  x-ray  (hux)  cross  section  for  stimulated  em.s 
ox,  x-ray  cross  section  for  stimulated  emission  for  radiatively-broadened  transition  ayp 
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X-Ray  Laser  Scaling 

Spatia!  Distribution/ Amplifying  Volume 
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G  -  p*crx  t 


Laboratory  Scale  Technology— 


small  6 


E  Xcjx  5  2 
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propagation 


(52  /  A 
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/.  |  Guided  Mode  of  PropagatiorT"] 


Fig  2  Spatial  distribution  of  energy  of  excitation  (E)  for  an  x-ray  (fiwx)  amplifier.  Param- 
are  the  same  as  in  Fig.  1  with  X  the  wavelength  of  the  excitation  energy,  ass^ 
longitudinally  delivered,  and  with  p,  p*.  and  pe  representing  the  particle,  inversion 
electron  densities,  respectively.  The  analysis  shows  that  optimization  of  G/E  require 
guided  mode  of  propagation  so  that  high  concentrations  of  power  can  be  organized 
high-aspect-ratio  spatial  volumes. 
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a  dynar.iical  theory2  yield  close  agreement  for  both  the  longitudinal  structure  and  the  radial 
extent  of  the  propagation  observed.  These  results  represent  a  profound  change  in  the  field  of  x- 

ray  laser  research  because  they  alter  drastically  the  fundamental  scaling  relationships  among  the 
relevant  physical  variables. 

The  implications  of  this  development  for  general  applications  to  x-ray  imaging  and  the  micro¬ 

characterization  of  condensed  matter  are  extremely  important  and  propitious.  In  terms  of  the  x- 
ray  source,  they  are  (1)  that  a  properly  controlled  energy  deposition  rate,  sufficient  for  tne 
production  of  stimulated  x-ray  sources  up  to  a  few  kilovolts  in  quantum  energy,  can  now  ce 

achieved  with  an  excitation  energy  of  -  1  J,  (2)  that  an  x-ray  output  energy  of  ~  1  mj  per  puise 
is  achievable  with  laboratory-scale  technology,  and  (3)  that  an  x-ray  beam  diameter  (-2-3  _mi 
arises  as  a  natural  consequence  of  the  physics.  These  parameters  represent  an  exceptionally  high 
peak  brightness  figure  that  permits  a  new  and  completely  unexplored  range  of  physical  measurements 
to  be  made.  Indeed,  a  high-brightness  source  of  this  nature  is  ideal  for  the  microimaging  of 
condensed  matter  In  particular,  an  x-ray  source  with  these  parameters  is  perfectly  matched  to 
the  requirements  for  holographic  imaging  of  biological  materials11-1-14  in  terms  of  all  its  relevant 
properties,  specifically,  wavelength  (10  -  40  A),  pulse  energy  (-  1  mj),  pulse  length  (-  10'"3  s, 

beam  diameter  (-2-3  pm),  and  divergence  (-  1  mrad). 


II.  Conclusion 


Recent  experimental  and  theoretical  results  on  electromagnetic  propagation  at  high  intensities 
in  plasmas  lead  to  very  favorable  scaling  relationships  for  high-brightness  x-ray  amplifiers.  Such 
x-ray  sources  will  have  important  applications  in  the  holographic  imaging  of  biological  materials. 
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